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ABSTRACT 


A method of classification of digitized mul ti spectral 
images is developed and experimentally evaluated on actual 
earth resources data collected by aircraft and satellite. 
The method is designed to exploit the characteristic 
dependence between adjacent states of nature that is 
neglected by the more conventional simple-symmetric decision 
rule. Thus contextual information Is incorporated Into the 
classification scheme. The principle reason for doing this 
is to improve the' accuracy of the classification. . For 
general types of dependence this would generally require 
more computation per resolution element than the 
simple-symmetric classifier. But when the dependence occurs 
in the form of "redundance", the elements can be classified 
collectively, in groups, thereby reducing the number of 
classifications required. Thus a potential exists for 
increased, rather than’ decreased, efficiency. 

Basically, the method can be thought of as an image 
partitioning transformation that delineates (extracts) the 
statistically homogeneous groups (samples) of elements and a 



X 


sample classifier that classifies them. Various 
possibilities are considered for both operations. 

The main result is that a combination of the two is 

» 

found 'Which consistently provided the lov;est error fates, 

rivaling those obta'in'ed when' ground observational data was 

used to delineate the samples manually. The relative 

✓ 

efficiency of this method depends largely on the complexity 
of the classification, task. For relatively complex 

classification/ the time saved by sample classification more 
than compensates for the extra time required for 
partitioning. But for relatively simple classification the 
simple-symmetric classifier is faster. Of course in the 
latter case/ efficiency is not as great a consideration 
since the total CPU time involved is much less than in the 


former case. 
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CHAPTER 1 
INTRODUCTION 

The general objective of this thesis is to advance the 
state of the art of pattern recognition as it is applied in 

I 

remote sensing technology. This chapter opens with a 
discussion of pattern recognition and remote sensing systems 
that leads up to the specific problem under investigation. 
In the process much of the prevalent terminology is 

introduced. Other work that Is related to this problem is 
discussed in Section 1.3. 

1.1 Pattern Recognition Systems 

Man's most abundant source o.f information about a scene 
is the radiant electromagnetic energy which emanates from 
it. The information is embodied in the spatial, spectral, 
and temporal variations (patterns) of the radiance. The 
general process of extracting information from patterns 
(radiance or otherwise) is known as pattern recognition. 
The most common form of pattern recognition is 

"classification", the assignment of an observed pattern to 
one of several prespecified categories (classes). This 
requires a certain degree of experience; i .e. the 
recognition system must know the possible classes and have 
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some sort of unique characterization for each one. 
Typically this experience is "learned” from representative 
"training" patterns (or sets of patterns) that are supplied 
as references for each class. In the simplest case, each 
set of patterns Is a complete characterization of the class 
it represents. Then classification Is a straightforward 
matter of comparison. More generally, . a statistical 
characterization might be the only adequate approach, and 
the training patterns might be used to estimate statistical 
quantities. Classification then becomes a problem in 
statistical decision theory. i 

Of course it Is not always possible to prespecify the 
categories that a pattern might belong to. This is often 
true in scene analysis, where the number of possibilities 
can be enormous. Then pattern recognition can take the form 
of "description". In general, pattern recognition can 
involve both classification and description. A’ complex 

scene composed of rel-atively simple objects is .often 

♦ 

described by classifying the objects and recording their 
re-lative positions and orientations - in the scene. , This 
description might be considered the final result, or It 
might in turn be used to classify the scene itself. 

All systems that extract information from a scene 
consist of. a data collection system and a data processor. 
The purpose of the data collection system Is to reduce the 
scene to a manageable number of measurements (features) 
vdthout losing the desired information. Further reduction 
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(feature selection) is often possible In the processor. The 
choice of features obviously depends upon the information 
that is desired^ and conversely the information that can be 
extracted depends upon the choice of features. Most 
collection systems are similar in many ways to the human 
eye, which forms features by "samplirig"' the spatial, 
spectral, and temporal dimensions, thereby converting a 
scene into series of electrical pulses. Spatial sampling, 
can be accomplished by formi’ng an image of the scene on an 
array of detectors (electrical or chemical) or by scanning 
the image with an electrical detector. The resolution 

element of such a system is the projection of the detector 
back through the optical system onto the scene. It is 
commonly called a "pixel", short for picture element. The 
overall system resolution depends on both the pixel size and 
the interval between samples, which are normally about 

equal . 

Spectral sampling is accomplished by measuring the 
radiance of each resolution element with detectors 
(channels) that are sensitive to different spectral bands. 
A prism, grating, or interference filter Is often used to 
separate' the radiant energy spectrally before detection. 
Temporal sampling i s accompl i shed merel y- by taking spatial 
and spectral samples at discrete times. 

Depending on the type of information that is desired, 
one can emphasize or de-emphasize a particular dimension by 
sampling it relatively many or relatively few times. 


A 
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single black-and-white photograph, for example, emphasizes 
spatial Information since it Is created by sampling only 
once spectrally and once temporally, A color photograph 
contains three spectral samples and thus emphasizes both 
spatial and spectral Information. The extent to which a 
pattern is sampled falls under the category of "measurement 
complexity". Under-sampling results In loss of information, 
but over-sampling results In an excess of data to process. 
Technically, the data dimensionally increases faster: than 
its intrinsic dimensionality. 

"Data dimensionality" refers to the dimension of the 
measurement or observ.ation space, in which a sampled pattern 
can be considered an observation of a mul t i -d imens i onal 
random variable. The probability density of this random 
variable Is a function of N variables (dimensions), where N 
is the number of measurements. The "intrinsic" 

dimensionality of a random variable (X) Is the minimum 
dimension that another random variable (Y) can have if X Is 


uniquely related to Y. 

Thus 

* a • 

It IS 

the minimum 

number of 

measurements 

that could 

be 

used 

to convey 

the same 

informat i on 

as X if 

the 

relationship vyere 

known . 


Over-sampling Increases the data dimensionality, but the 
individual measurements tend to be more highly correlated 
causing the information conveyed per measurement to 
decrease. 

The Information that can be extracted from an image is 
also limited by the sophistication' of the processor which 
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must handle the data. Just as the necessary measurement 
complexity depends on the information being sought, so does 
the method of processing. The human mind is often an 
extremely good processor, particularly when the information 
is of primarily a spatial nature. For this purpose the data 
is presented in visual image form, . which is known as an 
••image-oriented" processing system. By contrast. In a 
"numerically-oriented" system the decision-making element is 
a computer, and the visual image plays little or no part. 
Advantages of the computerized approach are its high load 
(volume) capacity, comparatively low cost under high load, 
and capacity to handle high measurement complexity. 

1.2 Remote Sensing of Earth's Resources 

An important subject before the engineering and 
scientific community at the present time is the processing 
of scenes which represent tracts of the earth's surface as 
viewed from above. A typical scene consists primarily of 
regular and/or irregular regions arranged in a patchwork 
manner and each containing one class of surface cover type. 
These homogeneous regions are the "objects" In the scene. A 
basic processing goal is to locate and classify the objects 
and produce a description of the scene in terms of tabulated 
results and/or a "type-map". As in other image processing 
applications, the locations and spatial features (e.g. size, 
shape, orientation) of objects are revealed by changes in 
average spectral properties that occur at boundaries. But 




unllKe most other appl Icat ions^ the spatial features of an 
object often have only a weak relationship to its class. 
Research has shown, however, that many classes can be 
distinguished reasonably well on the basis of their spectral 
features, using statistical pattern classification 
techniques. Current research is directed toward use of 
temporal features as well, but not fn this investigation. 

Our Interest is in . the numerically-oriented system 
approach to processing these scenes. The input to the 
system is in the form of digitized mul ti-spectraK scanner 
(MSS) data stored on magnetic tape. A typical 

multi-spectral scanner' samples the spectral .dimension and 
one spatial dimension. ' The second spatial dimension is 
provided by the motion of the platform which carries the 
scanner over the region of interest, generating a 
raster-type scan. The temporal dimension is provided by 
rescanning the region at different times. 

Computer cl ass I fii cat ion, of MSS data is typically done 
by applying a "simple symmetric" decision rule to each 
pixel. ■ This means that each pixel is classified 

Individually, on the basis of its spectral measurements 
alone. A basic premise of this technique is that the 
objects of interest are large compared to the size: of a 
pixel. Otherwise a large proportion of pixels woujd be 

composites of two or more classes, making statistical 
pattern classification unreliable; i.e. the prespecified 
categories woufd be inadequate to describe the actual states 
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of nature, (For later reference we shall call this "Premise 
A".) Since the sampling interval is usually comparable to 
the pixel size (to preserve system resolution), it follows 
that each object is represented by an array of pixels. This 
suggests a statistical dependence between consecutive states 
of nature, which the^simple symmetric classifier fails to 
exploit. To reflect this property, we shall refer to simple 
symmetric classification as "no-memory" classification. 

One method for dealing with dependent states is' to 
apply the principles of compound decision theory or 
sequential compound decision theory. Abend III points out 
that a sequential procedure can be implemented relatively 
efficiently when the states form a low-order Markov chain. 
However the prospect is considerably less attractive when 
they form a Markov mesh, which is a more suitable model for 
tv/o-d imens i ona 1 scenes. Furthermore, estimation of the 
state transition probabilities could be another significant 
obstacle to implementation of such a procedure, A short 
appendix on the compound decision approach is included in 
this thesis. 

The compound decision formulation is a powerful 
approach for handl ing very general types of dependence. 
This suggests that perhaps by tailoring an approach more 
directly to the problem at hand, one can obtain similar 
results with considerable simplification. A distinctive 
characteristic of the spatial dependence In MSS data is 
redundance; i.e, the probability of transition 
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from state i to state j is much greater if j = I than if jj^i/ 
because the sampl ing interval is small compared to the size 
of an object. This suggests the use of an "image 
partitioning" transformation to delineate the arrays of 
statistically similar pixels before classifying them. Since 
each homogeneous array represents a statistical "sample" (a 
set of observations from a common population)/ a "sample 
classifier" could then be used to classify the objects. In 
this v)av/ the classification of each pixel in the sample is 
a result of the spectral properties of its neighbors as well 
as its own. Thus its "context" in the scene is used to 
provide better cl ass i f i cat ion. The acronym ECHO (extraction 
and classification of homogeneous objects) designates this 
general approach., 

A characteristic of both no-memory and compound 
decision techniques is that the number of classifications 
which must be performed is much larger than the actual 
number of objects in the scene. VJhen each classification 
requires a large amount of computation/ even the no-meinory 
classifier can be relatively slow. An ECHO technique would 
substantially, reduce the number of classifications/ 
resulting In a potential increase in speed (decrease .in 
cost). Whether or not this potential is realized depends on 
the efficiency of the partitioning operation. 

The goal of the current investigation is to further the 
development of the ECHO concept. in particular/ various 
processing options are devised, implemented, and tested on a 
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wide variety of 
determi ne thei r 
using no-tnemory 
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data sets. Input parameters are varied to 
effect, and performance comparisons are made 
classification as a norm. 


beproducbili^J 
or'iginaIj page 


1.3 Related Work 

The recent literature contains numerous references to 
image partitioning algorithms. Robertson 12| divides them 
into two main categories. "Boundary seeking" algorithms 
character i St i cal 1 y attempt to exploit object contrast. 
These techniques include local gradient |3,4l, template 
matching 151/ two-dimensional function fitting 161, 
clustering 14,7.1, and gradients estimated from 
variable-sized neighborhoods i8|. Two of these have been 
implemented with digitized mul t i spectral Imagery. 

Anuta 141, investigated a multivariate extension of a 
two-dimensional gradient operator. The gradient operator of 
a unispectral Image maps each pixel into a number which 
reflects the average positive difference between that pixel 
and its neighbors. The multivariate operator sums these 
numbers over, all spectral features for each pixel. Since 
the differences are generally larger for boundary pixels 
than for non-boundary pixels, thresholding this sum (for 
each pixel) ’ at the "proper" level provides a boundary 
enhanced version of the original image. This technique is 
relatively fast, but it has several serious problems. 
First, It is .inherently noisy, which is typical of 
d i f f pr ent i at i on techniques. It is also very sensitive to 
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the threshold level used. Furthermore the boundaries 
derived by this technique often fail to close- upon 
themselves. For example, a boundary line may become 
discontinuous or fade out completely, leaving the objects 
ambiguously defined. In special cases v/here the object 
shape is restricted 13,91, the true boundaries can sometimes 
be deduced, but in general they cannot. This may not be a 
serious drawback for applications such as 'image 
registration, but closed boundaries are necessary for sample 
classification. This particular problem is common to all 
the boundary seeking algorithms mentioned above. 

• Wacker |7| developed an algorithm for MSS data'which 
performs a cluster analysis (unsupervised classification) of 
a small region of the image' and then scans the result for 
the oresence of a boundary. ‘ The estimated boundary 
structure for the entire image is obtained simply by taking 
the union of the boundaries found in all such regions.- This 
is a much more time-consuming process, but it is less noisy 
and less sensitive to input parameters. Of course it 
suffers from the same open boundary problem as the other 
boundary seeking algorithms. 

The other category of image partitioning algorithms can 
be called "object seeking" algorithms, , which 
character i s t i cal 1 y exploit the Internal regularity 
(homogeneity) of the objects. As the name implies, an 
object seeking algorithm always produces well-defined 
samples (and thus closed boundaries as well). There are two 
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opposite approaches to object seeking/ which we shall call 
conjunctive and disjunctive. A conjunctive algorithm begins 
with a very fine partition and simplifies it by 
progressively merging adjacent elements together that are 
found to' be similar according to certain statistical 
criteria 110,111. A disjunctive algorithm begins with a 
simple partition and subdivides it until each element 
satisfies a criterion of homogeneity. For example, 
Robertson's algorithm 12,121 Is based on the premise that if 
a region contains a boundary, splitting the region 
arbitrarily will usually produce two subregions with 

significantly different statistical characteristics. 

Early v;ork in the application of sample classification 
to MSS data v;as reported by Huang 1131. His method of 
"polling" requires classification of the individual pixels 
in the sample and is thus relatively inefficient. Vlacker 
and Landgrebe 1141 investigated the "minimum distance 
approach" using parametric and non-parametr i c methods. Both 
studies relied on manual definition of the object 

boundaries, based on actual surface (ground) observations, 
to locate the samples that were classified. 

We combined Rodd's conjunctive partitioning algorithm 
with a minimum distance sample classifier and observed an 
improvement in classification accuracy over conventional 
no-memory classification, but processing time was increased 
1151. Gupta and Wintz 1161 added a test of second order 
statistics to Rodd's first order test, but obtained 
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essentially the'same results as the first order test at 
greater cost In processing time, Robertson |2,12i 
implemented a disjunctive partitioning algorithm with the 

j 

same minimum distance classifier. He obtained about the 
same classification accuracy as conventional nd-memory 
classification with an order of magnitude increase in 
processing time. This points to one essential difference 
between the disjunctive and conjunctive approaches. With a 
disjunctive approach/ every time a region is divided new 
sample statistics must be calculated from raw data. With a 
conjunctive approach, every time' two regions are merged the 
statistics for the resultant region can be obtained merely 
by "pooling" the statistics of the original two subregions. 
This results -in a significant computational advantage for 
the conjunctive approach. 

The current investigation is devoted to further 
development of the conjunctive approach. A much faster 
sampl e- cl ass i f ier Is proposed and tested. This problem Is 
discussed in Chapter 2. New statistical criteria are 
proposed as' wel 1 as new object seeking logic in Chapter 3, 
Extensive test results appear in Chapter 4, comparing 
different algorithms against' each other and against 
conventional no-memory classification. The main result is 
that the stability, cl ass I f i cat i on accuracy, and speed of 
the ECHO technique have been greatly improved. Compared to 
the no-memory classifier, consistently lower error rates are 
observed using an ECHO approach, and for a reasonably 
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complex classification its efficiency exceeds that of the 
conventional method. 
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CHAPTER 2 
CLASSIFICATION 

The motivation for object extraction is to enable 
faster and more accurate classification of the pixels within 
the object. In Section 2.3 we discuss the classification 
algorithms that accomplish this. They are based on a 
certain model of the objects to be classified, which is 
described next. 

2.1 Statistical Model of Mul t i -Spectral Scanner Data 

As we have indicated, a typical scene consists 
primarily of objects whose boundaries form a partition of 
the scene. The partition is generally unknown at the 
outset, but we can at least assume that it is relatively 
coarse compared to the size of a pixel. Each object in the 
scene belongs to some class. For representation purposes, 
each class is divided Into one or more "subclasses". They 
are also called "spectral classes" (as opposed to 
"informational classes") to indicate that they can be 
distinguished spectrally although it may not be useful to do 
so. Let W|j denote the jJLb. subclass of the iJLh class. Let 
F denote an object (represented by an array of pixels), and 
let X denote a pixel in some object. (The underbar is used 
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to indicate a q-d imens ional variable v/here q 

henceforth denotes the number of spectral channels.) Then 
FeWjj denotes the event that F belongs to the subclass 
The a-priori probability of this event is denoted by 
PCFeWjj ). In accordance with Section 1.2, we ignore any 
statistical dependence of this event on the spatial features 
of F. If there were a strong, known dependence then it 
could be used to ^help cla.ssify F, but that is not our 
intention, A consequence of this assumption is that 
PCXeV/..) = PCFeW;.)', and we denote both quantities simply by 

— I j .• I j ' 

P(Wjj). 

The pixels within a given ' object of a given spectral 

class are completely character i zed ' by their class- 

conditional, Joint, probability d i str i'but ion function. For 

no-memory classification, such a complete model is 

unnecessary; only the marginal distribution of each pixel Is 

required. Furthermore, the pixels within a single object 

are usually assumed to have a common (i.e. stationary) 

marg-inai distribution, which is due to the homogeneity of 
1 

the types of objects typically encountered in remote sensing 
aopl I cat ions . Although the data is digitized, it is 
convenient to represent this q-variate distribution by a 
continuous-parameter probability density function Cpdf) 
which, for subclass W.., Is denoted by pCX=x.!2 ^ W; • ) ot 
simply by p(it|W.j). -CThe vertical bar indicates conditional 
probab i 1 i ty ) . 
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Two pixels in spatial proximity to one-another are 
unconditionally correlated, with the degree of correlation 
decreasing as the distance between them increases. Much of 
this correlation is attributable to the effect. of dependent 
states, discussed in Section 1.2, which is the effect we 
wish to exploit. For simplicity we shall ignore other 
sources of correlation. Thus we assume that pixels within 
the same object are cl ass-condl tional 1 y independent; i.e. 

I 

each object Is a "simple" sample from one of the spectral 
class populations. Then the joint pdf of the pixels can be 
expressed as' just the product of their marginal pdf's. This 
approximation leads to fast, effect iVe (though subopt imal) 
processing algorithms, but theoretical predictions based on 
this simplified model should be interpreted cautiously. 
This aspect of modeling is discussed at greater depth in 
Appendix A. 


It is possible to express other 

characteristics In terms of the ones above. 

the ijJl class, then 

PCileW, ) = PCU XeW. . ) = Z P(W, , ) 

1 . ij : ‘J 


stat i St i cal 


If W| denotes 


2 . 1.1 


where U denotes the union of events. The pdf of X, 
conditional on this event. Is given by 

pC>llWi) = - 1- Z pCxIWt -)P(W; f 2.1.2 

PCWj) j 'J »J 

This equation defines the representation of a class in terms 
of Its subclasses. The unconditional pdf can be written in 



2.1 


17 


two vjays: 

P(x) = 22 P(x|W..)P(W..) = 2 PCxlWj)P(Wj) 2.1. '3 

' i j I 

Within this Framework, all that is required to complete 
the statistical model for a given scene (or class of scenes) 
is to specify the spectral classes that are present and 
assign an a-priori probability and conditional pdf to' each. 
Of course the- true distributions are assigned by nature, and 
the accuracy of the model depends on how wel 1 ■ we can 
estimate them. Fortunately ive are usually able to obtain 
estimates of the class-conditional pdf's based on training 
samples taken directly from the data ^set. For this we 
usually rely on actual surface (ground) observations or 
manual photo- interpretat i on to locate areas representing 

each class of cover-type. For the purpose of classifier 
design, we assume that the -size of each sample is 
sufficiently large that the error in the correspond i ng 
distribution estimate is negligible. The subject of 

training is discussed further in Chapter 4., c 

The distribution estimates can be parametric or 

non-parametr Ic in general. It has been found that the 
mul t i -var iate normal (MVW) distribution is a reasonable 
model for MSS data 117|; i ,e. p(ii|Wjj) = N(Mij/Cjj;ii^/' where 
N(M,C;x) = ('!2 ttC| exp( (iL-M) ' C \>i-M) ) ) “ ^ 2.1.4 

(Note that (x.-M) ' denotes the transpose of vector ( x,-M') . ) 
it follows that if AeWjj/ 

E(D == Mjj 



2.1 


18 


where E( • ) denotes statistical expectation. Thus and 
C. . are the mean vector and covariance matrix of the 
subclass distribution. Note that In order to obtain a 
parametric estimate of a MVN distribution, it is only 
necessary to estimate its first and second order moments. 
This is the approach that we will use. 

2.2 Mo-Memory Classification 

In order to introduce certain concepts that will be 
useful later, we now review some common techniques of 
no-memory classification including Cin one case) a 
discussion of a bound on the probability of error. 

2.2.1 Maximum A Posteriori Probability CMAP) Strategy 

Let X be a pixel, as before. Under the hypothesis that 
X eW,, the pdf of' A is p(X=xlXeW. ), which is given by 
equation 2.1.2. ■ Assuming that this function is accurately 
known, the hypothesis is ‘'simple". The goal of 
classification is to devise a strategy for choosing one of 
the possible classes (hypotheses) based on i, the observed 
value of I .e,. we must specify a function, W(x), which 
maps X into the set of possible classes. We can maximize 
the probability of a correct' deci s ion by always choosing the 
class, Wj, which has the jifjaximum a posteriori probability, 
P(XeW; 1 show this we merely write the probability of 
a correct decision in the following form: 



PCX e W(X)) =JpC>i e WCa) I l=x ) p(l=x) d^s. ’ 2.2.t.l 

j 

It Is apparent that this quantity Is maximized with respect 
to the decision function by adopting the .MAP decision rule. 
To implement this strategy we use the mixed form of Bayes 
rule to write 

p(X=x i X e W. )P(Wj ) 

P(1 e Wr l^=x.) = 2. 2. 1.2 

' p(X.-2<) 

The denomi nator .is independent of i, so we need only to seek 

the' i which maximizes the numerator. In other wordS/ for a 

given observation, WCx.) is chosen such that 

d()<=><U^ e VI( X.) ) P(W(i<) ) = max p(i<=>L!i( e Wj)P(Wj) 2. 2. 1.3 

i 

This result can also, be obtained' as a special case of Bayes 
decision rule for i.iinimum risk v/hen a "zero-one" loss 
function is assumed (i.e. when' the risk equals the 
probability of error). Thus’ it is often referred to as 
"Bayes classifier". 

2.2.2 Maximum Likelihood (ML) Strategy 

When all the classes are equ i probab 1 e, the MAP decision 
rule reduces to 

p(Jl=x.l)L ^ V7(j<) ) = max e W: ) 2. 2. 2.1 

!- 

t \ 

As a function of i, the statistic pC)(=iclZ ^ W; ^ is called 

» 

« 

the likelihood function, so this decision rule is called the 


maximum likelihood strategy 
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The ML strategy is usually a- reasonable approach even 
when the classes are’ not equ i probabl e. In particular, the 
'•;AP strategy tends to discriminate against classes whose 
a-priori probability is low; i.e, it encourages a relatively 
large conditional probability of error when a "rare" class 
occurs In order to minimize the overall error probability. 
Thus when one is interested in classifying the less abundant 
classes (as well as the more abundant classes) with 
reasonable accuracy, the MAP strategy may not be as 
desirable as one which makes more errors but distributes 
them more equitably among the classes. With the ML 
strategy, the conditional probability of error when the iili 
class occurs depends only on the degree of statistical 
"separability" (or "distance") between class i and the other 
classes. It is independent of the a-priorl probability of 
class i.. 

2.2.3 Generalized Maximum Likelihood (GML) Strategy 

Often the a-priori subclass probabilities are unknown. 
Then the hypothesis that e Wj is a composite hypothesis; 
i.e. p(j<lVi.) = DA|. p(xlW. .) where the coefficients are 

I J 

unknov;n. Of course^we know that A^* ^ 0 and “1* ^ 

procedure that has been found to be useful in this situation 
is to form maximum likelihood estimates of the unknown 
parameters under each hypothesis. Then the unknov^ns are 
replaced by their estimated values, and a hypothesis is 
selected by the ML strategy. We will refer to this 
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procedure as the sen^ral ized maximum likelihood strategy. 
The; resultant decision rule can be simply expressed in the 
fol lowing form: 

• p(illV(^)) = max max pC>LiWj.) 2. 2.3.1 

• • ^ J 

I J . 

where V(x) maps x into the set of spectral classes.' Then 
iKx.) is simply defined to be the informational class 
containing VC^i). 

We note that. the. GML strategy is equivalent to a ML 
strategy over the set of spectral, classes. Thus when all 
spectral classes are equ i probab 1 e/ It maximizes the 
probability of classify-Ing the observation into the correct 
one . 

2.2.4 Probability of Error For The GML Strategy 

Let Vj denote the iiii spectral clasS/ and let E be the 
event that X is classified into the wrong spectral class. 
Then 

PCE)' = Z'PCEU e Vj)P(Vj) 2 I 2 . 4 .I 

j 

If E|j is the event that Vj, produces a larger likelihood 

V 

statistic than Vj , then 

% 

P(E|Vj) = P( U EjjIVj) 2 PCEfjlVj) ■ 2.2. 4. 2 

I I 

i i 

Thus it is of some interest to investigate the pairwise 
error probabilities. 

. Let F|j,(T) = PCRijCX) >T!XeVj)’= P(L|j(X) >1 n(T) | X eVj.) , 
where Rjj(l) and L] j Cl) are the random variables: 



2.2.4 


22 


RjjCX) = 


pCXIVj ) 


= 1 I kel i hood ratio 


2. 2. 4.3 


pCXiVj) 

L;j(X) = In Rjj(X) 

Then P(E|-iV?) = F:;(l). Unfortunately the conditional 

distribution functions of and Ljj(^) are not usually 

explicitly available. But/ if we can find the moment 

generating function/ <{>{j(u)/ corresponding to the 
conditional distribution of- LjjC^) given X^Vj/ then we can 


bound F|jCT) as follows: 
FjjCT) j,(u)/ 0 s< u 


2. 2. 4. 4 


Furthermore 

F.j(T) v< :(u)/ u 1 

This is known as the Chernoff bound 118|. 

By def i n i t ion: 

(J)..Cu) = EUxp(uL..(X))|X e V.) 2. 2. 4. 5 

I j t j J 

When the subclasses are MVN/ the expectation can be 
explicitly evaluated 1181. The result is: 

4) J . ( u ) = 



Substituting into 2. 2.. 4. 4 provides the desired bound. in 
particular/ for u = .5 we have: 


Fij.(T) (|>.j(.5)//f 


2 . 2 . 4. 7 


V/e note in passing that -In (})|j(.5) is simply the Bhatta- 
charyya "distance" between subclasses Vj and Vj. 
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Combining equations 2. 2, 4,1/ 2.2,4, 2/ and 2, 2, 4, 7 gives 

a 

an expression For a bound on P(E); 

A 

PCE) Z PCV.) Z 2. 2. 4. 8 

j i ■' 

\n 

By dropping the terms for v/hich Vj and Vj are in the same 
dass/ this becomes a bound on the total probability of 
error . 

i 

2.3’ Sample Classification 

For the purposes of this section we can assume that the 
partition of the scene is known and we simply want to 
classify the objects. (In Chapter 3 we discuss conjunctive 
partitioning algorithms for actually estimating the 
partition.)' We shall treat each object separately, thus 
ignoring any contextual Information resulting from spatial 
relationships of objects. So v;e observe a set (sample) of 
q-dimens iona1 random variables, X = (J(]^, . . . , )ip) / from a 

common population, and our goal is to classify them. 

2,3,1 Minimum Distance (MD) Strategy 

- A Structured Approach to Classification 

A structured approach is one In which the basic, form of 

the processor is simply assumed, perhaps leaving certain 

parameters or options to the discretion of the user. A 

reasonable procedure is to choose some characteristic that 

differs from class to class, measure it for the sample-to be 

classified, and select, the class whose character istic most 

closely matches this observation. Under our assumption of 
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simple samples, each class is completely character i zed by a 
known q-d imens i onal pdf. Therefore, in MD classification, 
the n data vectors are used to estimate the pdf of the 
population, and the class is selected whose pdf is closest 
to this estimate as measured by some appropriately defined 
"distance measure" on the set of density functions. Ideally 
one would like to choose the density estimator and distance 

i 

measure in some optimum manner, but in practice the best 
guidelines are provided by experimental investigations 1141. 
Wote that a possible drawback of the MD strategy is that the 
sample size (n) must not be too small to obtain meaningful 
density estimates. 

When spatial correlation is introduced into the model 
(Appendix A), each class is only partially characterized by 
a simple q-dimens ional pdf. Although perhaps not as 
effective as a higher dimensional pdf would be, it is still 
a reasonable and valid characteristic for distinguishing 
between classes. In fact if the spatial correlation Is 
cl ass- i nva r i ant (such as that induced by the scanner), the 
q-d imens i onal pdf night be just as effective as the higher 
dimensional one. 

2.3.2 M.A.P. and M.L, Sample Classification 

In contrast to the MD strategy, the MAP strategy is a 
completely non-structured approach. The decision rule Is 
determined solely by the criterion of minimum error rate 
v/i th no a-priori restrictions. Of course a greater degree 
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2.3.2 

of statistical information is also required .'(the a-priori 
class probabilities). We can obtain the MAP decision rule by 
direct extension of 2, 2, 1,3 if we consider X ?as a 
qn-d i mens i onal random variable to be classified. Let’ x be 
the set of variates event X=x be defined 
as the joint' event i^l^.^n. Then, under the 
hypothesis XeW., the pdf of X is 


p(X=x IW. ) 


1 2 .)p(X=xlW. .) 

P(W.) j 


2 . 3 . 2.1 


n 


1 2p(W..) ii PCXm=x IW.j) 

T. T , I J .. m I 1 


P(W.) 

I 


■m I J 


J ' m=l 

The MAP decision rule can be stated as follows: 


p(X=x l'W(x) ) P(V'Kx) ) = max pCX=x 1 V/. ) PCW. ) 

i ‘ ‘ 


2 . 3 . 2. 2 


There is no minimum sample size required to implement this 
strategy. For n=l it simply reduces to MAP no-memory 
classification (2. 2. 1.3). 

f'lote that we have represented the joint pdf of a sample 
in terms of the marginal pdf of one pixel. When spatial 
correlation Is present, this is no longer a fully adequate 
representation. But as in the case of MD classification. It 
still provides a useful statistic for distinguishing classes 
v/hile avoiding the complexities of more rigorous 
re present at ions. 

So far v/e have tacitly assumed that the decision rule 
must assign the same class to all the pixels in the sample. 
With this type of strategy, either all the pixels are 
classified correctly or all are misclassif ied. Thus the MAP 



26 


2.3.2 

decision rule maximizes the average number of times that aJJ. 
the pixels in X are classified correctly. But performance 
is generally measured by just the average number of pixels 
in X that are classified correctly each time. Vie can show 
that the MAP decision rule maximizes this criterion also. 
Any decision rule that we adopt must assign a class to 1| 
for any - event X=x. Vie denote this mapping by WjCx). Let 
Z(‘) be an indicator function, i .e. a zero-one random 
variable which assumes the value 1 if and only if the event 
specified in the arguement actually occurs. The number of 
elements correctly classified in the sample is given by the 
random variable 

M = i Z(X. eW.(X)) 2.3.2. 

i=l ‘ 

E(M) = 2 E(Z(i, £Wi(X))) - 2 PCI, eW,(X)> , 

i=l ' i=l 

= X f PCX- ^ Wj (x) lX=x)pCX=x) dx 

i=lJ ‘ 

xeR^^" 

The integration implied here Is a qn-dimens ional one. Note 
that the event X- eW. (x) is equivalent to X e W- Cx), so a.ll 
terms of this summation are identical, with the possible 
exception of the decision function. Thus the -decision 
function which maximizes one term also maximizes the others. 
This confirms that the optimum decision rule assigns the 
same class to al 1 the elements of the field. Denoting this 
decision function by WCx), we have 
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E(N) =. nJPCX e W(x)lX = x)p(X=x) clx 2.3.2. 4 

ThiS/ of course/ is maximized by the MAP strategy (2, 3, 2. 2). 

The ML strategy follows directly from the MAP strategy 

by dropping the a-priori probabilities. The result is 

pCX=xlV/(x)) = max pCX=x[W ) 2. 3. 2. 5 

? ■ i 


2_.-3.3 G.H.L, Sample Classification 

We can obtain the QML decision rule by direct extension 
of 2.2.3. 1. The result is 

pCX=xlVCx)> = max max p(X=x|Wf.) = max p(X=xlVj) 

I j i 

■ 2. 3. 3.1 


VJe can also bound the probability of error for classifying 
simple samples. The analysis of Section 2.2,4 carries over, 
directly when Ji is replaced by X and the moment generating 
function is recomputed as follows: 
n 


R 


fj(x) - n 


2 . 3 . 3. 2 


V 


= S L,j<V 

m=l 


This is a sum of independent/ identical-ly distributed random 
variables. Thus 


n 


ECexp(uL|j(X))!XeVj) = H ECexpC uLj j (X^,) ) | X^^ e' Vj ) 

m=l 

= CECexpCuLjj CX)) IXsVj))*^ 

= CcbjjCu))" 

Equation 2. 2. 4. 7 becomes 

FjjCT) f P(R|jCX) > TIXeVj) s< ( 4 > j j ( .5 ) ) ‘^/TT 


2 , 3 . 3. 3 


2 . 3 . 3. 4 
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It Is a property of moment generating functions that 

^ (Jj.jCC) = 1/ so this bound is an exponentially 
decreasing function of n when cj)|j(.5) f 1, or equivalently 
when the Bhattacharyya dis^tance is non-zero. Thus the 
probability of error for the GML sample classification 
strategy is bounded by a sum of exponentially decreasing 
functions of the sample size. 

To illustrate hov7 powerful this bound can be we novj 
consider a simple example. Suppose that the ItJi and jih. 
spectral class • densities are as depicted in Fig. 2. 3. 3.1. 
The mean vectors are equal/ whi ch 'resul ts in a high degree 
of "overlap". Therefore the Bhattacharyya distance is only 
G.ll/ and 4)jj(.5) =70^ = 0.8944. The actual conditional 

error rate/ Fjj(l)/ for no-memory classification (n=l)/ Is 
50%/ which represents very poor performance. This implies 
that a "polling" classifier also has a 50% error, rate 
regardless of the sample size. But Figure 2. 3. 3. 2 shows how 
the Gf^L performance improves as the sample size increases. 
For a sample of just 40 observations the error rate is 
practically insignificant. Although we probability cannot 
expect such dramatic performance In practice . (due to the 
idealizations of our model), this still provides a strong 
motivation for our effort to apply sample classification to 
MSS data. 

2.3.4 Maximum Likelihood vs. Minimum Distance 

Let X = (Xi/.../Xp) be a simple sample from a MVN 
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population/ and defi'ne tPie statistics 


^ i=r ' 
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The maximum likelihood estimates of the mean vector and 


2 . 3 . 4. 2 

' 4 

1 


covariance matrix are: 

id = 

n 

C = 1 T (X.' -M)(X. -M)' = S<,7n - MM' ■ 
n 1=1 

The corresponding density .estimate is .given by equation 

f . ^ < 4 - . . 

2.1.4. 


Two popular distance measures are the Bhattacharyya 
distance and the divergence. If W. is a class with density 
!!(Mj ,C| ;;i) then the Bhattacharyya distance between this and 
M(M/C;><) Is given by 

B = .25 ( In I (C4CD/2 1^ * ) ' CC+C . )^ (K-H ..) ) 2.3.1..3 

ICMC, I ~ 

and the divergence can be efficiently calculated from 
D = .5 tr((C"^+C"^ )(C-+C. + (M-M.KM-M.)' )) “ 2q 2.3. 4.4 

~'^l — ~1 I 

r 

Computationally/ D is faster than B, requiring about 2q(q+2) 

i 

multiplications plus 1 matrix inversion per class for each 

(S.^/S 2 ) pa i r cl ass i f i ed. In addition to thiS/ 3 requires a 

determinant and a logarithm. (This does not Include 

quantities such as ' |C.l which can be' computed once and 

1 » 
a 

saved.) However B appears to provide an advantage in terms 
of classification accuracy/ based on -experimental evidence 
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!14|. And Its direct relationship to the' Chernoff bound 

gives B some intuitive appeal as well. 

In order for the ML strategy to be computationally 

competitive with D and B, the likelihood function must be 

expressed in terms of and $2 as follows. 

PCXIW.) = n 2. 3. 4. 5 


= ( |2nC, !*^exp( 5) (X. (Xj -M 

j = l 


-Hj ))) 


In pCXlfi:) = 


-,5Cn lnl2nC|l +Z^^j’9i^-^J Xj +M 5 ‘Cj H|)) 

j=l 


The quadratic term yields 


ix.-cTlx. = 2 - trCcT^ txjlj') 

j=l -* J=1 


In pCXlW.) T -.5trCC’^S2) - . 5 n (M| ' C .^M , +1 n | 2nC , I ) 


2 . 3 . 4. 6 


which can be computed with just ,5qCq+5) multiplications, 
once the non-data-dependent quantities have been 
initialized. Thus the ML strategy requires only 25^-50^ as 
many multiplications as D and no matrix inversions or 


determinants. 

It is interesting to express equation 2. 3. 4. 6 in terms 
of M and C. Substituting for S^^ and from equations 
2. 3. 4.1 and simplifying, provides: 
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1 In p(XIV^j) = -.5(ln|2TTC.l + tr( cT^( C+ KM-H| > ' ) )) ‘ 

n 

2. 3. 4. 7 

which we shall denote by LjCM,C). By adopting the ML 
strategy^ one is essentially using this quantity as a 
measure of the "similarity" of sample X to class V/j, just as 
B and P are used to measure their "dissimilarity". 
Therefore, -rLtCM/C) can be interpreted as a measure of 
dissimilarity between the distributions M(M,C;>() and 

I 

IKII: ,C: ; . However it is not a distance measure in the 

sense of Wacker and Landgrebe 114|, because it satisfies 
none of the three basic properties of distance measures; 

i.e. if f(i/j) Is a "distance" between distributions 
N(M:,C;;x) and NCM;/Cr;)i) then 

1. f( I,j ) >,. 0 2. 3. 4. 8 

2. f( i ,I ) = f Cj,j ) = 0 

3. f(i,J) = f(j,i) 

One can force compliance with properties 1 and 2 by adding a 
bias term as follows: 

d(i,j) = -Lj(Mj,Cj) + Lj(Hj/Cj) 2. 3. 4. 9 

- .5(1n IC< I <■ trCCf^CCj +(iij -Mj)')) -q) 

ICjl 

which can .be ’ recognized as a form of one of the 

Ku 1 1 back- Lei bl er numbers 1 19 i . Since the bias term is 
independent of i (the class number), use of this criterion 
is still equivalent to the ML strategy, as long as iC*l > 0. 

^ J 

Also, the quantity d( i , j )+d(j , i ) is equivalent to the 

divergence, v^hich satisfies all three distance measure 
propert i es . 
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The ML strategy has other compell i ng properties besides 
computational efficiency. On theoretical grounds, for the 
idealized conditions we have stated, it is the optimum 
strategy (for minimum error rate) when the a-priori class 
probabilities are equal. Also, the Chernoff bound for ML 
no-memory classification can be extended to provide an error 
bound for ML sample classification. Experimentally, under 
non- ideal ized conditions, the ML strategy does appear to be 
slightly better than MD (using 5) on the whole, although it 
Is not consistently better. The experimental results appear 
in Chapter 4, Wacker's experimental results for 
Kul 1 back- Le i bl er numbers !14[ also lend some support to this 
observat i on . 

Another important property is that for small sample 
sizes the ML strategy does not break down as do the MD 
strategies. For a sample size of 1, it merely reduces to 
no-memory classification. Finally, the summation in 
equation 2. 3. 4. 5 is distributed as chi-squared with nq 
degrees of- freedom when ’Wj is the correct hypothesis for 
sample X. Therefore it can be used to construct a 
significance test of this hypothesis. This is useful for 
detecting samples that belong to none of the specified 


cl asses . 
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CHAPTER 3 

IMAGE PARTITIONING 

Once the partition is known, powerful techniques are 
available for classifying the individual objects. Thus, 
v/hen the partition is unknown, an image partitioning 
algorithm offers an attractive alternative to no-memory 
classification.., For reasons discussed in Section 1.3, the 
algorithms considered here are the type we refer to as 
conjunctive object-seeking. We have previously described 
this approach as a progressive merging of adjacent elements 
which are found to be similar according to some statistical 
criterion. Thus an algorithm consists of statistical tests 
applied in some logical sequence. The "logical sequence" is 
the subject of Section 3.1, and the rest of the chapter 
surveys some possible test criteria. 

3.1 Partitioning Logic 

In general It is not possible to design an error-free 
partitioning algorithm. First of all, there is a certain 
amount of ambiguity in defining the "true" partition due to 
real effects such as pixels that overlap physical boundaries 
or ambigui.ty in the physical boundaries' themselves. 
Secondly, two main types of decision errors can occur. 
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leading to: (1) false boundaries/ and (2) missed boundaries. 
Also the combined effect of these two errors can produce 
^approximate" boundaries/ which is not actually a 
well-defined' category due to the ambiguity of the true 
partition. Since object size and shape are not used as 
classification features/ Type-1 errors are generally' much 
less li'kely to lead to misclassif ications than are Type-2 
errors. This philosophy accounts for certain 
simplifications in the partitioning logic. 

. The basic approach that we have adopted (due to Rodd 
till) consists of two "levels" of tests. Initially the 
pixels are divided/ by a (hypothetical) grid/ into small 
groups of four (for. example). At the first level of 
testing/ each group becomes a unit called a "cell"/ provided 
that -it satisfies a relatively mild criterion of 
homogeneity. Those groups that are rejected are assumed to 
overlap a boundary and their individual pixels are usually 
classified by the no-memory method. These groups are 
referred to as "singular" cells. At this level it Is 
usually des i rabl e to ma intain a fairly 1 ow rej action rate to 
reflect the relatively high a-pr.iori probability of a' group 
being homogeneous. The goal at this level is essentially 
the same as the goal of the boundary seeking techniques 
discussed in Section 1.3, i.e. to detect as many pixels as 

I t 

possible that lie along -boundaries v;ithout requiring that 
the ones detected form closed contours or even be connected. 
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At the second level/ an individual cell is compared to 
an adjacent "field"/ which is simply a group of one or more 
connected cells that have previously been merged. If the 
two samples appear statistically similar by some appropriate 
criterion/ then they too are merged. Otherwise the cell is 
compared to another adjacent field or becomes a new field 

• I 

Itself. By successively "annexing" adjacent cellS/ each 
field expands until it reaches its natural boundaries/ where 
the rejection rate abruptly increases/ thereby halting 
further expansion. The field is then' cl ass I f 1 ed by a sample 
classifier/ and the classification is assigned to all Its 
pixels. 

This approach has the important advantage that it can 
be implemented "sequentially"; 1 .e. raw data need be 
accessed only once and in the same order that it is stored 
on tape. This is important for practical/ rather than 
theoretical/ considerations. The flow chart in Figure 3.1.1 
indicates how it can be done. In this chart/ the top of the 
scene is referred to as north/ and the general processing 

sequence is from north to south. 

A possible drawback of the approach described above is 
that in certain hypothetical situations/ Type-1 errors are a 
certainty. For example/ a U-shaped object would develop as 
two separate' fields which expand southward and eventually 
meet at the base of the U. But since no provision is made 
for merging such fields/ a false boundary between them will 
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Figure 3 , 1«1 Basic Flow Chart for a Tsjo- L evel, Conjunctive 
Partitioning Algorithm 
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3.1 

result. Such a provision can of course be made, but only at 
the cost of additional complexity. This does not appear to 
be warranted by the relatively harmless nature of an 
isolated Type“l error. Thus the false boundary actually 
results from a design simplification rather than a true 
decision error. 

Many modifications (both large and small) to the basic 
flow chart are, ■ of course, possible. For example, the 
Level-1 test can be removed from the loop If. performed in 
advance and intermediate results saved on tape. (Tnis Is 
particularly useful in a research environment.) Another 
modification is described in Section 3.3.3. It involves 
comparing a cell to as many as three different fields at 
once, instead of one*’at“a"time» 

3.2 Unsupervised Mode 

In order to implement the sequential approach we must 
specify two test criteria corresponding to the two levels. 
In this section v;e consider ways to do this "unsupervi sed"; 
i.e. the test. cr i ter i a are independent of specific knowledge 
of the spectral class distributions. Note that our usage of 
this term Is analogous, but not identical, to the 
conventional usage. 

3.2.1 Unsupervised Annexation 

Let X = represent the pixels In a group of 
one or more cells which have been merged by successive 
annexations. Let Y = ( 1,^) represent the pixels in an 
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adjiacettt, non-singular cell. Since both X and Y v/ill have 
satisfied certain criteria of homogeneity, we assume that 
each is a sample from a MVN population. ' Let f and g 
represent the corresponding density functions. It is 
desired to test the (null) hypothesis that f = g. This is a 
composite hypothesis, since it does not specify f and g. 
The "likelihood ratio procedure" 1201 provides an effective 
statistic for testing this hypothesis. Van Trees •I21'| 
refers to it as the "generalized likelihood ratio". Let 
Hg(x,y) ** {p(x,y{f,g): g=f, feP} 

H^(x,y) = {p(x,y|f,g): fefi , ge«, gjff } 

where pCx,y!f,g) is the conditional joint density of X and Y 
evaluated at and yeR*^? and R is a set of MVN density 

functions. The assumption of cl ass-cond i t i onal . i ndependence 
enables us to express the joint density of pixels as the 
product of their marginal densities. Thus: 

n m 

p(x,ylf,g) = p(x]f ) p(ylg) = ( n ))( n sC^r )) 

i=l i=l 

The generalized likelihood ratio is defined by: 

sup HqCX,Y) max p(X|f) p(Ylf) 

A = = feR 3. 2. 1.1 

sup Hi(X,Y) max pCXj.f) p(Y|g) 

feP 

geft 

For an "unsupervised" approach to partitioning V'/e take n to 
be ^ the' following set of functions of 

MeR*^, C = symmetric, positive-def ini.tej- 
Since for any feJl there exists a geR that Is arbitrarily 
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close to f/ the condition ’ can be dropped from the 

t 

denominator of 3. 2. 1,1. Therefore: 


A = 

max N(M/C;X) l^(M^C;Y) 

1 



max 

h’(My/Cy;Y) 




max M(Ex/C;X) NCMy,C;Y) 

t 

max N(il/C^X) 

NCM/C;Y) 


max 

N(Mj^^Cx;X) N(My,Cy;Y) 

max NCiilj(,C;X) 

M(My,C;Y) 

= 

^2 

* Ai 



v;here 


n 

= HnCM/C;!,) 
i =1 



N(M,C 

;X) 

i 



iKM/C 

;Y) 

iTi 

= riNfil/C;!,) 
1=1 



and in each case the maximization 

is with respect to the 

mean 

vectors and covariance matrices. 



Anderson |221 shows that: • 



Ai = 

C 

lAI/lBI 


3. 2. 1.2 

hz ~ 

( 


,)•= 

3. 2. 1.3 

■where 





F! = 

n + 

m 




X =i;Xi/n Y =2Yj/m 

i=l i=l 

Ax = Z (li -I) Cli -I) ' Ay = I { Y j - Y) ( Y i - Y) * 

■''1=1 ~ i = l 

(In order to assure non-singular matrices with pr 
we need n > q < m.) [221 

A ™ Ax Ay 

H = (nl + ml)/N 
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n _ __ 

Bx ■= 2 = Ax + n(l-M)(X-M)' 

i =1 

ill _ _ 

By '■= 2 = Ay + 

~ i =1 

B ■= Bx + By = A + mn d-Y) (X-'Y).' 

: N 

Anderson also suggests the following modification: 

X ) ~ ^1*^2 

where Xx and ^2 are obtained from Axiand ^2 by repl aci ng* the 
number of pixels In each sample by the number of degrees of 
freedom; i.e. replace n by n-1^ m by m-1, and N by M-2 in 
formulas 3.'2,1.2 and 3. 2. 1.3. In either case; the 

f 

statistics are invariant with respect to a linear 
transformation on the data vectors.. It follows that their 
distributions’ under the null hypothesis are independent of 
the actual MVN population from which the samples are drawn. 

The test procedure is to compare X with some decision 
threshold T < ' 1 , which depends in general upon n and m. The 
hypothesis is. accepted if X ^ T and rejected if X < T. In 
the unsupervised mode/ T is determined by specifying the 
desired "size" (significance level) of the test, because the 
power of the test is indeterminate. ' in order to. do. this, 
however, the distribution of X must be tabulated, f Under 
the null hypothesis Xx and X 2 are ’ i ndependerit ly distrabuted 
1221, so we can. simplify the distribution theory and 
accomplish the same objectiv.e by the following procedure. 
Pick significance levels s-j and S 2 


such that 
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s = l'(l-si)Cl-S 2 )/ where s is the desired significance 
level. Test X| using a threshold T{ such that PCX. < TjlH) 

= sj, i=l/2, where H denotes the event that the null 
hypothesis is true. The null hypothesis Is rejected if 
either test produces a rejection. Thus the effective size 
of the test is given by 

I _ p(Xi T]^,X2 >y T 2 IH) = 1- PCXi >, TilH) P(X2 T 2 IH) 

= 1 - (l-si)(l-S2) = s 

as desired. This procedure gives us complete freedom to 

pick the ratio s^/s^ >. 0. As this ratio increases, the 

power of the test against the alternative '*Mx=My/ £x^9y" 

increases, and the power against the alternative 

C„=C '* decreases. 

A y 

We now review the distribution theory that is needed to 
implement these tests. There is a transformation of 
which, given, H, has an F-d i st r i but ion withq and (M-q-1) 
degrees of freedom 1 23 | . It is given by 

Fi = Cili ' 3. 2. 1,4 

lAl q 

Thus the test for a significant difference between the mean 

vectors can be implemented by computing Fj and comparing it 

to a threshold tj determined by the relation PCFj > t^lH) = 

5 -^. Alternatively, it can be shown 1 181 that 

C 131 -1)(N~2) = (N-2)nm(X-Y)'A"^(I-Y) = T^ 3. 2. 1.5 

{Al N ~ 

which (given H) has a T^ distribution with N-2 degrees of 
freedom 1221. T^ is Hotelling's generalization of the 
Student-t statistic, which is commonly used to test the 
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hypothesis that the means of two univariate normal 

distributions are equal given that the variances are equal. 

The following transformation of ^2 has the 
F-d i str i but ion (given H) 1 23 | : 

G = -2 In X 2 • _ 3.2.1. 6 

g = 2q^^-3q-ir 1. . + 1 " JL,.\ 

6(q + l)ln-l m-1 N-2j 

(q-l)(qf2)f 1 __ 1 ?1 2 

^ ^ l(n-l)2 (m-l)2 (N-2)2J ^ 

V = .qjj-q-tU. 

2 . • ■ 


w 1 = 


v+2 


u 


Z = 




( l-g-v/V7)G 
V 


w 


{ 


Z 

wZ 


u >0 


( I-g-^27 w>C ^ u <0 


u >0 
u <0 


v(l-Z) 

* 

F2' has an F-distribution with v and w degrees of freedom. 
Thus the test is implemented by computing Fz and comparing 
it to a threshold p 2 determined by the relation 


P(F2 > t2lK) =S2. 

Due to the complexity of the test for different 
covariance matrices, it may be desirable to rely only on the 
difference in mean vectors (assuming that a difference 
exists). A common approach is to simply assume that all 
covariance matrices are equal, thereby eliminating the need 
to test the statistic F 2 ; i.e. let S2=0. The test of Fj^ is 
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probably fairly robust under departures from this assumption 
anyway 1231. An ‘ al ternat ive approach, which does not 
require this assumption, is to use the Behrens-F i scher 


statistic defined as follows |22l: 

Z. = Xj” / 1 = 1 , 2, ...,m^ n. 


3. 2.1. 7 


Z 


Zii/m 

I =1 


D = S C2.,-Z)CZ.-l)V(m-l) 

1=1 

= mCl-X)’D'^Cl-Y) 

This has a T^-distri bution with m-1 degrees of freedom, or 

equivalently the statistic 

9, N 3. 2. 1.8 

F, = T^(m-q). 

^ (m-l)q 

has an F-dl st r i but ion with q and m-q degrees of freedom 
under the hypothes i s that r^=Ey. Thus the test is imple- 
mented by comparing Fg to a threshold tj which satisfies the 

relation P(F 3 > t3lMx“i^-y^ " 

These multivariate tests all have the same weakness as 
^';D classification, namely the problem of estimating a MVN 
density from a relatively small sample (sometimes known as 
the "dimensionality" problem). This led to the constraint 
n > q, a condition which is often not met. Even when the 
condition is net, poor estimates can result, leading to 
decislon'errors. One approach to this problem is to reduce 
q by deleting features. it is well-known, for example, that 
a subset of features used to train a classifier from small 
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training samples can sometimes produce better cl ass I f.l cat ion 
results than the- full set |24I. With this approach, 

hovvever, one is faced with the problem - of choosing the 
subset . 

Another approach is to'base the decision on the q, 
univariate, marginal distributions; i.e. simply consider the 
data in one spectral channel at a time. This has been 
termed a "multiple univariate" (MUV)' approach. In each 
channel v/e test the uni var i ate hypothes is that the means and 
variances of- the two samples are equal. Since the 

boundaries may be strong In some spectral channels and weak 

in others, we accept the null hypothesis . only if the 

univariate hypothesis is accepted in al 1 q channels. 

Besides avoiding the dimensionality problem, the MUV 

procedure requires less computation and simpler distribution 
theory. Hov/ever, it must be pointed out that in situations 
v^here class separability is primarily a multivariate effect, 
the MV procedure may be more advantageous. 

In order to obtain the univariate tests we can follow 
the same development that led to the multivariate rtests 
except • that q = l and Aj^, Ay, A, and B are just 
one-dimensional matrices (scalars). Thus equations 3. 2. 1.4 
and 3. 2. 1.5 simplify as follows: 

Fi - T^ = (N-..2.).rt.m 

' N A 3. 2. 1.9 

I ’ . ' ' 

This has an F-distribut Ion with 1 and. (N-2) degrees of 
freedom, under the null hypothesis. Equivalently we can say 
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that the statistic 


J- 


(M-2)nm (X-Y) 
MA 


3.2.1.10 


has a Student-t distribution with M-2 degrees of freedom. 

The statistic ^2 simplifies to 

“ 3.2.1.11 


X2 
where 




(l+K)^“^f(r) 


K = (m-l)/(n-l) 


r .= 
f(r) 


K Ax/ Ay 

„n-l • 


(r+K) 


M-2 


The statistic r has an F-di str i but ion with n-1 and m-1 
degrees of freedom/ and it is independent of under the 
null hypothesis |221. But since f( • ) is not monoton ic, two 
thresholds must be determined in order to implement a test 
on this statistic. For a s i gn i f i cance level the 
thresholds T' and T" must satisfy 


PCr < T’ or r > T"!H) = S 2 3.2.1.12 

f(T’) = fCT”) 

Alternatively one could resort to the transformation in 


3 . 2 . 1 .. 6 . 

3.2.2 Unsupervised Cell Selection 

"Cell selection" refers to the Level-1 test, which is 
used to detect cells that apparently overlap boundaries. 
Such cells frequently exhibit abnormally large sample 
variances. Thus a possible criterion for a cell is to 



require that the sample variance in each spectral channel 

fall below some reasonable threshold. A similar approach is 

to form the ratio of the square root of the sample variance 

to the sample mean and compare I t to a threshohd (which we 

shall call "c"). This criterion has the advantage of being 

independent of the scale of the data, 

A possible multivariate approach is to place an upper 

limit on the sarnp.le generalized variance, lAv/ml, that any 
; c 

cell (Y) can have. This is equivalent to placing a lower 

limit on the value of the statistic max M(M,C;Y). But again 
■ 11/ C ~ . 

we mention that the dimensionality problem seriously weakens 

the iW approach. .It can cause very poor estimation of the 

generalized variance and increase the chance of a decision 

error . 

3.3 Supervised Kode 

In this section we -develop a way to "supervise" the 

sequential partitioning process, using the known spectral 

class’ distributions. Our approach is based on the same 

composite hypothesis testing procedure as. the unsupervised 

approach. The effect of the spectral class distributions Is 

to greatly, simplify each hypothesis, but paradoxically the 
\ • 

resultant test criterion is much more complicated. 
Fortunately, much of - the computation can be done 

I * 

"sequentially", i.e. relying on previous saved results, 

3.3,1 Supervised Annexation 

Let X and Y be samples from a field and an addacent 
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cell as In Section 3.2,1, Vie follow the same development as 
in that section/ except that for a supervised approach to 
partitioning v/e take n as: 
n = {p(2c|Vj): 1=1/2, .../k} 

where k Is the number of spectral classes. Note that this 

Is a considerably more restrictive condition than before. 

* 

The corresponding generalized likelihood ratio statistic is: 
max (p(XlVj) pCYlV;)) 

A = _i [ ^ 3. 3. 1.1 

max (p(XlVj) p(Y|V,)) 

■ i.J • ■’ 

Note that this is a multivariate statistic without the 
constraint m > q that was necessary in the unsupervised 
mode. Ilov/ever the maxima in formula 3. 3, 1.1 cannot be 
expressed in a simple analytic form as in 3. 2. 1.1. They can 
only be obtained by exhaustive search. Furthermore, the 
distribution of 3. 3. 1.1 is unknown under either hypothesis, 
because it depends on the true classes of X and Y. But in 
return we gain a statistic which should be more "sensitive** 
to the presence or absence of a boundary. This should 
produce better performance and make the specification of a 
decision threshold less critical. In fact, the experimental 
results in Chapter 4 indicate that the threshold need not be 
a function of ,n, the current size of sample X, in order to 
obtain good results. Furthermore, the results tend to be 
fairly stable -over several orders of magnitude of threshold 
variation. Thus we will find it convenient to represent the 
decision threshold as 
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T = 10"^^ t 0 3. 3. 1.2 

Unlike the unsupervlsed approach/ a constant decision 
threshold of T=l' does not imply that the null hypothesis is 
alv/ays rejected. But it does lead to the same final result 
when the GHL strategy is used to classify the objects. This 
is because A can exceed 1 only if X and Y would be 
classified the same by the GML rule anyway. Consequently/ 
the only practical values of T are those between 0 and 1, 
Lacking any distribution theory to provide guidance in. 
choosing a suitable threshold in this range, we shall rely 
instead on an empirical approach. 

Calculation of the generalized likelihood ratio 
criterion’ can be greatly simplified by the following 
measures: 

1. Change the denominator of A to maxC p( X | V • ) p( Y I V ? ) ) . . The 

i / J 

only effect' of this change, Is to cause the value of A to 

saturate at an upper limit of 1. it doe's not affect the 

/value of A when A< 1, Since T is alv/ays less than 1, 

the change cannot affect any - decisions. The 

simplification that it affords is that A can npw be 

-written as follows: 

i ' max( p< XI V| ) p( Y1 Vj ) ) 

A = i 

(max p(Xl V; ) )Cmax p(Y|V;)) . 

I J 

I 

which is -simpler to compute. 

2, Compare In(A) to In(T) instead of A to T. 
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In A = maxCln p(X!Vj) +ln p(Y|V[)) 3, 3. 1.3 

i 

-maxCln p(X|V;)) -maxCln p(Y|V;)) 
i i 

Besides, converting mul 1 1 p.l Icat I ons Into (faster) 
additions, the quantity In'p(YlVj) can be efficiently 
computed b.y formula 2. 3. 4. 6. The quantity In pCXlV^j) can 
be obtained by other means, .as v;e shall see. 

3. Assume for the moment that the log-likelihood function of 
X Is already available In a storage array G (say); i.e. 
GCI) = In p(XiVp, I=l,2,.;.,k . 

and J I s’ an 1 nteger such that G(J) = max G(i). Then 

I 

compute the first term of ln(A) using an intermediate 

storage array ,g (say), as follows: 

g(D- = G(I)'+ 1-n p(YlV]), ' 1=1,2, .. .,k. 

and j Is an Integer such that g(j ) = max g( I ) . So; 

I 

ln(A) ■■= g(j) ’-G(J) -max In p(YlVj) 

I 

■If the Indica'ti,on is to merge Y and X, then the new 
log-likelihood function of the field is just the sum of 
the log-likelihood functions of X and Y. Therefore we 
can simply update G and J as follows: 

G( i ) = ■ g( i ) , i =1,2, . . . ,k 
J =-j ■ . 


and the preliminary assumption is justified. Thus v;e 
avoid using formula 2. 3, 4. 6 to compute In pCXl.Vj). 


When the 

po i nt 

is reached that 

the 

field 

X stops 

expanding, it 

must ■ 

be’ class i f i ed. 

This 

woul d 

no rmal 1 y 

require the 

sample 

mean vector and 

autocorrelation (or 
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covariance) inatrix/ which would have to be continually 
updated as cells are added to the field. Recall however 
that the GML strategy is 

pCXIV(X)) =. max pCXlVj) = max exp(G(i)) = exp(GCJ)) 

i • i 

Therefore V(X) = Vj , so no additional updating or 

computation is required to classify the field if the GML 

strategy is used. 

1 


3.5.2 


Qj(Y) 


Supervised Cell Selection 
A useful statistic for cell 
-.1 m 

= f(£3 - 2Mj’ 

I =1 


sel ect ion 
m 


III 


i«i 


i s 

m^j 


C 



where j is such that 

In p(YlV') = max In p(YlVt) = max - .5 (m« 1 n | 2 ttC; 1 + Q:(Y)) 

j i 1 . , 

The decision rule is to accept the hypothesis that Y is 
homogeneous if QjCY) < c, where c is a prespecifted 
threshold. Otherwise the hypothesis is rejected. This 
criterion has the particular advantage that it tends to 
reject not only inhomogeneous cells, but "unrecognizable" 
cells as v/el 1 , (Unrecognizable cells are those which 

represent spectral classes that the classifier has not been 
trained to recognize.) Another advantage of this criterion 
Is that its use of the' log- 1 ikel i hood function make-s it 
especially compatible with the supervised annexation 

criterion and the GML sample classifier. 

As a final note, the distribution function 
PCQj(Y) > ciYeVj) is chi -squared with mq degrees of freedom. 
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This can be used to provide initial guidance in choosing c. 

I 

3.3.3 A1 ternat ive ■ Part i tioni ng Logic, 

The logic of Fi’gure 3. 1.1.1 compares a cell to the 
norths west/ and east-adj acent f iel ds Clf necessary) seeking 
a "match”. If a match is found/ the merge takes place 
immediately without regard to whether it is the "best" match 
or not. Another approach that is used is to compare the 

cell to all three fields at once (if that many distinct 
adjacent fields exist) and attempt to determine the best 
match. In the supervised mode a match is determined by 
comparing the likelihood ratio to a fixed threshold/, so a 
reasonable definition for the best one is the field for 
v;hich this ratio is largest. Normally the east-adjacent 
field would' not exist at the time the other two comparisons 
are made/ so its likelihood ratio is supplied by looking 
ahead"; i.e. the east-adjacent cell is compared to its 
north-adjacent field and if they match/ the current cell is 
compared to their union to obtain the likelihood ratio. 

This approach has not been used in the unsupervised 
mode/ mainly because of the difficulty of determining the 
best match. A logical approach would be to choose the field 
■for which the null hypothesis is "least fejectable"; e.g 
choose the field which maximizes the minimum significance 
level at which the null hypothesis would not be rejected. 
•In other words/ if Xjj represents the observed value of 
for the jjLh field/ then the field is chosen for which 
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min PCX( < XjjlH) Is maximum. The difficulty in actually 
i=l,2 

doing this is that the complete distribution function of Xj 
would be required. Generally it is available for only a few 
isolated s ign i f i cance- 1 evel s . 
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CHAPTER 4 

CLASSIFICATION RESULTS 

Experimental results obtained in the investigation of 
inul ti spectral image partitioning techniques are presented In 
this chapter. Several different data sets of markedly 
different characteristics are classified by these 
techniques. In many respects they represent a cross-section 
of MSS data. Both 1 ov/ altitude aircraft data and 930 km 
high LAMDSAT-1 data are included. The ground resolution 
varies from 4.6m to 80m, and the size of phys ical . objects 
varies from just a few pixels to thousands of pixels. Data 
representation is 8 bits for aircraft data and 6 bits for 
LAfIDSAT. Spectral resolution varies from 0.02 to. 2,40 
micrometers, while the number of spectral channels 
available varies from 12 to 4. The actual number of 
channels used for analysis varies betv/een 3 and 6. The 
number of spectral classes represent.! ng ground cover types 
varies from 5 • to 17, and the number of informational classes 
varies from 5 to 11. 

The results are grouped by data set rather than by 
analysis technique to facilitate the comparison of different 
analyses of a given data set. In order to provide a quanti- 
tative measure of comparison, only data sets are used for 
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which a substantial number of "test areas" are available, 
by comparing the results of a given analysis on a point~by“ 
poi'nt basis with the desired result in each test area/ one 
can obtain an estimate of the accuracy Cor inaccuracy) of 
the analysis. The larger and more numerous these; test 
areaS/ the better this estimate will be. Thus one analysis 
technique is regarded as being better than another ! if it 
tends to achieve fewer mi scl ass i f i cat ions in the test areas. 

Relative error rate is an important measure of a 
classification scheme, but It Is not the only consideration. 
Obviously speed is a desirable attribute. Although CPU 
times are compared in this chapter, it is important to 
remember that efficient coding has a lot to do with speed, 
and no claim is made that the research programs used here 
are optimized. A less tangible consideration is the amount 

of effort and experience required to use a particular 

analysis scheme. The schemes considered in- this 

invest! gat ion were designed with si mp licity In* mind, 

requiring a minimum of user input. The results in this 
chapter will help to assess the degree of exper t ence cneeded 
to, provide this input, and they provide a data base of 
experience from which to draw. 

4.1 Analysis Schemes 

Within the framework of Chapter 3, an analysis scheme 
is specified by choosing: 

1. A Level-1 option and associated parameters 
(threshold and cell size) 
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2. A Level-2 option and as.sociated parameterCs) 
(significance levels or threshold) 

3. A sample classifier option. 

It would be a hopeless, and probably pointless, task to try 
to investigate all the possible comb Lnat i ons of these three 
options. Instead the less logical combinations were 
arbitrarily eliminated' in order to concentrate more effort 
on evaluation of the remaining ones. Consequently, only 
wholly unsupervised and v^holly supervised methods are used 
For the partitioning phase of processing. (No "iiybrid" 
combinations of Level-1 and Level-2 options are considered.) 
Thus' partitioning is done in either the unsupervised or 
supervised ’’mode". Furthermore, in the unsupervised mode, 
no hybrids of MV and MUV tests are used to test first and 
second order statistics at Level-2. At Level-1, only the 
MUV ratio test (described in Section 3.2.2) is used in the 
unsupervised mode. Although the Behrens-Fi scher test 
requires only one significance level to be soecified by the 
user, it tests only First order statistics, provides 
unattractively few degrees of freedom, and requires a 
substantial amount of computation. Consequently It was 
eliminated as an option. Due to the advantages (enumerated 
in- Chapter 2) of maximum likelihood sample classification 
over minimum distance classifiers, it was the logical choice 
for the classifier option. The results of this chapter are 
based on the ML strategy. When subclasses are necessary the 
generalized ML strategy is actually used, although it is 
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refered to simply as a ML classifier. 

’ After all the above simplifications/ we basically are 
left with four schemes to evaluate and compare; 

-1. unsupervised MUV partitioning and ML sample 
cl ass i f i cat i on 

2, unsupervised MV partitioning and ML sample 
cl ass i f icat ion 

i 

'3. supervised partitioning and ML sample classification 

:4 , conventional ML no-memory classification. 

Furthermore the cell size for the first three schemes, was 
elLminated as a variable by fixing it at a constant 2x2 

pixels/ which Is the minimum size that can be used in the 

unsupervised mode of partitioning. This choice appears to 
provide a reasonable compromise between speed and resolution 
for MSS data. 

A common element of all four schemes is the process of 
"training". Jhis is the process by which each main class is 
modeled statistically with the aid of data vectors 

(patterns) known to belong to that class. If the training 
data for a class exhibits a multimodal structure/ then It is 
usually divided into two or more subclasses/ each 

corresponding to a mode. This serves two purposes. l(l) It 
enables each subclass to be modeled approximately by a MVN 
distribution which is completely characterized by a mean 
vector and covariance matrix. These can be estimated easily 
from the data vectors assigned to that subclass. (2) Data 
from a single physical object is usually reasonably un^i-modal 
and symmetrical In distribution. Often/ those obj ects- v/h ich 
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are multimodal can be divided into a few smaller objects 
which are unimodal , Although multimodal training data may 
be representative of a particular main class as a whole# it 
is. not representat i ve of the individual objects v;hich 
compose that class. Since it is the individual objects that 
must be dealt with# the definition of unimodal- subcl asses is 
a logical step to take. in other v-^ords# pCxjWj ) (eqn. 
2, 3, 2,1) cannot be expressed in terms of p(x,|V/| ) (eqn. 
2.1.2), Each component (mode) of p(i<iV/| ) must be knov/n. 

The training and test data for a given scene compose a 
set of labeled observations which we shall refer to as 
"reference data". There are many possible methods of using 
a finite amount of reference data to train a classifier and 
estimate its error rate. Theoretically the best training 
(I.e. the lov/est error rate) is obtained by using al 1 the 
available reference data for training. If the same data is 
used for testing# this is called the "C-method". 
Theoretical and experimental results indicate that# for the 
Sayes classifier at least# the C-method produces an 
optimistic (negatively biased) estimate of the error rate; 
but the bias and variance of the estimate decrease roughly 
as the reciprocal of the number of observations used 1181. 
In contrast# the "U-method" requires test data to be 
independent of training data. The most common procedure 
(called "sample partitioning") is to use a relatively small 
proportion (p) of the reference data for training and the 
remainder for testing. In this case the error estimate is 
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unbiased and its variance decreases as the reciprocal of the 
number of test data; but the actual error rate tends to be 
larger than vvith the C-method, and its variance is times 
larger than the error variance by the C-method. 

The interpretation of results is usually somewhat 
easier for the C-method, because the question of whether or 
not the training is "representative" of the test datb does 
not arise. For comparative purposes our interest 'is in 
relative (rather than absolute) performance, so the bias 
Induced by the C-method tends to cancel out. There is no 
reason to believe that the bias v/ould be significantly 
greater for one scheme than for the others. ? 

On the other hand the U-method is routinely used in 
conventional analysis work w'here absolute performance is 
emphasized. Effective representation is- obtained in most 
cases by — using a fairly large training data set that 
consists of observations drav/n from the same general regions 
as 'the test data. For some of the data sets used In our 
Investigation, reasonably good training statistics are 
available from previous conventional analyses. By i using 
available training we obtain results with minimum effort, 
and the results relate directly to those obtained by 
conventional methods. 

For those data sets where previous training is 
unavailable or inadequate, neither the C-method nor the 
U-method Is used. Instead, the available test areas are 
sampled at an interval sufficient to provide a reasonably 
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large training set. Like the C-method this method is 
simple, it produces representative training, and it 

eliminates human bias in selecting the training set. It 

also Induces considerably less bias Into the error estimate 

i 

than the C-method does. Of course, once the training set is 
obtained, feature selection and subclass definition may have 
to be done before training Is complete. An example of this 
process is described in the next section. 

4.2 Run 71052800 - Crop Identification 

? 

MSS data collected over a particular region at a 
particular time and stored on digital magnetic tape is 
catalogued by "run” number. Run 71052800 Is a set of 12 
channel data collected over flightline 221 in Indiana on 
August 12, 1971 during the 1971 Corn Blight Watch Experiment 
[251, The correspondence between channel numbers and 

spectral bands Is indicated in Figure 4.2.1. Channels 1-7 
cover the visible portion of the speptrum, 8-11 lie in the 
reflective IR portion, and 12 is a thermal IR sensor. It is 
evident from the figure that there is a considerable amount 
of redundancy in the coverage of the visible spectrum. In 
other words, the data in channels 1-7 will tend to be 
strongly correlated, causing the information in these 
channels to be rather redundant, at least more so than in 
the IR channels. 

The area covered by this run is a rectangular strip of 
agricultural land about 1.6 km wide and 13.8 km long. It is 
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sampled 222 times along its width and 1374 times along its 
length. The scanner was carried by aircraft at an altitude 
of > 1524m with an instantaneous f i el d-of-vi ew of three 
mi n i radians . 

This data set vsas chosen for analysis' for several 
reasons, .Cl) A large number of test areas, containing 
84,855 pixels, were available from a previous crop 
identification study |25[, (2) The complexity of the 
classification is high , providing an opportunity to see how 
well the new techniques perform in such a situation. (3) 
The data set contains a combination of some very challenging 
classes and intermediate classes, as well as some easy 
classes to Identify. The 11 main classes arec corn, soy, 
v/heat (mostly harvested), rye, hay, . lespedeza (a grass), 
pasture, wooded pasture, forest, idle fields, and non-farm. 
The latter two categories tend- to be "catch-^aH" types which 
are characteristically difficult to identify by conventional 
methods. The reason for this will be discussed later. 

No previous training statistics were available for this 
data set, so training data was obtained by sampling the test 
areas as per Section 4.1. The resultant ratio of the amount 
of test data to the amount of training data is approximately 
5:1. The LA.RS system (LARSYS) STATISTICS processor 126! was 
used to compute the statistics of each main class, and 
SEPARABILITY 1261 was used to compute transformed divergence 
values for every pairwise combination of main classes and 
every combination of 6 channels. The purpose of this is 
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“feature selection";' i.e. due to high correlation between 
channels^ it is usually possible to find a subset of the 
available channels that discriminates between the’ main 
classes almost as well as the- complete set. Typically only 
3 or 4 of the 12 available channels are used to analyze 
aircraft scanner data, resulting in a large time savings.- 
Fon the present study, on the basis of the transformed 
divergence results, the best set of 6 channels is 

(2,.4,10,11,12) plus either channel 8 or channel 9. 

(Recalling the discussion of Figure 4.2.1, this result is 
not surprising.) Channel 8 maximizes the r^ge. 

transformed divergence (averaged over all class pairs,), and 
channel 9 maximizes the mini mum transformed divergence for 
any pair of classes. The difference between the two is 
slight, so channel 8 was selected arbitrarily. Based on 
histograms of the training data it Vifas decided to sub,divide 
some of the training classes. The LARSYS CLUSTER processor 
1 26 1 was used to cluster these classes into 2 or 3 modes and 
SEPARABILITY was used to determine the divergence between 
modes. Histograms, cluster quotients, and d i vergence wal ues 
were examined to determine if the modes of each clasis were 
distinct, and if two modes were not distinct, they were 
recombined' (pooled) into a single mode. The final result is 
2 subclasses each for corn, soybeans, lespedeza, and Ldle, 3 
subclasses of pasture, and only 1 "subclass" for each ,'of the 
remaining classes, a total of 17 spectral classes. The 

LARSYS MERGESTATIST! CS processor was used to merge them into 
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a single LARSYS statistics deck. Expecting to further 
reduce the number of channels needed/ SEPARASiLITY v;as again 
applied. it was found that although most class pairs can be 
distinguished on the basis of some set of 4 or fewer 
channels/ there is no one set of 4 channels which can 
adequately do this for all class pairs. Thus It was decided 
to do 6-*channel classification. This ended the training 
phase of the analysis. 

Next the data set was classified by a number of differ- 
ent schemes. The LARSYS CLASS I FYPO I NTS processor I26-I was 
used to perform ML no-memory classification. The LARSYS 
SAMPLECLASSI FY processor -1261 was used to perform minimum 
distance ( Bhattacharyya ) sample classification of the test 
areas. In the latter case the processor is essentially 
given a-priori knowledge of the boundaries. Maximum 
likelihood sample classification of the test a-reas was 
accomplished by modifying the SAMPLECLASS I FY software. To 
avoid confusion this processor will be referred to as 
'‘SAMPLECLASS I FY (ML)"/ and the minimum distance version will 
be. referred to as "SAMPLECLASS I FY (MD)". 

Unsupervised/ MUV partitioning and supervised 
partitioning schemes were implemented using 
LARSYS-compat ible processors that • were developed 
specifically for this investigation. The unsupervi sed, MV 
version cannot be used for this data set/ because it 
requires that the number of channels be less than the cell 
size/ which is not the case. The results of all these 
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analyses are presented in the remainder of this section. 

Figure 4.2.2 shows ci ass if i'cat ion performances achieved 
by ! the various processors indicated above. Note that 
processor #2 Is equivalent to cell selection without 
annexation. Thus comparing the results of #2 with the 
results of #1 (CLASS I FYPO I NTS ) gives a good indication of 
the effectiveness of Level -1 alone. And comparing the 
results of #2 with the results of #3 and #4 indicates the 
effectiveness of just the annexation CLevel~2) phase of 
processing. \ 

Also note that processor #5 should give about the same 
results as if . the entire partitioning phase were done 
flawlessly. Thus one can think of the results of #5 as a 
performance "goal”. This goal, however, is not a strict 
bound (more on this later). >’ 

Both "average" and "overall" error rates are shown in 
Figure 4.2.2. The former Is just a straight average (over 
all classes) of the observed class-conditional error rates. 
The latter Is a weighted average, where the error rote of 
each class is weighted by the proportion of test data in 
that class. These proportions are given in Table A. 2.1. 
Assuming that these proportions coincide roughly with th.e 
actual porportions of the classes in the data set, then the 
overall error rate can be taken as an estimate of the 
unconditional probability of error. 
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PROCESSOR TYPE 


Processor Key: 

#1 CLASSIFYPOINTS 
#2 Supervised Partitioning, t“0 
#3 Optimum Unsupervised Partitioning 
#4 Optimum Supervised Partitioning 
#5 SAMPLECLASSIFY (ML) 

#6 SAMPLECLASSIFY (MD) 


Figure 4.2,2 Classification Performance vs. Processing Scheme 
- Kun 71052800 
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Table 4.2.1 Relative .Influence of Each Class on Overall 
Performance - Run 71052800 


Percentage of Total 



Corn 

43.1 

Soy 

22.4 

Wheat 

17.6 

Idle 

4.9 

Non-Farm 

4,3 

Lespedeza 

3.0 

Pasture 

1.6 

Hay 

1.5 

VJooded Pasture 

0.8 

Rye 

0.5 

Forest 

0.4 


J 
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The class-conditional error rates are given' in Figure 
4.2.3. The -.results are grouped by class for easy 
comparison. Several observations are worthy of mention at 
this point. 

Observation 1 

Both the unsupervised and supervised modes are 
effective at reducing the error rate. As expected, the 
supervised mode has a fairly consistent advantage. It 
performed bet.ter' for 8 of the 11 classes, and its ‘average 
and overall error rates are lower. The actual reduction in 
error rate due to the supervised mode is 10.9^ (average) and 
8.1% (overal 1 ) ., 

Observation 2 

As one would expect, the ’ rel at i ve effectiveness of the 
ECHO approach is highly class dependent. The effect varies 
from slight degradation for some classes to vast improvement 
for others. 

. Observation 3 

The classes where the greatest gains are made are 
wheat, wooded pasture, idle, and non-farm. It has already 
been obse.rved that the latter two are "catch-a.ll” categories 
which are, typically difficult to identify using 
CLASS 1 FYPOI NTS . The reason for this is that such classes 
tend to have relatively broad probability density functions 
V(/hich overlap with, those of other classes but’ at a lower 
likelihodd level. Recalling the case of Figure 2. 3. 3.1, the 
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in CLASSIFYPpiNTS 
in Supervised Partitioning, 
in Optimum Unsup eryised 
#4 Optimum Supervised 
in SAMPLECLASSIFY (ML) . 
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Figure 4.2.3, continued 
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PROCESSOR TYPE 


processor type 
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conditional error rate for the "broad" distribution Is 50 %, 
v/hil>e the conditional error rate for the "narrow" 
distribution is just 17. 4?^. If other classes overlap the 
"tails" of the broad distribution/ then this discrepancy 
becomes even greater. But when a sample of data from the 
broad distribution is made available for classification/ it 
usually consists of a mixture of .values both near and far 
from the mean. This makes it possible for the classifier to 
determine the correct classification of the sample. 

; The wooded pasture class also has a relativel y:broad 
d is-tr i but ion due to its composition and spatial texture. 
Note that this does not necessarily imply that wooded 
pasture is statistically i nhomogeneous .or bimodal. We refer 
to it as a "compound" class. (See Appendix C.) In this 
panticular case it might at first appear that the method of 
using test~area's to evaluate performance is biased in= favor 
of:sample classifiers. In other words/ the large error rate 
observed for CLASS I FYPO ! NTS may be due to the assignment of 
many test points to the classes forest and pas tu re/ " wh I ch 
may actually be accurate .labels for those particular pixels. 
This arguement is invalid on several counts. First of all/ 
the number of test points classified as forest or pasture 
accounts for only 12 of the 44 percent error rate. The 
reduction in error rate brought about by supervised 
partitioning is 40 percent/ or 5.3 times as great as the 
maximum possible error attributable to this cause. 
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Secondly^ whenever CLASS I FYPO I NTS classifies an area in a 
"salt and pepper" manner^ the Information is highly 
unreliable. If the area actually were that way# Premise A 
(Section 1.2) v;ould be violated. Thirdly, even if valid 
po i nt-by- po i nt classification were possible, most analysts 
are not interested in the actual classification of each 
individual pixel. Instead their goal is to produce a "type 
map" which consists of a partition of the region with a 
general label assigned to each element of the partition. An 
element containing a mixture of trees and pasture for 
example would be labeled "v^ooded pasture". 

These points are illustrated in Figures 4.2.4, 4.2.5, 
and 4.2.6. Figure 4.2.4 shows a section of Run 71052800 
(lines 101-300 ) that has been classified by CLASS I FYPO I NTS . 
Each class has been assigned a gray level and displayed 
electronically to form the image. The "classification 
noise" is readily apparent. In contrast to this. Figure 
4.2.5 shows the same section as classified by ECHO 
(supervised). The random errors have, for the most part, 
been eliminated. This map is much closer to the desired 
result than is the CLASS I FYPO I NTS output. Figure 4.2.6 
shows the centers of these two maps in greater detail. Each 
class is represented by an assigned symbol (or blank), and 
each symbol represents one pixel. The four rectangular 
areas are test areas designated as wooded pasture. This 
class is displayed as a blank space to emphasize the 



Figure 4.2.4 Gray-Scale-Coded Classification Map 
- Produced by CLASSIFYPOINTS 



Figure 4.2.5 Gray- Scale- Coded Classification Map 
- Produced by ECHO 
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Figure 4.2.6 Logogrammatic Classification Maps 
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contrast between it and the others. The diversity of 
symbols in the test areas testifies to the inadequacy of 
CLASS I FYPO 1 NTS for classifying such textured regions. Most 
ofithis confusion is avoided by the ECHO technique. 

The wheat class too haS' a broad distribution/ probably 
due to the fact that the wheat is mostly iiarvested., 
I'Jhatever the cause, it adds further support to the arguement 
tha't classes with broad distributions tend to benefit the 
iiios't by sample classification. To clarify this 'point 
further, the classification improvement is plotted in Figure 
4.2,7 vs. the common logarithm of the generalized variance. 
In the case of a class with subclasses the average 
generalized variance is used. For this data, the 
correlation between these quantities is 0.81. 

Of course a broad distribution does not necessarily 
Imply that partitioning and- sample, classification: will 
produce dramatic improvements over CLASS I FYPO I'NTS . For 
example, another class may have about the same distribution. 
In which case no classification scheme can reliably 
distinguish between them. Or the class may be so unlil<;e any 
other class that CLASS I FYPO I NTS leaves no room for 
Improvement. Also the broad distribution may be caused by 
inadequate training (i.e. not representative). In which case 
accurate classification may be impossible until the training 
Is corrected. Obviously the mechanisms which , affect 
classification performance in a multidimensional. 
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Figure 4.2.7 Class Improvement vs. Log- Generalized Variance 
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mu'li't i cl ass/ mul t i -subcl ass situation- are very complex. 
Obsservation 3 provides only limited insight Into the overall 
process . 

’ Observation 4 

The supervised ECHO results for class idle actually 

surpass the performance "goal" set by processor #5. A 
! ■ / 
conceivable, explanation for this Is that idle test areas may 

4 

^ _ ' 
actually consist of several physical objects containing 

* _ c 

different subclasses of Idle. Since the ECHO processor can 

I 1 

classify such objects separately, it can actually provide an 

f 

advantage over SAMPLECLASS IFY (ML) which must classify each 
test area as a whole. 

Observation 5 - 

As expected, processor #2 CLevel-1 partitioning) can 
provide a fairly significant degree of improvement' on Its 
own. Again the effect* Is strongly class dependent./ The 
effect would probably be much greater if not for the 
correlation that exists between adjacent pixels. 

to 

The main parameter that is required for the supervised 
mode is the annexation threshold, t. Figure 4,2.8 shows how 

t . : 

the average error varies for seven values of- t.- Of these, 
the optimum value Is t=5, although all values tried gave 
significantly better performance than CLASS I FYPO I NTS ,, The 

Level-1 threshold, being of much -lesser importance, was not 

* * » 
varied in this study. It v/as previously established at c= 

90 by processing a small subregion of the data set several 
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Figure 4.2.8 Effect of Annexation Threshold (t) on Average Performance 
- Run 71052800 
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times while varying c. This value provides a sufficiently 
low rejection rate that the occurrence of singular cells is 


limited mainly to patterns 

resembl i ng 

boundary 

1 i nes 

(as 

des i red) . 

Thus the 

class i f i cat i on 

■ results 

are 

not 

necessar i 1 y 

optimized 

wi th 

respect to 

C/ but 

they 

are 


believed to be near that optimum, - 

Figure 4.2,9 shows the behavior of the overall error 
vs, t. It Is very similar to the average error'except that 
the results are shifted downward due to the heavy influence 
of the corn. class. Its minimum also occurs at t-5 . : 

The analogous results for the unsupervised mode cannot 
be presented as easily because the performance is a function 
of two variables/ the significance levels. For the same 
reason/ the optimum performance cannot be determined as 
easily. Tables 4.2.2 and 4,2.3 give average and overall 
error rates for eight different combinations of significance 
levels. The Level-1 threshold was .maintained at a constant 
c*=;0.25. Cells found, to be singular at this level were 

classified as small samples rather than as individual 
pixels; i.e, "cel 1 -spl i tt i ng” was not in effect. The best 
results occur at about S 2 = .005 or ,001 and si= .001. 
Possibly a lower value of si would produce better results/ 
but this is beyond the capability of the current processor. 
Figure 4.2.10 shows how the processors compare with re- 

i 

gard to both error rate and CPU time. In terms of time, the 
unsupervised mode is the fastest by far because it performs 
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Figure 4.2.9 Effect of Annexation Threshold (t).on Overall Performance 
- Run 71052800 




Table 4.2.2 Matrix of Average’ Error Rates iZ) for Eight 
Combinations of Significance Levels 
Run 71052800 
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4.2.3 Matrix of Overall Error Rates C%) for Eight 
Combinations of Significance Level.s 
- Run 71052800 
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Figure 4.2.10 Error Rate and CPU Time for Four Classification Schemes ‘ 
- Run 71052800 
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the fewest classifications. CLASS !'FYP0 1 NTS is the slowest 
because it performs the most. The supervised mode is in 

j ' 

between and provides the lowest error rate. 

One of- the factors Influencing CPU time is the size (in 
pixels) of an average object, since the larger- the objects 
the greater the number of pixels that can be classified at 
one time. A rough indication of this factor is obtained by 
dividing the number of test areas into the total number of 
test pixels. As indicated. on Figure 4.2.10, for this data 
set an average, test area (object) Is equivalent to a square, 
17 pi xel s wi de . 

( 

4.3 Run 72064412 - Classification Of Satellite Data 

Three LANDS AT passes over a region In Indiana on 
different dates were combined to produce this data set. 
Only data from the first date, August 25,1972, is used for 
analysis in this study. Four spectral channels are 

available on LANDSAT-1. The spectral bands are indicated in 
Figure 4.3.1. The instantaneous field of view for the three 

j 

visible band channels is 86 microradians. The region 
covered by the data set is a rectangle 45,1 km wide and 53.1 
km long. it iS' sampled 804 times along its width and 675 
times lengthwise, for a total of 541,092 pixels. A 21.4 km 
by 43.5 km subregion. (containing 2-10,100 pixels) was 
analyzed. This region was previously the subject of a study 
of strip mine activity (unpublished). The analyst provided 
both training statistics and test areas.* Briefly/ the 


* Courtesy of John Berkebile, LARS. 



Figure 4.3.x Correspondence Between Channel Number and Spectral Band for LANDSAT-1 Scanner 
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metihod of training was to perform both a manual analysis of 
designated training areas, using maps and aerial 
photography, and an unsupervised clustering analysis of the 
MSS data corresponding to each- training area. The manual 
analysis was used to associate each cluster class with an 
informational class, and the statistics of the cluster 
classes became the training statistics. Training and test 
areas have no pixels in common. The number of test pixels 
is i9512, a fairly large number for a LANDSAT analysis. The 
main classes are agriculture, forest, recent-mine; pit 
(containing -various' grades of water), revegetated ,mi.ne,i 
residential, and clouds. The numbers of subclasses are 3,2, 
4,2,1, 1, and 1 respectively, a total of 14. Test areas were 
supplied only for the first 5 of 7 main classes. All 4 
channels are used in the analysis. 

Figuf^4.3.2 shows the average and overall error:rates 
ofithe various processors. The weights for overal 1 : error 
are listed in Table 4.3.1. The class-conditional error 
rates are given in Figure 4.3.3. As for Run 71052800, the 
number of channels used is too large to permit unsuperVi sed, 
!^:V processing. On the whole, the results appear : quite 
similar to those of Run 71052800, in spite of the 
considerable differences between the two runs. Both the 
average and overall error rates are significantly reduced by 
the ECHO techniques, with the supervised mode providing a 
cons i stent -advantage over the unsupervised mode. 
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PROCESSOR TYPE 


Processor Key: 

#1 CLASSIFYPOINTS 
#2 Supervised Partitioning, t=0 
#3 Optimum Unsupervised Partitioning 
#4 Optimum Supervised Partitioning 
#5 SAMPLECLASSIFY (ML) 

#6 SAMPLECLASSIFY (MD) 


Figure 4.3.2 Classification Performance vs. Processing Scheme 
- Run 72064412 




90 


Table 4.3.1 

Cl ass 

Forest 

Agr icul ture 

Recent Mine 

Revegetated 

Pit 


Relative Influence of Each Class on Overall 
Performance - Run 72064412 


Mi ne 


Percentage of Total 
Test Pixels 

42.1 
24.6 

19.1 

11.0 

3.2 



ERROR RATE (%) 
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Figure 4.3.3, continued 
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The behavior of the pit class- is misleading and 
requires further discussion. Normally water is one of the 
easiest classes to identify using CLASS I FYPO I NTS, yet the 
17.6?^ error , rate is. by far the highest of any class. 
Therefore it is apparent that something is wrong with either 
the original training statistics or test areas. As it turns 
out. Premise "A" (Section 1.2) has been violated, and this 
has caused the pit test areas to contain many pixels that 
overlap bo^h pit and mine or revegetated mine classes. 
Consequently the pit class results are not truely indicative 
of performance, but they are included here for the sake of 
completeness. This accounts for the relatively high average 
error rate of the unsupervised mode, which is still lower 
than the average CLASS I FYPO I NTS error rate. With the 
exception .of this class, it can again be said that no class 
is significantly degraded, while some are greatly improved. 

The recent-mine class is another that bears comment. 
Notice that processor #5 performs only slightly better than 
CLASS! FYPO I NTS. Consequently one cannot expect partitioning 
to improve the accuracy, since the classifier seems unable 
to use the information effectively. The root of the problem 
IS that about half of the test areas in this class are 
actually labeled "partially revegetated mine", and some of 
these appear as revegetated mine to the sample classifier. 
Apparently these categories are spectrally too simil.ar to 
distinguish reliably even on a sample basis. Confusion of 
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thiis sort is also the source of, almost all the recent-mine 
miscl ass i f I cat ions that CLASS 1 FYPOI NTS makes, so there is 
very little other type of error for the sample classifier to 
correct. 

The performance of the supervised mode is again plotted 
for seven values of the annexation, threshold, t (Figures 
4. 3,. 4 and 4.3.5). Again the optimum value Is t=5. Notice 
that the overall error approaches both SAMPLECLASS I FY 
results quite closely. 'For this study the Level— 1 threshold 
v^as held at a constant c*= 55. This choice reflects the 
lower number of channels and higher incidence of singular 
cells as compared to the aircraft data. 

The effect of c was briefly investigated when for t= 4, 
the values c= 55 and c= 80 were compared. The higher value 
causes slight improvements In the agriculture and forest 
classes and slight declines in the others. The overall 
error rate is unchanged however. The effects are noti large 
enough to be of any serious concern here. 

The performance of the unsupervised mode is given in 
Tables 4.3.2 and 4.3,3 for seven different combinations of 
significance levels. Of these, the best combination is 52=^ 
.001 and Sj.= ,005 . The Level-1 threshold was , maintained 

. ‘ i 

between 0.20 and 0.25, and cel 1 -s pli tt i ng v/as no.t In effect. 
For. comparison, the data was also processed v/ith cellrsplit- 
ting in’ effect using the "optimum” sj and S 2 above. The pit 
class performance improved to about the same' level as the 




Figure 4.3.4 Effect of Annexation Threshold (t) on Average Performance 
- Run 72064412 
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Figure 4.3.5 Effect of Annexation Ihreshold (t) on Overall Performance 
- Run 72064412 


97 


Table 4,3.2 Matrix of Average Error Rates C%) for Seven 
Combinations of Significance Levels 
- Run 72064412 



Table 4.3.3 Matrix of Overall Error Rates C%) for Seven 
Combinations of Significance Levels 
- Run 72064412 
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oth'er processors/ but this is overshadowed by degraded per- 
formance in 'the other classes. In contrast to this,ia non 
celil -spl i tt I ng version of the superv i sed mode was tried with 
t= 4 and c= 55, and the class-by-class performance was un i - 
f orml V worse than with the cel 1 -spl i tt i ng version. . This 
might imply that the supervised Level-1 test is more 

' i 

effective than the unsupervised Level-l test or simply that 
a better Level-1 threshold exists (for the unsupervised 
test) than the one- used here. The evidence is inconclusive 

. t 

on this minor point. . 

Figure 4.3.6 shows how the processors compare with 
regard to both overall error rate and CPU time. As ‘before/ 
the unsupervised mode is fastest/ and the supervised mode is 
most accurate. Due to the greatly- reduced number of pixels 
per physical object compared to aircraft data> the 
difference in speed for the three ^processors is much less 
significant . 

4.4 Run 71052501 - Forest Cover Mapping 

Corn Blight Watch Experiment flightline 218 is a 1,6 x 

16.1 km strip of land in southwestern Indiana. in contrast 

> 

to the relatively flat flightline 221/ it Is on a "maturely 

r 

dissected/ westv^ard sloping plateau characterized by 

T 

abundant stream valleys and a wel 1 - i ntegrated drainage 
system." "Most of the land area Is in: slope/ with flat/ 
narrow ridge tops and steep valley walls." 1271 

Consequently row crops (corn and soybeans) are in the 
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CPU Time (minutes) 

Figure 4.3.6 Error Rate and CPU Time for Four Classification Schemes 
- Run 72064412 
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4.4 

fninority compared to forage crops. About 60% of the land Is 
covered by hardwood forest with a feW/ small stands of white 
pine. MSS data was collected over this region on the same 
date and by the same means as Run 71052800 and labeled 
71D52501. The region was sampled 1605 times lengthwise and 
222 times across its width. In contrast to 71052800/ a sun 
angle correction transformation v/as applied to the data, 
Thjs data set v/as previously the subject of a forest* cover 
mapping study 1271/ and the analyst’s training statistics 
and test areas are used- in the present investigation. Six 
main classes are considered: deciduous forest/ coniferous 

forest, water, forage, corn, and poy, A composite class 
Cforest and forage) was deleted ■ by the analyst, who 
recognized the inability of CLASS 1 FYPO I NTS to handle such 
data adequately. Our previous results on wooded pasture 
Indicate “that this v/ould have been unnecessary if ECHO 
techniques had been available to him. In contrast to the 
Previously described analyses, the training is as simple as 
possible. Involving no clustering or subclasses,- The 
available reference data (53,516 pixels) was simply divided 
into non-overlapping test and training areas at a ratio of 
about 13:1. Based on the transformed divergences, the best 
set of 3 channels is (6,10,12), which is used in the present 
investigation. 

Figure 4,4,1 shows the average and overall performances 
of the various processors. The v/eights for overall 
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#2 Supervised Partitioning, t=0 #6 SAMPLECLASSIFY (ML) 

#3 Optimum MV Unsupervised #7 SAMPLECLASSIFY (MD) 
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penformance are given in Table 4.4,1. The original analyst 
has’ indicated that these are in roughly the same proportion 
as :the actual occurrence of the classes in the’ data set. 
Thus the overall error rates are estimates of the actual 
probabilities of error. The class-conditional error:rates 
are' shown in .Figure 4.4.2, Several observations can be 
made: 

Observation 1 

Again the average and overall error rates are 
s ignif icantly reduced by the ECHO techniques, with the 
supervised mode providing the greatest advantage. 

if 

Observation 2 

This is the only data set analyzed to which the MV 
version of the unsuperv'i sed mode (processor #3) Is 

applicable for the number of channels deemed necessary-. Its 
performance can be described as erratic. it performed much 
better than CLASS I FYPO t NTS for the deciduous class and much 
better than all processors (except SAMPLECLASS I FY) for the 
forage class. However it had an especially difficult time 
distinguishing between soybeans and corn. • 

Observation 3 

With the exception , of the coniferous class, the super- 
vised mode (processor #5) again equaled or greatly improved 
upon the. performance of CLASS I FYPO I NTS , And' the 

t 

unsupervised, MUV mode did almost as well. In regard to the 
"apparently" poor coniferous performance, there are 
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Table 4.4.1 Relative Influence of Each Class on Overall 
Performance - Run 71052501 

Percentage of Total 


Class 

Test Po 

Deciduous Forest 

64.8 

Forage 

23.6 

Corn 

5.4 

Soy 

5.4 

Water 

0.7 

Coniferous Forest 

0.2 


Table 4.4.2 Classes Ordered By Generalized Variance 

Common Logarithm 

of General ized Improvement Over 
Class Variance CLASS I FYPO 1 NTS 


Forage 

8.02 

+9.7 % 

Deciduous Forest 

6.06 

+9.9 % 

Water 

6.01 

0 % 

Soy 

5.78 

+ 0.5 % 

Corn 

5.01 

+ 

o 

• 

Coniferous Forest 

4.33 

-4.6 % 



ERROR RATE (%) 
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Figure 4.4,2 Performance By Class (Run 71052501) 
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extenuating circumstances to consider. It has already been 
mentioned that only a sma 1 1 stands of v,/hlte pine are 

knov/n of in the entire data set. The original analyst used 

♦ 

most of these pixels to ensure good training statistics, 
leaving only four areas, averaging only 22 pixels (or 5 1/2 
cells) each, for testing purposes. V/ith such an extremely 
small sample, a difference of 4.6% is entirely 

insignificant. it represents the mi scl ass i f ication of only 
1 cell. Thus the coniferous results are Included here only 

t 

for completeness. Mo conclusion can be dravvn from them. 
Observation 4 

The classes which benefit the most by the new 
techniques are deciduous forest and forage. Once again 
these are the classes with the largest generalized 

variances, as can be seen from Table 4,4.2. The correlation 
coefficient For this data Is .81, which Is coincidentally 
the same as for Figure 4.2.7. 

The performance of the supervised mode is plotted for- 
five values of the annexation threshold, t (Figures 4.4.3 
and 4.4.4). The optimum performance occurs at the value t= 
4. At thi-s value, the overall error rate is reduced from 
12.9 percent to just 4,2 percent. For this study the 
Level-1 threshold was held at c= 60, so again the results 

r 

are'not necessarily optimized v/ith respect to c. 

The performances of the unsupervised modes are given in 
Tables 4.4.3 through 4.4.6 for various combinations of 
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Figure 4.4.4 Effect of Annexation Threshold (t) on Overall Performance 
- Run 71052501 
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Table 4.4.3 Matrix of Average Error Rates C%) for Four 

Combinations of Significance Levels (MV mode) 
- Run 71052501 



table 4.4.4 Matrix of Overall Error Rates for Four' 

Combinations of Significance Levels (MV mode) 
- Run 71052501 



* result of "cel 1 -spl i tt i ng’* 
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significance levels. The study was done using a constant 
Level-1 threshold of 0.45. This unusually large value is a 
reflection of the narrow autocorrelation functions of the 
forest class (Appendix A). Cell splitting was found to 
reduce the average error rate some/ but it has little or no 
effect on the overall error rate. The best results for the 
f'UV version occur at about 0 and Sj^= .100, while the MV 
version Is best at S 2 = 0 and s^= .050. 

Figure 4.4.5 shows how the processors compare with 
regard to both overall error rate and CPU time. The small 
number of classes and channels has reduced the time required 

t 

to do a single classification to such, a low point that a 
processor's speed depends on considerations other than the 
number of classifications that it must perform. Thus the 
speed advantage of the unsupervised modes has disappeared, 

4.5 RUN 72032803 - Classification of Satellite Data 

This LANDSAT-1 data set covers a 66.3 km long by 111 km 
wide region in Illinois on August 9/ 1972. The region was 
sampled 837 times lengthwise, and 1920 times along its width. 
The subregion analyzed is 40.3 km by 46,3 km and contains 
406/400 pixels. This area was chosen because it contains a 
large set of reference data (15/067 pixels)/ corresponding 
to actual surface observations |28[. The 5 main classes are 
corn/ sov/ other-agriculture (including alfalfa/ grasS/ 
oatS/ pasture/ and wheat)/ woodS/ and town. Training data 
was obtained by sampling the test areas as indicated in 
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Figure 4.4.5 Error Rate and CPU Times for Five Classification Schemes 
- Run 71052501 
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4.5 

Section 4.1. The resultant ratio of test data to training 
data is about 5:1. Examination of the training data 
indicated that no subclasses v>(ere needed. All four channels 
were used in the analysis, because the separability of the 

classes is fairly low. 

Figure 4.5.1 shows the average and overall performances 
achieved by the various classification schemes. The weights 
for overall error are given in Table 4.5.1. The 
class-conditional error rates appear in Figure 4.5.2. 
Several observations are v/orthy of mention: 

Observation 1 

Once again supervised ECHO greatly reduced the average 
and overall error rates and performed as well or much better 
than CLASSIFYPOINTS on every class. On three of the five 
classes, both the supervised and unsupervised modes 

out-performed even SAMPLECLASS I FY CML) . 

Observation 2 

The unsupervised mode did very well in terms of overall 
error rate, but Its average error rate Is slightly greater 
than that of CLASSIFYPOINTS due to unusually poor 

performance in the "woods" class. The reason for this 

behavior is unknown, but it suggests that the unsupervised 
mode tends to be less stable than the supervi sed mode. 

Observation 3 

This is the only data set in v/hich a "town" class has 
We would anticipate a fairly broad 


been included. 
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Processor Key 

#1 CLASS IFYPOINTS 
#2 Supervised Partitioning, t*=0 
#3 Optiraum Unsupervised Partitioning 
#4 Optiraum Supervised Partitioning . 
#5 SAMPLECLASSIFY (ML) 

#6 SAMPLECLASSIFY (MD) 


Figure 4.5.1 Classification Performance vs. Processing Scheme 
- Run 72032803 
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Table 4,5,1 Relative Influence of Each Class on Overall 
Performance - Run 72032803 

Percentage of 


Class 

Total Test 

Corn 

67,1 

Soy 

18. '8 

Other-Ag, 

10,2 

Town 

2,9 

Woods . 

0,9 


Table 4.5.2 Classes Ordered According to Generalized 
Variance ” Run 72032803 


Class 

Other-Ag. 

Town 

Soy 

Corn 


Common Logarithm of 
Generalised Variance 

4.38 

3.63 

3.50 

2,61 

2.27. 


Improvement (sup. mode) 
Over CLAS SIFYPOINTS 

+ 14.0 

+ 21,7 

+ 1,6 

+ 21,3 


'Woods 


0.0 




Figure 4.5.2 Performance By Class (Run 72032803) 
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Figure 4.5.2, continued 
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4.5 

distribution for this class and therefore a potentially 
large i mprovernent In classification accuracy by the ECHO 
technique. The results bear this out. Table 4.5.2 shows 
that the town class ranks second in terms of generalized 
variance and first in terms of classification improvement. 
But unlike the aircraft data sets., the overall .correlation 
between these quantities is fairly low. 

‘ The performance of the supervised mode is plotted' vs. t 
in -Figures 4.5.3 and 4.5.4. The best results occur at t=3 
and t='4. The Level -1 threshold was held constant at c=50 
for this study, and cel 1 -s pi i tt i ng was not used. The effect 
of c was briefly investigated for t=3 by trying the values 
c=30 and c=80. Only class "other” was significantly 
affected. For c=30 its error rate • i ncreased 5.4% and for 
c=80 It Increased just 0.3%. Thus performance again appears 
to be quite stable with respect to c. 

Table 4.5.3 shows the average and overall error rates 
of the unsupervised mode for seven combinations of parameter 
values. Cell splitting was not used. The best performance 
occurs for $2=0, Sj^ = .005 and c=.25 or .20. c=.25 produces 
the .lowest overall error, and c=.20 produces the lowest 
average error. The -value .25 was finally judged as "best" 
for this data for two reasons. (1) The value .20 produces 
an I excessive number of singular cells. The value .25 
produces a pattern of singular cells more closely resembling 
boundaries, as desired. C2) The high average error rate for 
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Figure 4.5.3 Effect of Annexation Threshold (t) on Average Performance 
- Run 72032803 
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Table 4.5.3 

Effect of Parameters on 
( un superv J sed mode) 

Performance 

^2 

c 

• H 

Avg. 

Error 

Overal 1 
Error 

.001 

.25 

.005 

40.3% 

38.8% 

0 

.25 

.001 

36,2 

35 .5 

0 

.25 

.005 

42.8 

29.7 

0 

.25 

.025 

41.7 

36.1 . 

0 

,20 

.001 

42.5 

32.3 

0 

.20 

,005 

35.2 

32.1 

0 

.15 

.005 

38.9 

35.2 
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c=*.25 is due solely to the behavior of the mi nor i ty' cl ass 
"v>/oods'‘. (All the other classes are at their optimum 
performances for this value.)' This behavior is believed to 
be an anomalous effect attributable mainly to the 'small 
amount of test data available for woods (.9?;). This belief 
is 'supported by the unusual observation that for si?. 025, 
.00,5, and .001, the error rates for this class are 56.9^, 
89. 2?^,, and 48.5% respectively, !\ similar effect also 
occurred in' the supervised mode. For t=2, 3, 4, and 5, the 
error rates are 14.6%', 28.5%; 31.5%, and 9.2% respectively. 

Figure 4.5.5 comoares the overall error and CPU times 
of the various processors. In this case the partitioning 
schemes require more CPU time ,than the non-partitioning 
schemes. This can be attributed to the lov/ number of 
spectral classes combined with a small number of pixels per 
physical object. 

We also note that processor #2 performed^ only 
one-fourth as many actual classifications as #1 (since the 
cell size v^as 4), and yet they required the' same CPU time. 
The same effect occurred on Run 71052501. This indicates a 
significant margin for Improvement of the efficiency (speed) 
of the supervised mode. The cause of the inefficiency is 
that the ECHO processors are coded in FORTRAN, whereas the 
classification subroutine - in CLASS ! FYPO I NTS is- coded in 
assembler language for optimum efficiency. 
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CPU Time (minutes) 

Figure 4.5.5 Error Rate and CPU Time for Four Classification Schemes 
- Run 72032803 



124 

4 .& 

4,6 General Observations 

■ Having studied four significantly different data sets, 
it is now possible to fnake' some general observations. 
Strictly speaking, our' observations (and therefore our 
conclusions) apply only to these data sets. But due to the 
consistency of the results, it is reasonable to expect-other 
data sets to exhibit similar behavior as long as they are 
reasonab 1 y s im i 1 ar to these. This .assumes, of course, that 
representat i.ve training data is available and that training 
is done in accordance with the assumptions of the model. 
Observation 1 

First of all it is apparent that the new techniques are 

effective for a wide range of variables including 

classification parameters (number of classes, subclasses, 

channels, etc.) and data ' parameters (spatial, spectral, and 

measurement -tesol ut Ion and spatial correlation). . 

Observation 2 f. 

Equally important is the stability of the perfoirmance 

with respect to the processor's input parameters. For the 

supervised inode the performance is very stable. The main 

input parameter Is the annexation threshold, t. On all four 

data sets all values tried gave better performance than 

CLASS I FYPOi MTS. Furthermore the value t=4 gave optimal, or 

near optimal, performance in every- case. There is- some 

/ 

evidence to suggest that as the number of channels 
increases, so does the optimum value of t. 
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For the onsupervised MUV mode the main input parameters 
are the significance levels Si and S 2 / in that order. The 
results seem to be fairly stable v/ith respect to them, but 
their optimum values tend to be elusive. For the four data 
sets the optimum value of Si varies among .001, .005, and 

.100. S 2 is usually effective at around .001 or 0. 
Presently' there is no way to predict the optimum 
combination, which means that in order to use this mode one 
will generally have to settle for suboptimal results. in 
most cases however, the results are still somewhat better 
than CLASS iFYPOl NTS. 

A secondary parameter is the cell selection threshold 

/ 

c. A suitable value can generally be obtained by processing 
a small subregion of the data set with different thresholds 
until one is. found which produces a pattern of singular 
cells that resembles object boundaries (as opposed to random 
noise). In the unsupervised mode c=,25 most frequently 
produced this effect, but when spatial correlation was low, 
a larger value was needed. In the supervised mode, c 

depends on the number of channels, q. As a rule of thumb, 
c= 15q appears to be a reasonable empirical guideline (at 
least for 3' ><: q x< 6). But a larger value (c= 20q) was 

required when the spatial correlation was low. 

Observation 3 

Another type of stability is the sensitivity of the 
processor to the particular characteristics of a class. 


For 
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example/^ CLASS I FYPO I'NTS Is sensitive to the variance of the 
class pdf. It tends to become "unstable" when a class with 
a relatively broad distribution Is encountered. Since a 
single pixel contains’ no Information about the variance of 
Its class,. CLASS I FYPO I NTS Is better suited for 
distinguishing classes that differ In mean rather than in 
var lance. 

The partitioning schemes, on the other hand, use both 
the mean and variance of the data In an object In order to 
classify It. Consequently they tend to be much more stable 
In this respect. But In return they should be somewhat 
sensitive to the typical size of the objects of a particular 
class. One would expect them to become unstable as the size 
of the objects approaches the cell size. Some evidence of 
this v/as detected (e.g. pit class), but the problem was 
avoided for the most part by choosing a relatively small 
cell width (2 pixels). If the object size approaches this 
value, then Premise A Is violated and CLASS I FYPO I NTS also 
becomes unstable, as observed for the pit class. 

Observation 4 

The main advantage of the unsupervised mode appears to 

be speed v;hen classification time, as opposed' to 

partitioning time. Is the limiting factor. In other words, 

when the classification Is relatively complex. Involving 

* 

many spectral classes and. many channels, one can save a 
significant amount of time by extracting and classifying 
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objects rather than classifying on a point-by-po! nt basis. 
This is particularly true when the number of classes is 
large, because the time required for unsupervised 
partitioning is independent of the number of classes. The 
largest number of spectral classes considered in this 
investigation is 17, for which the unsupervised mode 
required less than half the time that CLASS I FYPO I NTS 
required. It is not uncommon to see analysts using many 
more than this, a trend attr j butabl e' to the increasing size 
of areas classified (encompassing more classes) and to the 
increasing use. of unsupervised clustering techniques for 
subclass definition. 

It seems advisable at this point to inject a word of 
caution against the use of clustering beyond that which is 
necessary to achieve reasonably unimodal training classes. 
"Creating" more subclasses reduces the number of pixels 
available per subclass for training. Thus the chance that 
each subclass is representative of an actual spectral class 
in the data set Is reduced. Also the narrow distributions 
that result cause the classifier performance to be more 
sensitive to non-stat ionary class statistics over a large 
area-. Finally, as the number of subclasses increases, the 
effective size of the objects In the scene decreases, thus 
reducing the 'potent ial advantage of a sample classifier. 

Observation 5 

The MV version of the unsupervised mode is not generally 
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applicable to all situations because of the relation between 
cell size and number of channels. In one case in wh'ich it 
v;as applicable it performed erratically^ but Its overall 
performance was as good as the MUV version. Although' cases 
probably exist for which the MV version provides a 
significant advantage over the MUV version, none have yet 
been encountered. 

; Observation 6 

The supervised mode consistently provided the lowest 
average and overall error rates (excluding SAMPLECLASS I FY), 
and for complex classifications it is more efficient than 
CLASS I FYPO I NTS . (This is in spite of- the fact that It is 
programmed considerably less efficiently than 
CLASSI FYPOI MTS . ) For the four data sets studied, the 
average reduction in class-average and overall error rates 
from the rates provided by CLASS I FYPO I NTS is 7.1% and 9.6% 
respectively. For comparison, the corresponding values for 
SAMPLECLASSIFY (ML) are 10.8% and 9.9%. For SAMPLECLASS I FY 
(MD) they are 9.7% and 8.9%. Thus we conclude that the ML 
strategy is an effective means of sample classification for 
MSS data, and that supervised partitioning is an effective 
means of applying it when the true' part i t i on is unknown. 
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CHAPTER 5 
CONCLUSION 


5.1 Summary 

The general problem that we have investigated is the 
application of sample classification techniques to MSS data. 
The purpose of this is to incorporate "memory" into the 
classification process^ thereby improving performance and 
reducing the number of items that must be individually 
classified. To begin/ we modeled the objects in the scene 
as simple samples from multivariate normal populations. 
Then we motivated the investigation by showing that the. 
class If icati^ scheme which achieves the minimum error rate 
v/hen object boundaries are known/ is a particular type of 
sample classification CMAP), Furthermore/ the closely 
related strategy (GML) that we use in practice has an error 
rate that is upper bounded by a sum of exponentially 
decreasing functions of the sample size. 

When object boundaries are unknown, the technique that 
is commonly used is to classify each resolution element 
independently of the others Cno-memory classification). 
This of course is suboptimal, because spatially adjacent 
states of nature are usually strongly dependent. However an 
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optimal strategy CAppendix B) would require a complete 
statistical description of this dependence as well as a 
large number of classifications and a large amount of 
computation per classification. As an intermediate approach 
we chose to exploit the particular nature of the dependence 
by applying an image partitioning transformation to the data 
set prior to actual classification. In particular we 
focused our attention on what we call the conjunctive 
obj ect**seek i ng approach. The - basic algorithm that we 
implemented requires two “levels" of statistical tests that 
are applied in a logical sequence in order to "merge" 
adjacent elements of the scene that are spectrally similar. 
The likelihood ratio procedure led us to multivariate tests 
of first and second order statistics using criteria whose 
distribution functions are known under the null hypothesis. 
This enables us to relate the size of the test to the 
decision threshold. The power of the test depends on the 
alternative hypothesis. Due to the "dimensionality problem" 
the "mul t ipl 6“uni var late" versions of these tests are 
actually more useful in practice. - 

The likelihood ratio procedure can also be adapted to 
provide a supervised mode of operation. The test statistic 
is more complex, but it is multivariate and yet free of the 
dimensionality problem. Experimental results indicate that 
it adds a significant measure of stability to the 
processor's performance. Consequently it consistently 
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provided the lowest class-average and overall error rates, 
approaching those attained by direct sample classification 
of the test areas. Also, visual inspection of the results 
indicated that the classification map produced by this 
method is much closer to the "type-map” form that is usually 
desired than Is the map produced' by no-memory 

class i f i cat i on . 

In terms of efficiency, the partitioning schemes varied 
from better to worse than the no-memory classifier to which 
they were compared, depending on the complexity of the 
classification. But given comparable programming' 

efficiency, it appears that this balance would be shifted 
significantly In favor of the partitioning approach. 

5.2 Recommendations for Further Work 

Parameter Selection 

One would expect that the optimum values of input 
parameters such as si, S 2 , and t are statistically dependent 
in some complex way on factors such as class separability, 
spatial correlation, number of channels, average object 
size, and cell size. Therefore it could be beneficial to 
investigate the possible use of such information to predict 
the optimum input parameters. As we have noted, this 
appears to be needed more in the unsupervised mode than In 
the supervised mode. 

Use of Texture 

In this investigation we have restricted ourselves to 
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the assumption that pixels are class-conditional ly 
uncorrelated. However in Appendix A we ‘have shov/n not only 
that they are correlated but that the correlation function 
is class-dependent. We have also indicated how a sample 
classifier could be designed to efficiently exploit this 
dependence for Improved discrimination between classes. 
This effect could be investigated through direct' sample 
classification of test areaS/ similar to the investigation 
done by V/acker and Landgrebe. if the degree of improvement 
proves to be significant/ then it is likely that the 
performance of the partitioning schemes can also be 
significantly improved by redesigning them to exploit 
spatial correlation. 
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APPENDIX A 
SPATIAL CORRELATION 

In Chapter 2 we assumed/ for simplicity/ that pixels 

within the same object are statistically independent 

observations from some subclass population* Then the joint 

pdf of the pixels can be expressed in terms of the marginal 

pdf of a single pixel. A more general approach is to allow 

for some correlation to exist between each pixel (A/ say) 

and the other pixels in the object that lie within some 

“neighborhood" of X* Spatial correlation can be inherent in 

the object/ and it can also be induced by the scanner. For 

example/ a line scanner induces a certain amount of 

correlation along each scan line because Its bandwidth Is 

constrained to reduce detector noise. Also/ one commonly 

finds that adjacent pixels actually overlap. In LANDSAT-1 

data/ for example/ the overlap Is 29%. Due to the position 

invariance of these effects and to the homogeneity of the 

objects, it is reasonable to assume that the 

class-conditional/ joint pdf of the pixels in the 

neighborhood of ^ Is invariant with respect to the position 

of X within the object Cas we previously did for the 
% 

marginal pdf of alone). This of course neglects any 
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non-stat ionary effects of the scanner such as banding/ 
variable sun-angle, non-uniform sampling, and d.c. drift, 
much of which can, in’ principle, be corrected for in the 
data. Another reasonable assumption for most types of 
objects is "transpose symmetry"; i .e. the joint pdf of a 
spatial array of pixels is invariant with respect to the row 
or column transpose of the array, A possible consequence of 
this is discussed later. 

{ 

Assuming (as before) that the pixels in an object are 
> ^ 

jointly MVN, then all that is required to specify their pdf 

Is a mean vector and covariance function matrix Ci.e. 

interchannel covariance matrix as a function of 

displacement). if J((m,n) represents the data vector for the 

scene element whose spatial coordinates are m and n (line 

and column numbers), then we can estimate this matrix for a 

given object by computing the following average: 

, C(m,n) = 1. Z (j((m+k,n+l) -MCm,n)) (XCk,l ) -M(0,0))* A1 
~ N k,l 

= i Z lCm+k,n+l )X' (k,l ) - M(m, n )M' (0, 0) 

H k,l 

where 

M(m,n) = 1. Z j<Cm+k,n + l) 

N ^k,l 

and N is the number of terms in the summation. To measure 
"local" characteristics (e.g, for a particular spectral 
class), the summations include only pixels from a single 
object. To measure the characteristics of a larger region, 
the summations extend over that region. These measurements 
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t , 

have been made for' several data setS/ and some typical 
results are presented in Figures A1-A5 (aircraft data) and 
Figures A6-A10 (satellite data). The quantity Rjj(m/n) 
which appears in • these figures is the correlation 
coefficient/ which is related to the covariance by 


Rjj(m/n) = Cjj(m/n)/ v/C j | ( 0 , oT Cjj ( 0 / 0 ) A 2 

where C;.(m/n) is the i,j 1^ element of matrix C(m/n). This 

i J 

normalization provides easier comparison of the various 
functions. Local correlation was measured within the test 
areas to ensure that only pixels from the same object were 
used. The. results were averaged over all the test- areas for 
a given class to obtain the final estimate for that class. 
Regional correlation is generally greater than local 
correlation for a given displacement/ because it includes 
the effect of dependent states as wel 1 as the 

class-conditional (local) correlation effects. 

The measurements indicate that intraclass spatial 

correlation is a significant effect. Naturally the 

strongest correlation occurs between adjacent pixels/ and 
the effect diminishes rapidly to a fairly low level- The 
"knee" in the curve generally occurs at a displacement of 
about 2 pixels. An important point is illustrated by 
Figures A2 and A3, which compare the classes "deciduous 
forest" and "forage". We observe a definite class 

dependency for the spatial- correlation functions. The 


relatively narrow correlation functions for forest indicate 
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a broader spatial frequency spectrum/ thus a faster rate of 
change/ than the forage class. The impl.lcation Is that this 
characteristic can be measured and used by a sample 
classifier to distinguish these two classes. Thus the 
"texture” of a sample contains potentially useful 

information about its identity. This comes as no surprise/ 
but it is not always obvious how to exploit such information 
in 'a numerical 1 y-oriented pattern recognition system. When 
spectral information is also available and both are observed 
In -a multivariate measurement space/ even the human system 
may be unable to use all the information contained In the 
sample. We shall now briefly discuss how the classifiers 
described in Chapter 2 can be generalized to accomplish such 
classification. 

To design a true MAP or i4L classifier for a spatially 
correlated sample, X = CJlX' • • conceptually a 

straightforward matter, given the mean vector and covariance 
function matrix of each class. These can be used to 

construct the nq-dlmens ional matrix: 

C| = E(QfCX) IXeWj) l A3 

vjhere 

Qj(X) = ( (Xx-Mj )• / . . ./(il„-M5 )' )* (CXx-Mi)'^.../CXp-Mj)‘) 
bearing in mind the spatial arrangement of the pixels. Then 
for the hypothesis X e the log-1 i kel i hood function 

evaluated at X is given by 
p(XlXeWt) = -.5( lnl2nC|l + trCC-^Qj(X)) ) 


Ak 
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For large samples^ this rapidly becomes too cumbersome to be 
practical. Also, a covariance matrix would have to be 
constructed for each class for each unique spatial 
arrangement of pixels'to be classified. 

For large samples with relatively small correlation 
distances (small neighborhoods), a more practical approach 
would be to implement a minimum distance decision rule based 
on just the marginal joint pdf of a neighborhood, which 
contains all the information that characterizes a class 
population. In fact> assuming transpose symmetry holds, 
this Information is contained in the marginal, joint pdf of 
just one quadrant of the neighborhood. To estimate this pdf 
from a sample, we estimate the mean vector and covariance 
function matrix according to formula Al.' But we construct 
only the covariance matrix of an array of pixels 
corresponding to one quadrant of a neighborhood, thereby 
avoiding both problems associated with the maximum 
likelihood approach. 

We note that when assuming transpose symmetry, the- 
covariance matrix estimate can be improved for any given 
displacements, (m,n), by averaging CCm,n) and C(-m,n) to 

nJ isj 

obtain the final estimate. 

Also, the above strategy can be easily modified Cfor 
simplicity) by truncating the tail of the intraclass 
correlation function as desired, thereby reducing the 
dimension of the required covariance matrix. Since the knee 
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in the correlation function generally occurs at about 2 
pixels/ this might be a reasonable correlation distance to 
use. 

Finally, the choice of distance measure is arbitrary. 
However, an interesting possibility is the use of -LjCM/C) 
(from Section 2.3.4), appropr ! atel y modified for the higher 
dimensional space. We have seen that It is essentially the 
same as using a Kul 1 back- Lei bier number, but computationally 
it is much more efficient. Due to its relationship to the 
likelihood function and its lack of distance measure 
properties, a strategy that uses this criterion can perhaps 
best be described as a modified maximum likelihood strategy 
(modified to avoid an nq-dimens i onal matrix). 
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Figure A1 Spatial Correlation Coefficients - Run 71052501 
Region: lines 501-800, columns 1-100 
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DISPLACEMENT 

J 

Figure A2 Spatial Correlation Coefficients - Class FORAGE 
- Run 71052501 
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Figure A3 Spatial Correlation Coefficients - Class DECIDUOUS 
- Run 71052501 
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Figure A4 Spatial Interchatmel Correlation Coefficients 
- Class FORAGE - Run 71052501 


SPATIAL INTERCHANNEL CORRELATION COEFFICIENT 



Figure A5 Spatial Interchannel Correlation Coefficients 
- Class DECIDUOUS - Run 71052501 
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e®eoducbilw« 



Figure A6 Spatial Correlation Coefficients - Run 72032803 
Region: lines 260-559, columns 998-1097 
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Figure A7 Spatial Correlation Coefficients - Class CORN 
- Run 72032803 
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Figure A8 Spatial Correlation Coefficients - Class TOWN 
: - Run 72032803 
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Figure AlO Spatial Interchannel Correlation Coefficients - Class TOWN 
- Run 72032803 
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APPENDIX 8 • 

THE COMPOUND DECISION APPROACH 

Let denote a set of pixels (q-d imens ional 

random variables) with dependent states Cclasses), As in 
Section 2,1 we assume class-conditional independence. The 
goal is to classify each pixel into the set of classes W = 
Ci.e., classify x" into ) such that the 
expected number of mi scl ass i f i ed pixels is minimized. As 
shown In Section 2.3.2, this is accomplished by decision 
functions Wj(x^) such that 

P(Xi e Wf Cx'^) ix'^=x") = max P(l| e Wj IX" = X^), i=l,...,n B1 

J 

A simpler, though subopt imal procedure is to classify the 

pixels "sequentially" without "looking ahead"; I.e, let the 

i t h decision function be based only on the observations x* ~ 

( xj^ , , , ,,21i ) • This is equivalent to setting n=i in equation 

81. Making this change and applying Bayes rule provides: 

pCX' =x‘ IX; (x‘ ) ) PCWjCx')) = max p(X‘ =x MX; ^Wj ) PCWj ) B2 

j 

Before proceeding we shall shorten the notation by 

expressing the right hand side of B2 as; 

max p( X* 1 A; ) PC A| ) B3 

A; e W 

where A; denotes both a class and the event that X; is a 
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random variable from that class. Similarly A* will denote a 
vector of i classes and the event that these are the true 

‘ i 

classes of the set X- . . Defining j-i-l and applying the law 
of total- probability to B3 provides: 

max y p(x‘ } Aj ) PCAilA^) PCA-^ ) BA 

Aj ' , . , ■ ■ ■ . ■ 

A-’ e W-' • , . ■ ' 

= max 2 p(xi | A| ) PCA; lA-> ) pCxJ l AJ )/ P(AJ ) 
Ai , 

■ A-i 

Defining the quantity Qj CAi,# . , A| ) =. pCxMaV) P(A* it is 

apparent that 54 can be computed recursively as: 

max 2 Qj (Ai A| ) !• B5 

Ai j - ■ ■ ‘ ‘ 

A 

vihere 

Qi.CAi, .. .,A| ) = pCiii IA,j ) ,P(A| lA-^ ) QjXAi , . . , , Aj ) 86 

but the number of terms in. the summation grows as m"* , 

If the states' form a first order Markov chain, B4 
reduces to simply; 

. '■ 

max R| C A| ) ■ • ' •' ■ ' ' • B7 

Aj - • ■ ■ ■ 

where 


RjCAj) = pUilAj) 2 PCx-’iA'-') P(a') 

A-* 

but it can be computed recursively as; 

R;CA;),= p(iL||Aj)2 PCA||Aj) RjCAj) 

'A; 


B8 


. For a -kjLh order Markov chain C'k>l). we define j = I-k, and 
B4 reduces to; 

max p(Ai-iAf ) 2 P ^ A] - iIA{- 2_^2 ••••2 Rj(A{/».., A j + ]_ I 

A? A{_i ' ■ Aj- 2 Aj+i 
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v/here 

Ri(A|,. 2 PCxJ|A-J) PCAO 

A-* 


BIO 


P(A.|A. ./•••/A) R. 

I 1-1 j 1-1 1-1 j 

The number of summation terms in B9 and BIO Is only m^. 

To obtain this relatively efficient classifier we have 
assumed a sequential approach, Markov chain dependence 
between states, and cl ass-cond i t ionall y independent data 
vectors. In spite of these simplifications it is clear that 
the computational and memory requirements of the compound 
decision approach are considerably greater than the 
no-memory approach. 
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APPENDIX C 
COMPOSITE CLASSES 

Let p(x|W^) and p(xlW 2 ) be two uni variate normal class 
densities. Let be a third class that is a composite 
(spatially)’ of the two ground cover-types represented by 
these densities. The problem is to determine the density of 
this class. We shall assume that the conditional density of 
a pixel containing 100*a^ of class ' 1 and lOOCl-a)-^ of class 
2 is also normal with mean M(a) = a*M^^ + (l-a)M 2 and 
variance V(a)= a*V^ + (l-a)V 2 / where M. and V. represent the 
mean and variance of the ilh class. Thus 
’,.1 

p(x|Vl») = ! N(M(a)/ V(a) ;x) pCa) da Cl 

for some distribution p(a) which depends on the overall 
proportion of each class and the **texture“ of the composite. 

Assume for example that the overall proportion of each 
class is 0.5/ and consider three different cases of texture: 

Case 1 “ Maximum Variance 

If the texture of the composite is very coarse compared 
to the size of a pixel/ then a pixel will usually contain 
either one class or the other/ not both. pCa) is 
approximated by .’5( 6 (a) + 6 Ca’-l) )/ where gC*) represents the 
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Dirac delta function. Thus pCxjWj) is approximately 
.5 ( p(x lWi)+pCxlW 2 ) ) / which of course^ is bimodal In general. 
We - refer - to. this as a "hilxture" class/ because the 

constituents retain their individual characteristics. The 

normal way to handle such a class i s to cluster the data and 

treat- each mode as a subclass. Thus a mixture of W^- and W 2 

canfiot be treated as a d i st i net ‘ th i rd class by a Gaussian 
classifier. (Of course/ post-processing could be applied to 
the classifier output to search for such a mixture/’ if 
desjired.) 

j Case 2 - Minimum Variance 

If the texture of the composite is very fine/ then pCa) 
is approximated by (a- . 5 )/ ■ and . p(x|Wj) is approximately 
normal' with mean (Mj^+M2)/2 and variance (V2+V2)/2. This can 
(and should) be treated as a distinct third class. We> refer 
to [it as. a —compound*' clasS/ because the constituents lose 
their individuality. 

Case 3 - Intermediate Variance 

' When the texture of the composite is on the- same- order 
as the pixel size (e.g. Figures Cl and C2)/ a randomly 
selected pixel can contain any proportions of the two 
classes. As a first-order approximation to pCa) we consider 
the ' strictl y. periodic pattern in Figure Cl and let the 
coordinates (X/Y) of the pixel be uniformly distributed 
random variables over the area covered by the pattern. 
Assume for simplicity that the orientation of the pixel 
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square is the same as the square elements in the pattern. 

Symmetry considerations require that p(a) be unchanged if 

t 

the point (X,Y) is uniformly distributed over the region 0 
X ^ 1, 0 y min(x,l-x). In this region aCx,y) = x-2xy+y. 
By transformation of random variables it follows that: 
pCa) = - In 12a~ll , 0 4 a 4 1 C2 

which is shown in Figure Cl, The variance of this 

distribution is intermediate between Cases 1 and 2. 

Figure C3 shows two hypothetical class densities, 
pCxlW^) and p(x|W 2 ), and the density of their composite, 
p(xlW 3 ) (obtained by numerical evaluation of formula Cl 
using formula C2). As in Case 2, pCx[W 3 ) forms a distinct, 
unimodal class, but its variance is larger due to the spread 
of p(a) . 

Another interesting case is the pattern shown in Figure 
C2 . Here the texture is primarily one-dimensional. (This 
tendency can be observed for instance in a few of the 
objects in' Figure 4.2.4.) When the detector size matches 
the line width, p(a) is uniform on the interval 0 a 1. 
Its variance is greater than that of formula C2, which 
results in the composite class p(xlW 4 > shown In Figure C3 . 
Important points to notice are that it is unimodal, 
relatively broad, and it forms a distinct compound class. 
As increases relative to the variance of Wj and W 2 , 

p(xlW 4 ) tends toward uniform on the interval (Mi,M 2 ). Thus, 
if other classes lie in this interval, observations from the 
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compound class will tend to be mi scl ass i f 1 ed at a high rate 

I 

by a no-memory classifier. Of course/ much less confusion 

* I » 

should result if classification is done on a sample basis. 
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APPENDIX D 


FIELD LISTINGS 
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LAST 

LINF 

FIRST 

LINF 

LINE 

INT 

COL 

113 

124 

1 

20 

• 127 

131 

1 

34 

136 

139 

1 

22 

14? 

152 

1 

2 1 

227 

239 

1 

130 

2 33 

239 

1 

111 

242 

262 

1 

111 

2 42 

248 

1 

139 

254 

262 

1 

127 

265 

288 

1 

111 

291 

317 

1 

196 

292 

319 

1 

2 3 

292 

298 

1 

181 

301 

308 

1 

190 

311 

319 

1 

181 

322 

352 

1 

23 

324 

352 

1 

83 

364 

383 

1 

99 

365 

384 

1 

157 

386 

408 

1 

104 

388 

417 

1 

47 

389 

396 

1 

40 

42 5 

429 

1 

190 

432 

439 

1 

180 

442 

446 

1 

166 

452 

463 

1 

157 

453 

475 

1 

30 

453 

455 

1 

139 

458 

463 

1 

148 

466 

482 

1 

38 

466 

480 

1 

140 

486 

492 

1’ 

140 

486 

503 

1 

120 

495 

502 

1 

145 

50 5 

510 

1 

151 

515 

520 

1 

193 

517 

541 

1 

154 

519 

524 

1 

190 

522 

526 

1 

196 

524 

532 

1 

18i3 

52 9 

53? 

.1 

190 

536 

539 

1 

184 

538 

574 

1 

145 

539 

603 

1 

66 

579 

601 

1 

20 

609 

628 

1 

68 

613 

659 

1 

147 

634 

739 

1 

60 

66? 

729 

1 

149 

663 

729 

1 

165 

663 

697 

1 

183 

7 07 

739 

1 

100 

732 

740 

1 

149 

745 

801 

1 

45 

74 5 

761 

1 

68 

746 

800 

1 

29 

753 

761 

1 

84 

7 55 

765 

1 

96 

755 

807 

1 

121 

756 

792 

1 

141 

764 

784 

1 

86 

768 

772 

1 

9 5. 

768 

77? 

1 

11? 


LAST 

COL 

COVE 

COL 

INT 

TYPE 

45 

1 

CORN 

43 

1 

CORN 

29 

1 

CORN 

43 

1 

CORN 

140 

1 

CORN 

117 

1 

CORN 

124 

1 

CORN 

145 

1 

CORN 

136 

1 

CORN 

145 

1 

CORN 

207 

1 

CORN 

93 

1 

CORN 

193 

1 

CORN 

193 

1 

CORN 

184 

■ 1 

CORN 

54 

1 

CORN 

85 

1 

CORN 

116 

1 

CORN 

160 

. 1 

CORN 

136 

1 

CORN 

94 

1 

CORN 

44 

1 

CORN 

199 

1 

CORN 

204 

1 

CORN 

210 

1 

CORN 

178 

1 

CORN 

33 

1 

CORN 

142 

1 

CORN 

154 

1 

CORN 

42 

1 

CORN 

179 

1 

CORN 

145 

1 

CORN 

134 

1 

CORN 

149 

1 

CORN 

153 

1 

CORN 

196 

1 

CORN 

178 

I 

CORN 

192 

1 

CORN 

198 

1 

CORN 

186 

1 

CORN 

193 

1 

CORN 

186 

1 

CORN 

146 

1 

CORN 

69 

1 

CORN 

22 

' 1 

CORN 

79 

1 

CORN 

170 

1 

CORN 

97 

1 

CORN 

162 

1 

CORN 

180 

1 

CORN 

192 

1 

CORN 

114 

1 

CORN 

172 

1 

CORN 

61 

1 

CORN 

81 

1 

CORN 

40 

1 

CORN 

92 

1 

CORN 

1L5 

1 

CORN 

126 

1 

CORN 

145 

1 

CORN 

92 

1 

CORN 

104 

1 

CORN 

118 

1 

CORN 
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TABLE Dl» 
71052800 

CONTINUED 

0-6 

775 

784 

1 

97 

118 

1 

CORN 

71052800 

0-21 

7 86 

801 

1 

79 

83 

1 

CORN 

71052800 

0-6 

787 

800 

1- 

88 

118 

1 

CORN 

71052800 

0-6 

789 

799 

1 

129 

133 

1 

CORN 

71052800 

0-10 

803 

846 

1 

27 

37 

1 

CORN 

71052800 

0-21 

803 

823 

1 

68 

98 

1 

Corn 

71052800 

0-21 

804 

811 

1 

52 

65 

1 

CORN 

71052800 

0-19 

804 

817 

1 

142 

147 

1 

corn 

71052800 

0-10 

817 

846 

1 

40 

54 

1 

CORN 

710 52800 

0-19 

819 

823 

1 

133 

140 

1 

CORN 

71052800 

0-14 

839 

842 

1 

142 

147 

1 

CORN 

71052800 

0-10 

840 

846 

1 

74 

85 ■ 

1 

CORN 

71052800 

0-14 

84 8 

850 

1 

150 

1 54 

1 

CORN 

710 52800 

0-10 

849 

859 

1 

28 

94 

1 

CORN 

71052800 

0-14 

854 

858 

1 

159 

162 

1 • 

CORN 

7 1052800 

0-14 

861 

865 

1 

1.65 

168 

1 

CORN 

71052800 

0-18 

881 

891 

1 

104 

118 

1 

CORN 

71052800 

0-18 

894 

944 

1 

88 

129 

1 

CORN 

71052800 

0-17 

903 

913 

1 

136 

148 

1 

CORN 

71052800 

0-22 

905 

943 

1 

75 

83 

1 

CORN 

71052800 

0-22 

916 

943 

1 

70 

72 

1 

CORN 

710 52800 

0-17 

925 

943 

1 

134* 

172 

1 

CORN 

71052800 

0-17 

93 7 

941 

1 

173 

182 

1 

CORN 

71052 BOO 

R-4 

941 

944 

1 

28 

58 

1 

CORN 

71052800 

R-B 

949 

961 

1 

44 

.54 

1 

CORN 

71052800 

R-3 

950 

1005 

1 

31 

37 

1 

CORN 

71052800 

R-2 

950 

1005 

1 

40 

41 

1 

CORN 

710 52800 

R-8 

964 

973 . 

1 

44 

49 

1 

CORN 

71052800 

AA-2 ■ 

965 

969 

1 

156 

164 

1 

CORN 

7 10 52800 

AA-1 

976 

985 

1 

178 

.187 

1 

CORN 

71052800 

V-6 

989 

997 

1 

140 

144 

1 

CORN 

71052800 

V-i 

1008 

1023 

1 

173 

184 

1 

CORN 

71052800 

1-1 

1014 

1037 

1 

121 

1 34 

1 

CORN 

71052800 

S-15 

1041 

1049 

1 

37 

55 

1 

CORN 

71052800 

S-13 

1052 

1065 

1 

37 

57 

1 

CORN 

71052800 

S-13 

1053 

1066 

1 

63 

94 

1 

CORN 

71052800 

S-6 

1074 

1101 

1 

130 

143 

1 

CORN 

710 52800 

S-2 

1083 

1098 

1 

1 55 

184 

1 

CORN 

71052800 

S-6 

1105 

1121 

1 

139 

142 

1 

CORN 

71052800 

S-6 

1105 

1160 

1 

145 

146 

1 

CORN 

71052800 

S-4 

1120 

1134 

1 

166 

185 

1 

CO R N 

71052800 

S-6 

1126 

1165 

1 

152 

L57 

1 

CORN 

71052800 

S-10 

1156 

1198 

1 

106 

127 

1 

CORN 

71052800 

BB-1 

1204 

1217 

1 

95 

178 

1 

CORN 

71052800 

BB-4 

1219 

1231 

1 

184 

187 

1 

CORN 

7 10 52800 

BB-1 

1220 

1230 

1 

163 

172 

1 

C ORN 

71052800 

BB-1 

1226 

1240 

1 

73 

88 

1 

CORN 

71052800 

1-2 

112 

128 

1 

153 

181 

1 

SOY 

71052800 

J-1 

113 

127 

1 

96 

130 

1 

Soy 

71052800 

H-8 

207 

224 

1 

80 

89 

1 

SOY 

71052800 

H-2 

229 

288 

1 

19 

34 

1 

SOY 

71052800 

G-9 

262 

274 

1 

160 

175 

■ 1 

SOY 

71052800 

D-6 

2 63 

275 

1 

181 

208 

1 

SOY 

71052800 

G-9 

277 

284 

1 

171 

176 

1 

SOY 

7 10 52800 

G-8 

278 

288 

1 

149 

157 ■ 

1 

SOY 

71052800 

K-1 

32 3 

355 

1 

141 

163 

1 

SOY 

710 52800 

E-6 ■ 

482 

510 

1 

199 

209 

1 

SOY 

71052800 

E-3 

482 

488 

1 

172 

178 

1 

SOY 

71052800 

N-9 

511 

514 

1 

102 

109 

1 

SOY 

71052800 

CC-4 

532 

602 

1 

73 

76- 

1 

SOY 

71052800 

N-4 

581 

631 

1 

12 

16 

1 

SOY 

71052800 

X-2 

633 

714 

1 

117 

142 

1 

SOY 

710 52800 

X-5 

635 

708 

1 

28 

55 

1 

SOY 

71052800 

X-7 

712 

740 

1 

13 

55 

1 

SOY 

71052800 

X-2 

717 

722 

1 

136 

140 

1 

SOY 

71052800 

X-2 

717 

739 

1 

118 

126 

1 

SOY 
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I 




table DIt continued 


71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
710-52800 
71052800 
710-52 800 
710-52800 
71052800 
71052800 
71052-800 
710-52800 
710-52 800 
710(52 800 
71052800 
71052800 
71052800 
710,52‘800 
71052800 
710,52800 

710,52800 
71052800 
710,52800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
. 71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
71052800 
7i;0 52800 
71052800 

71052800 
7 LO 52 800 
710 52-800 


0-25 
M-4 
M-4 
M-4 
0-12 
R-5 
R-5 
R-5 
AA-7 
A A -4 
Z-4 
S-11 
S-11 
S-9 
S-9 
S-7 
S-7 
. S-9 
S-7 
S-5 
BB-3 
BB-3 

I-l 
H-19 
H-19 
H-19 
H-19 
H-19 
H-25 
H-25 
H-25 
H-1 0 
D-5 
G-7 
H-28 
H-5 
D-5- 
G-7 
G-5 
K-2 
T-2 
T-6 
U-1 
D-8 
F-3 
F-3 
E-5 
C-2 
A-1 
A-1 
0-6 
CC— 1 
0-3 
N-3 
\'}-2 
W-2 
0-8 
■M-1 
0-9 
M-1 
M-6 
M-6 
M-6 
0-13 
0-13 
0-13 


755 
7 87 
787 
805 
864 
867 

867 

868 
950 
950 

101 ? 

1041 

1061 

1073 

1073 

1077 

1093 

1136 

1138 

1139. 

1237 

1251 


784 

795 

802 

807 

874 

884 

896 

934 

1004 

977 

1033 

1056 

1067 

1080 

1091 

1084 

1104 

1153 

1157 

1196 

1245 

1263 


112 

113 

114 
120 
133 
.151 
153 
163 
174 
190 
228 
22 8 
231 
245 

249 

250 
292 
32 3 

358 

359 
373 
421 
437 
437 
483 

488 

489 

490 
547 
558 
589 
610 
613 
639 
745 
745 
765 
768 
822 
82 7 
829 
846 
8‘55 
864 


128 

172 

117 

130 

172 

172 

203 

171 

203 

203' 

246 

247 
288 
288 
259 
258 
319 
353 
363 
371 
383 
42 9 
447 
439 
510 
519 
504 
50 4 
566 
63.0 
6.08 
631 
628 
659 
751 
784 

783 

784 
846 
83 5 
843 
853 
8 62 
869 


1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 

1 

1 

1 . 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

. 1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 


130 
171 
150 
153 
133 
81 
67 
28 
15- 
129 
70 
127 
122 
49 
61 
110 
118 
99 
137 
175 
127 
112 

139 
75 
69 ■ 
64 

.57 

50 

132 

125 

116 

34 

181 

149 
78 
90 

180 

155 

140 
168 

150 
181 
201 

50 

139 

155 

191 

47 

15 

8 

159 
37 

167 
20 

174 
174 
96 

168 
78 

174 
151 
139 
164 
-138 
147 
157 


132 
181 
168 
163 
145 

91 

78 
64 

112 

133 
96 

137 

127 

58 

79 
115 
122 
103 

139 
187 

140 
14.5 

148 
89 
73 
73 
72 
54 
185 
129 
129 
41 
208 
176 
85 
105 
208 
176 
175 
192 
161 
198 
205 
55 
152 
160 
197 
52 
29 
12 
173 
44 
■181 
22 
181 
■ 181 
124 
171 
82 
180 
161 
148 
169 
146 
154 
161 


1 SOY- 
1 SOY 
SOY 
SOY 
SOY 


1 

1 

1 


1 SOY' 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


SOY 

SOY 

SOY 

SOY- 

SOY , 

SOY 

SOY 

SOY 

SOY 

SOY» 

SOY 

SOY.; 

SOY 

SOY 

SOY. 

SOY- 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


WHEAT 
WHEA T 
WHEAT 
WHEAT 
. WHEAT 
1 WHEAT 
W HEAT 
WHEAT 
WHEAT 
WHEAT 
1 WHEAT 
1 WHEAT 
WHEAT 
WHEAT 
WHEAT 
WHEAT 
WHEAT 
WHEAT 
. WHEAT 
1 WHEAT 
■ WHEAT 
WHEAT 
WHEAT 
WHEAT 
WHEAT 
WHEAT 
W HEAT 
WHEAT 
WHEAT 
WHEAT 
WHEAT 
WHEAT 
1 WHEAT 
1 WHEAT 
WHEAT 
WHFA T 
W HEAT 
WHEAT 
WHEAT 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


i W H £,A T 

. • < ■ It-:: A T 


WHEAT 
WHEAT 
W H.E A T 
WHEAT 
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TABLE Dl» 
71052800 

CONTINUED 
A A.- 3 

948 

961 

1 

146 

178 

1 

WHEAT 

71052800 

AA-3 

949 

960 

1 

181 

1 86 

1 

WHEAT 

71052R00 

M-2 

1018 

1030 

1- 

160 

166 

1 

WHEAT 

71052800 

S-3 

1103 

1115 

1 

154 

184 

1 

WHEAT 

71052800 

D-7 

279 

285 

1 

181 

207 

1 

IDLE 

71052800 

K-5 

322 

352 

1 

202 

208 

1 

IDLE 

710 52800 

D-3 

514 

520 

1 

182 

185 

1 

IDLE 

71052800 

CC-2 

551 

629 

1 

49 

62 

1 

IDLE 

71052800 

X-11 

765 

801 

1 

65 

74 

1 

IDLE 

71052800 

0-5 . 

816 

823 

1 

102 

121 

1 

IDLE 

71052800 

0-5 

826 

836 

1 

87 

118 

-1 

IDLE 

71052800 

0-11 

853 

870 

1 

115 

126 

■ 1 

IDLE 

71052800 

0-15 

878 

881 

1 

127 

132 

1 

IDLE 

71052800 

0-15 

884 

888 

1 

135 

140 

1 

IDLE 

71052800 

0-15 

898 

900 

1 

151 

154 

1 

IDLE 

71052800 

AA-6 

950 

1003 

1 

117 

125 

1 

IDLE 

71052800 

AA-8 

962 

1006 

1 

66 

69 

1 

IDLE 

71052800 

AA-8 

97 2 

979 

1 

60 

63 

1 

IDLE 

710 52 800 

R-7 

983 

1006 

1- 

56 

63 

1 

IDLE 

71052800 

S-8 

1073 

1088 

1 

87 

90 

1 

IDLE 

■71052800 

S-8 

1073 

1105 

1 

93 

101 

1 

IDLE 

71052800 

S-8 

1106 

1107 

1 

97 

100 

1 

IDLE 

71052800 

S-8 

1124 

1147 

1 

114 

120 

1 

IDLE 

71052800 

S-8 

1126 

1138 

1 

108 

1.11 

1 

IDLE 

71052800 

G-16 

32 4 

352 

1 

77 

80 

1 

HAY 

71052800 

G-I8 

340 

343 

1 

57 

73 

1 

HAY 

71052800 

F-5 

40 6 

411 

1 

192 

199 

1 

HAY 

71052800 

N-6 

497 

501 

1 

78 

85 

1 

HAY 

71052800 

N-8 

565 

603 

1 

81 

82 

1 

HAY 

71052800 

X-3 

634 

703 

1 

101 

112 

1 

HAY 

71052800 

X-3 

72 6 

740 

1' 

142 

145 

1 

HAY 

71052800 

H-20 

147 

152 

1 

94 

107 

1 

LESPEDEZA 

71052800 

H-20 

155 

159 

1 

94 

100 

1 

LESPEDEZA 

71 052800 

H-26 

207 

223 

1 

118 

186 

1 

LESPEDEZA 

7 1052800 

H-3 

229 

288 

1 

38 

46 

1 

L ES PEDF ZA 

71052800 

M-2 

746 

■ 784 

1 

150 

164 

1 

L ES PEDEZA 

71052800 

M-8 - 

822 

827 

1 

173 

180 

1 

LESPEDEZA 

71'0 52 800 

M-8 

830 

83,2 

1 

179 

184 

1 

LESPEDEZA 

71052800 

M-7 

852 

855 

1 

166 

171 

1 

L ES PED E Z A 

71052800 

H-22 

183 

202 

1 

■ 93 

96 

1 

PASTURE 

71052800 

H-12 

188 

206 

1 

5 

10 

1 

PASTURE 

71052800 

H-22 • 

205 

22 3 

1 

93 

100 

1 

PASTURE 

71052800 

H-6 

228 

240 

1 

97 

106 

1 

PASTURE 

71052800 

H-27 

228 

288 

1 

■ 67 

73 

1 

PASTURE 

71052800 

D-9 

368 

384 

1 

23 

26 

1 

PASTURE 

71052800 

C-3 

488 

507 

1 

56 

73 

1 

PASTURE 

71052800 

C-3 

510 

512 

1 

56 

64 

1 

PASTURE 

71052800 

E-2 

484 

491 

1 

149 

169 

1 

RYE 

71052800 

E-2 

5 04 

511 

1 

157 

169 

1 

RYE 

71052800 

D-4 

518 

52 6 

1 

206 

210 

1 

RYE 

71052800 

V-5 

985 

995 

1. 

148 

156 

1 

RYE 

71052800 

H-23 

133 

150 

1 

123 

127 

1 

WOODPAST 

71052800 

• H-23 

155 

162 

1 

112 

118 

1 

WOODPAST 

71052800 

H-23 

174 

202 

1 

101 

111 

1 

WOODPAST 

71052800 

H-2 3 

205 

22? 

1 

104 

113 

1 

WOODPAST 

71052800 

X-1 

719 

740 

1 

130 

■133 

1 

FOREST 

71052800 

R-6 

862 

866 

1 

99 

107 

1 

FOREST 

71052800 

V-7 

97 5 

977 

1 

148 

155 

1 

FORES T 

71052800 

V-7 

978 

980 

1 

139 

145 

1 

FOREST 

71052800 

V-4 

1018 

1022 

1 

151 

154 

1 

FOR ES T 
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TABLF DIt 
71052800 

CONTINUED 

S-14 

1021 

1032 

1 

58 

■61 

1 

FOREST - 

71052800 

Z-2 

10 30 

1039 

1 

104 

115 

1 

FOREST 

71052800 

H-16 

156 

158 

1 

21 

37 

1 

NON-FARM 

71052800 

H-17 

161 

172 

1 

43 

46 

1 

NON-FARM 

71052800 

H-21 

166 

173 

1 

93 

95 

1 

NON-FARM 

71052800 

H-11 

176 

181 

1 

31 

35 

1 

MON- FARM 

7 1052800 

H-11 

184 

189 

1 

28 

29 

1 

NON-FARM 

71052800 

H-11 

192 

201 

1 

18 

22 

1 

NON-FARM 

71052800* 

H-7 ' 

227 

235 

1 

89 

91 

1 

NON- FARM 

71052800 

G-12 

252 

258 

1 

146 

151 

1 

NON- FARM 

71052800 

G-13 

347 

353 

1 

71 

'73 

1 

NON-FARM 

71052800 

G-1 

349 

354 

1 

97 

100 

1 

NON-FARM 

71052800 

G-1 

355 

359 

] 

100 

115 

1 

MON- FARM 

71052800 

P-2 

422 

42 8 

1 

22 

44 

1 

NON-FARM 

7 1052800 

P-2 

431 

449 

1 

37 

46 

1 

NON-FARM 

71052800 

P-2 

454 

473 

1 

48 

56 

1 

MON- FARM 

71052800 

P-4 

478 

483 

1 

31 

34 

1 

NON-FARM 

71052800 

B-1 

50 8 

561 

1 

9 

29 

1 

NON-FARM 

71052800 

B-1 

516 

525 

.1 

67 

73 

1 

NON- FARM 

71052800 

B-1 

52 3 

535 

1 

55 

64 

1 

NON-FARM 

71052800 

8-1 

5 32 

5 39 

1 

45 

52 

1 

NON- FARM 

71052800 


540 

541 

1 

47 

52 

1 

NON-FARM 

71052800 

e -1 

540 

545 

1 

33 

46 

1 

NON- FARM 

71052800 

0-7 

546 

566 

1 

178 

180 

1 

NON-FARM 

71052800 

0-5 

548 

574 

1 

151 

153 

1 

NON-FARM 

71052800 

B-1 

563 

631 

1 

27 

32 

1 

NON- FARM 

71052800 

B-1 

5 64 

568 

1 

9 

22 

1 

NON-FARM 

71052800 

0-1 

569 

576 

1 

178 

182 

1 

NON- FARM 

71052800 

0-2 

579 

586 

1 

178 

182 

1 

MON- FARM 

71052800 

V^-1 

631 

636 

-] 

179 

181 

1 

NON-FARM 

71052800 

0-1 

702 

709 

1 

• 184 

190 

1 

NON- FARM 

71052800 

0-10 

733 

739 

1 

176 

183 

1 

NON-FARM 

71052800 

0-7 

744 

751 

1 

128 

133 

1 

NON-FARM 

71052800 

X-8 

744 

749 

1 

85 

92 

1 

NON-FARM 

71052800 

0-4 

8 37 

845 

1 

121 

126 

1 

NON-FARM 

71052800 

DD-1 

990 

995 

1 

193 

194 

1 

NON-FARM 

71052800 

DD-1 

998 

1004 

1 

194 

202 

1 

NON-FARM 

71052800 

S-1 

1072 

1077 

1 

151 

166 

1 

NON-FARM 
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TABLE D2 

TEST AREAS 

FOR RUN 72064412 




RUN 

AREA 

FIRST 

LAST 

LINE 

FIRST 

LAST 

COL COVER 

NUMBER DESIGNATION 

l^ne 

LINE 

INT 

COL 

COL 

INT TYPE 

7206441? 

AG 

59 

62 

1 

, 30.7 

317 

1 

AGRICULTURE 

72064412 

AG 

68 

77 

1 

426 

432 , 

1 

AGRICULTURE 

72064412 

AG 

74 

78 

1 

369 

375 

1 

AGR IC ULTURE 

72064412 

AG 

79 

87 

.1 

2 62 

265 

1 

AGRICULTURE 

72064412 

AG 

84 

90 

1 

273 

285 

1 

AGRICULTURE 

7206A412 • 

AG 

100 

110 

1 

274 

283 

1 

AGRICULTURE 

72064412 

AG 

106 

124 

1 

234 

244 

1 

AGRICULTURE 

72064412 

AG 

111 

114 

1 

245 

257 

1 

AGR I CULTURE 

7206441? 

AG 

113 

119 

1 

276 

288 

1 

AGRICULTURE 

72064412 

AG 

12 6 

130 

1 

252 

267 

1 

AG RICULTURE 

72064412 

AG 

164 

174 

1 

586 

591 

1 

AGRICULTURE 

72064412 

AG 

180 

187 

1 

369 

373 

1 

AGRICULTURE 

72064412 

AG 

194 

199 

1 

264 

267 

1 

AGRICULTURE 

72064412 

AG 

197 

205 

1 

233 

238 

1 

AGRICULTURE 

72064412 

AG 

202 

209 

1 

247 

253 

1 

AGRICULTURE 

72064412 

AG 

212 

217 

1 

296 

302 

1 

AGRICULTURE 

720644 12 

AG 

217 

221 

1 

275 

280 

1 

AGR IC ULTURE 

72064412 

AG 

219 

222 

1 

■443 

454 

1 

AGRICULTURE 

72064412 

AG 

229 

232 

1 

395 

403 

1 

AGRICULTURE 

72064412 

AG 

247 

251 

1 

268 

274 

1 

AGRICULTURE 

72064412 

AG 

250 

253 

1 

488 

497 

1 

AGRICULTURE 

72064412 

AG 

257 

264 

1 

236 

240 

1 

AGR I CUL TURE 

72064412 

AG 

257 

263 

1 

291 

298 

1 

AGR IC ULTURE 

72064412 

AG 

283 

289 

1 

501 

506 

1 

AGRICULTURE 

7206441? 

AG 

300 

302 

1 

514 

521 

1 

AGRICULTURE 

72064412 

AG 

30 8 

315 

1 

314 

319 

1 

AGR ICULTURE 

72064412 

AG 

319 

326 

1 

302 

306 

1 

AGRICULTURE 

72064412 

AG 

32 0 

324 

1 

366 

369 

1 

AGRICULTURE 

720644 12 

AG 

350 

354 

1 

233 

238 

1 

AGRICULTURE 

72064412 

AG - 

376 

380 

1 

400 

404 

1 

AGRICULTURE 

72064412 

AG 

3 84 

388 

1 

524 

530 

1 

AGRICULTURE 

72064412 

AG 

400 

404 

1 

538 

542 

1 

AGRICULTURE 

72064412 

AG 

406 

411 

1 

353 

356 

1 

AGR ICULTURE 

72064412 

AG 

416 

418 

1 

372 

375 

1 

AGRICULTURE 

72064412 

AG — 

422 

424 

1 

■ 524 

536 

1 

AGR ICULTU RE 

72064412 

AG 

430 

433 

1 

385 

394 

1 

AGRICULTIJRE 

72064412 

AG 

460 

464 

1 

265 

274 

1 

AGR I CU LTU RE 

720'64412 

AG 

460 

467 

1 

495 

500 

1 

AGR I CULTURE 

7206441? 

AG 

461 

464 

1 

293 

298 

1 

AGR ICULTURE 

72064412 

AG 

464 

470 

1 

320 

327 

1 

AGRICULTURE 

72064412 

AG 

479 

482 

1 

452 

455 

1 

AGRICULTURE 

72064412 

AG“ 

488 

491 

1 

397 

402 

1 

AGRICUL TURE 

72064412 

AG 

499 

503 

1 

487 

492 

1 

AGRICULTURE 

72064412 

AG 

50 8 

510 

1 

328 

333 

1 

AGRICULTURE 

7206441? 

AG 

541 

545 

1 

483 

490 

1 

AGRICULTURE 

72064412 

AG 

563 

567 

1 

431 

440 

1 

AGRICULTURE 

72064412 

AG 

573 

577 

1 

352 

358 

1 

AGR ICULTU RE 

72064412 

AG 

57 8 

585 

1 

326 

331 

1 

AGRICULTURE 

720 644 1? 

AG 

.580 

5 86- 

1 

348 

352 

1 

AGRICULTURE 

72064412 

AG 

582 

588 

1 

438 

443 

1 

AGRICULTURE 

72064412 

AG 

586 

591 

1 

510 

515* 

1 

AGR ICULTURE 

720644 12 

F-S 

66 

68 

1 

236 

247 

1 

FOREST 

72064412 

F-T 

85 

96 

1 

367 

403 

1 

FOREST 

72064412 

F-T 

86 

91 

1 

322 

326 

1 

FOREST 

72064412 

F-T 

89 

98 

1 

313 

321 

1 

FORES T 

72064412 

,F-T 

89 

101 

1 

424 

434 

1 

FOREST 

72064412 

F-T 

90 

103 

1 

330 

357 

1 

FOREST 

72064412 

F-T 

92 

101 

1 

2 66 

270 

1 

FOREST 

72064412 

F-T 

94 

104 

1 

< 232 

236 

1 

foreS t 

72064412 

F-T 

94 

98 

1 

275 

283 

1 

FOREST 

72064412 

F-T 

97 

107 

1 

361 

375 

1 

FO RES T 

72064412 

. F-T 

106 

112 

1 

489 

499 

1 

FOREST 
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TARLF D2, 

CO NT I MU ED 

7 20-64412 

F-S 

72064412 

F-T 

72064412 

F-S 

72064412 

F-T 

7 20^64412 

F-T 

72064412 

F 

72064412 

F 

7 2 0644 12 

F-T 

72064412 

F 

72064412 

F-T 

72064412 

F 

72064412 

F 

72064412 

F-S 

72064412 

F-T 

72064412 

F-T 

72064412 

F-T 

72064412 

F-S 

72064412 

F-T 

72064412 

F-S 

72064412 

F-T 

72064412 

RM 

72064412 

RM 

72Q64412 

PRVGM 

72064412 

RM . 

72064412 

RM 

72064412 

PRVGM 

72064412 

PRVGM 

72064412 

RM 

72064412 

PRVGM 

72064412 

PRVGM 

72064412 

RM 

72064412 

PRVGM 

72064412 

RM 

7 20,64412 

PRVGM 

72064412 

PRVGM 

72064412 

PRVGM 

72064412 

PRVGM 

72064412 

PRVGM 

72064412 

PRVGM 

72064412 

PRVGM 

72064412 

PRVGM 

72064412 

PRVGM 

72064412 

RM 

72064412 

RM 

72064412 

RM 

72064412 

RM 

72064412 

RM 

72064412 

RM 

72064412 

RM 

72064412 

RM 

72064412 

PRVGM 

72064412 

PRVGM 

72064412 

PRVGM 

7206441? 

TIP 

72064412 

w-s 

72064412 

w— s 

72064412 

TIP 

72064412 

W-S 

72064412 

W— S 

72064412 

w-s 

72064412 

W“C 

72064412 

w-c 

72064412 

w-c 

72064412 

TIP 

72064412 

TIP 


132 

136 

1 

340 

132 

156 

1 

5 5 7 

139 

147 

1 

48 6 

157 

164 

1 

557 

170 

181 

1 

542 

202 

205 

1 

379 

213 

221 

1 

380 

254 

257 

1 

270 

321 

335 

1 

579 

33 9 

348 

1 

320 

379 

388 

1 

572 

412 

42 9 

1 

572 

436 

443 

1 

30 5 

456 

469 

1 

236 

456 

469 

1 

584 

460 

467 

1 

509 

514 

522 

1 

384 

547 

551 

1 

433 

556 

565 

1 

268 

581 

588 

1 

'477 

57 

72 

1 

484 

68 

72 

1 

436 

81 

87 

1 

496 

153 

157 

1 

354 

166 

174 

1 

382 

170 

176 

1 

518 

172 

174 

1 

317 

175 

179 

1 

388 

181 

189 

1 

498 

185 

188 

1 

224 

189 

191 

1 

327 

192 

194 

1 

.230 

236 

242 

1 

325 

331 

335 

1 

453 

339 

346 

1 

460 

362 

369 

1 

302 

366 

373 

1 

315 

378 

383 

1 

343 

385 

389 

1 

333 

404 

413 

1 

267 

404 

408 

1 

30? 

407 

415 

1 

254 

41? 

415 

1 

289 

414 

416 

1 

279 

42 3 

42 8 

1 

316 

440 

445 

1 

323 

44 5 

465 

1 

388 

445 

447 

1 

515 

466 

469 

1 

385 

480 

485 

1 

294 

492 

49 4 

1 

297 

497 

500 

1 

299 

506 

521 

1 

306 

58 

62 

1 

465 

62 

63 

1 

519 

64 

65 

1 

534 

76 

79 

1 

32 9 

80 

80 

1 

407 

80 

80 

1 

534 

81 

81 

1 

409 

151 

154 

1 

298 

153 

156 

1 

516 

186 

190 

1 

315 

236 

239 

1 

416 

2 40 

243 

1 

410 


1 


360 

1 

FOREST 

577 

1 

FOREST 

493 

1 

FOREST 

570 

1 

FOREST 

553 

1 

FOREST 

386 

1 

FOREST 

385 

1 

FOREST 

275 

1 

FOREST 

598 

1 

FOREST 

325 

1 

FOREST 

578 

1 

FOREST 

591 

1 

FOREST 

310 

1 

FOP.EST 

240 

1 

FOREST 

590 

1 

FOREST 

514 

1 

FOREST 

391 

1 

FOREST 

447 

1 

FOR'EST 

285 

1 

FOREST 

480 

1 

FOR.es T 

499 

1 

RECENT-MINE 

439 

1 

RECENT-MINE 

501 

1 

RECENT-MINE 

359 

1 

RECENT-MINE 

391 

1 

RECENT-MINE 

525 

1 

RECENT-MINE 

324 

1 

RECENT-MIME 

394 

1 

RECENT-MINE 

505 

1 

RECENT-MINE 

229 

1 

RECENT-MINE 

334 

1 

RECENT-MINE 

236 

1 

RECENT-MINE 

329 

1 

RECENT-MINE 

459 

1 

RECENT-MINE 

469 

1 

RECENT-MIME 

309 

1 

RECENT-MINE 

321 

1 

RECENT-MINE 

349 

1 

RECENT-MINE 

337 

1 

RECENT-MINE 

270 

1 

RECENT-MINE 

315 

1 

RFCENT-MINE 

259 

1 

RECENT-MINE 

293 

1 

recent-mine 

285 

1 

RECENT-MINE 

320 

1 

RECENT-MINE 

330 

1 

RECENT-MINE 

393 

1 

RECENT-MINE 

521 

1 

RECENT-MINE 

392 

1 

RECENT-MINE 

313 

1 

RECENT-MINE 

315 

1 

RECENT-MINE 

310 

1 

RECENT-MINE 

311 

1 

RECENT-MIME 

467 

1 

pit! 

520 

1 

PIT 

535 

1 

PIT 

332 

1 

PIT, 

408 

1 

PIT- 

536 

1 

PIT 

410 

1 

PIT. 

299 

1 

PIT. 

518 

1 

PIT 

320 

1 

PIT 

418 

1 

PIT 

413 

1 

PIT' 
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TABLE D2t 

CONT INDED 

72064412 

TIP 

72064412 

TIP 

72064412 

TIP 

72064412 

W-C 

72064412 

TIP 

72064412 

W-C 

72064412 

TIP 

72064412 

TIP 

72064412 

V'J-C 

72064412 

W-C 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM' 

72064412 

•RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

7206441? 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 

72064412 

RVGM 


245 

249 

1 

416 

24 8 

249 

1 

315 

251, 

2 52 

1 

324 

251 

256 

1 

447 

390 

396 

1 

288 

391 

393 

1 

320 

392 

397 

1 

292 

•394 

396 

1 

303 

469 

471 

1 

419 

565 

568 

1 

408 


78 

1 

391 

, 

76 

79 

1 

383 

115 

127 

1 

354 

149 

152 

1 

488 

150 

159 

1 

467 

161 

166 

1 

• 521 

164 

170 

1 

471 

176 

181 

1 

444 

192 

196 

1 

438 

2 09, 

215 

1 

52 8 

217 

220 

1 

480 

221 

225 

1 

342 

22 4 

228 

1 

323 

240 

244 

1 

349 

263 

266 

1 

405 

453 

455 

1 

465 

459 

461 

1 

338 

464 

467 

1 

345 

464 

467 

1 

458 

467 

471 

1 

464 

532 

538 

1 

345 

536 

548 

1 

357 

543 

550 

1 

302 

549 

553 

1 

331 

55 6 

559 

1 . 

321 

5 65 

567 

1 

312 

570 

573 

1 

264 

575 

579 

1 

308 


417 

1 

PIT 

320 

1 

PIT 

326 

1 

PIT 

448 

1 

PIT 

291 

1 

PIT 

322 

1 

PIT 

302 

1 

PIT 

308 

1 

PIT 

420 

1 

PIT 

411 

1 

PIT 

398 

1 

REVEG. MINE 

388 

1 

REVEG. MINE 

355 

1 

REVEG. MINE 

496’ 

1 

REVEG. MINE 

470 

1 

REVEG. MINE 

526 

1 

REVEG. MINE 

477 

1 

REVEG. MIME 

453 

1 

REVEG. MINE 

44 5 

1 

REVEG. MINE 

531 

1 

REVEG. MINE 

487 

1 

REVEG. MINE 

349 

1 

REVEG. MINE 

333 

1 

REVEG. MINE 

352 

1 

REVEG. MINE 

417 

1 

REVEG. MIME 

470 

1 

REVEG. MINE 

343 

1 

REVEG. MINE 

348 

1 

REVEG. MINE 

462 

1 

REVEG. MINE 

467 

1 

REVEG. MINE 

351 

1 

REVEG. MINE 

362 

1 

REVEG. MINE 

303 

1 

REVEG. MINE 

336 

1 

REVEG. MINE 

328 

1 

REVEG. MINE 

318 

1 

REVEG. MINE 

282 

1 

REVEG. MINE 

321 

1 

REVEG. MINE 
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TARLF D3 

REFFRFNCF 

AREAS 

FOR RUN 71052501 



RliW 

AREA 

FIRST 

LAST 

LINE 

FIRST 

LAST 

COL AREA 

NUMRFR DESIGMATinN 

LINF 

LINF 

I NT 

COL 

COL 

IMT TYRE 

7 10'5?501 

TRAINING 

173 

188 

1 

186 

192 

1 

DEC IDlinUS 

7105P501 

TRA IN IMG 

519 

52 5 

1 

026 

032 

1 

DECinunus 

71052501 

TRAINING 

046 

049 

1 

067 

070 

1 

COM I FERGUS 

71052501 

TRAINING 

048 

048 

1 

076 

078 

1 

C UN i F b KUDS 

71052501 

TRA INING 

052 

053 

1 , 

066 

067 

1 

CON i F E ROUS 

7105? 501 

TRAINING 

055 

055 

1 

073 

077 

1 

C UN 1 F b KUUS 

71052501 

TRA INING 

068 

072 

1 

052 

052 

1 

CUN 1 FbKDUS 

71052501 

TRAINING 

333 

339 

1 

083 

088 • 

1 

C ON 1 FE RUUS 

71052501 

TRAINING 

396' 

399 

1 

109 

126 

1 

CUNl FbKSJlJS 

7105250 1 

TRAINING 

146 

148 

1 

108 

112 

1 

WATER, 'pond 

7 10 52501 

TRAINING 

151 

155 

1 

113 

119 

1 

WATER , POND 

71052501 

TRAINING 

155 

159 

1 

120 

129 

1 

WA 1 fcK , PI.JND 

71052501 

TRAINING 

155 

156 

1 

023 

024 

1 

WA 1 bR , PUl'ID 

71052 501 

TRAINING 

227 

228 

1 

052 

053 

1 

WA 1 bK , P UNl) 

7 1052501 

TRAINING 

• 337 

338 

1 

108 

111 

1 

WATER , POND 

71052501 

TRA INING 

831 

834 

1 

Oil 

014 

1 

WA 1 bK , P UND 

7 1052’501 

TRAINING 

843. 

846 

1 

004 

008 

1 

WATER , POND 

71052501 

TRAINING 

847 

849 

1 

00? 

005 

1 

WATER , POND 

71Q52501 

TRAINING 

877 

879 

1 

041 

045 

1 

WATER , POND 

71052501 

TRA IN ING 

1174 

1174 

1. 

197 

202 

1 

WA 1 b K , K 1 V^b K 

710 52.501 

TRAINING 

1175 

1175 

1 

Oil 

020 

1 

WA 1 bKjKiVbK 

71052501 

TRA INING 

1184 

1185 

1 

030 

032 

1 

WATER T R 1 VbK 

71052501 

TRAINING 

1216 

1217 

1 

068 

073 

1 

WAI bR,KlVbK 

71052501 

TRAINING 

1217 

1219 

1 

075 

080 

1 

WA T ER , K 1 \/E R 

7 1052501 

TRAINING 

1239 

1245 

1 

077 

080 

1 

WAI bR,KlVbK 

71052501 

TRAINING 

1273 

1282 

1 

074 

076 

1 

WA' 1 b K , R i VER 

7 1052501 

TRAINING 

1290 

1294 

1 

08? 

085 

1 

WAIbRjKlVbK 

71052501 

TRA INING 

1291 

1293 

1 

140 

143 

1 

WA 1 EK , R 1 VE R 

71052501 

TRAINING 

32 

35 

1 

50 

57 

1 

FORAGE, PAST 

71052501 

TRAINING 

56 

67 

1 

188 

195 

1 

F UKAlib , P A S T 

71052501 

TRAINING 

10 5 

111 

I 

IRQ 

197 

1 

FORAGE , STUD 

7 10 52 5 01 

TRAINING 

115 

127 

1 

187 . 

. 193 

1 

fukAuE ,S 1 Ub 

71052501 

TRAINING 

196 

204 

1 

180 

192 

1 

F U K A (.T E , PAST 

71052501 

TRAINING 

335 

339 

1 

197 

204 

1 

F U K A G b , P A S f 

71052501 

TRAINING 

405 

419 

1 

211 

215 

1 

FORAGE , PAS t 

71052501 

TRAINING 

506 

510 

1 

073 

086 

1 

FORAGE ,PAS1. 

71052501 

TRAINING 

564 

574 

1 

205 

217 

1 

FORAGE, PAST 

71052501 

TRAINING 

650 

654 

1 

048 

056 

1 

FORAGE , PAST 

71052501 

TRAINING 

678 

681 

1 

T73 

185 

1 

FORAGE ,STUb 

7105250] 

TRAINING 

817 

825 

i 

ion 

118 

1 

FORAGE ,HAY 

71052501 

TRAINING 

853 

859 

1 

1 44 

171 

1 

FORAGE, HAY 

71052501 

TRAINING 

1143 

1152 

1 

047 

058 

■1 

FORAGE ,HAY 

71052501 

TRAINING 

1189 

1193 

■ 1 

184 

208 

1 

FORAGE , STUB 

7 1052 501 

TRAINING 

1202 

1209 

1 

184 

208 

1 

FORAGE , STUB 

71052501 

TRAINING 

1275 

1284 

1 

109 

125 

1 

FORAGE, STUB 

71052501 

TRAINING 

1449 

1452 

1 ‘ 

127 . 

142 

1 

FORAGE ,HAY 

7 1052501 

TRAINING 

661 

670 

1 

174 

181 

1 

CORN ■ 

71052501 

TRA INING 

681 

685 

1 

137 

143 

1 

CORN • 

71052501 

TRAINING 

683 

685 

1 

153 

167 

1 

CORN ' 

71052501 

TRAINING 

698 

703 

1 

141 

143 

1 

CORN 

71G52501 

TRAINING 

1250 

1257 

1 

93 

100 

1 

CORN ' 

71052501 

TRA INING 

1275 

1282 

1 

93 

100 

1 

CORN 

71052501 

TRAINING 

1403 

1412 

1 

194 

213 

1 

CORN 

71052501 

. TRAINING 

1481 

1485 

1 ■ 

033 

042 

1 

CORN f 

71052501 

TRAINING 

1487 

1512 

1 , 

049 

051 

1 

CORN 

71052501 

TRAINING 

1492 

1506 

1 

033 

042 

1 

CORN ’ 

7105? 50 1 

TRA INING 

498 

500 

1 

099 

117 

1 

SOY 

710525O1 

TRAINING 

1027 

1031 

1 

00 5 

025 . 

1 

SOY 



TARLF D3, CONTINUED 


71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052.501 

71052-501 

710-52501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052.501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 
71052501 
71.0 52 5 01 
71052501 

71052501 

71052501 

71052501 

71052501 

71052.501 

71052501 

71052501 

71052501 

71052501 

71052501 


TRAINING 

1045 

1048 

1 

TRAINING 

1191 

1193 

1 

TRAINING 

1253 

1274 

1 

TRA INING 

1473 

1479 

1 

TEST 

1 

025 

1 

TEST 

43 

051 

1 

TEST 

48 

088 

1 

TEST 

72 

83 

. 1 

TEST 

77 

107- 

1 

TEST 

88 

.110 

1 

TEST 

91 

122 

1 

TEST 

‘9 8 

134 

1 

TEST 

116 

161 

. 1 

TEST 

123 

144 

1 

TEST 

149 

158 

1 

TEST 

150 

190 

1 

TEST 

191 

209 

1 

TEST 

220 

261 

1 

TEST 

245 

286- 

1 

TEST 

246 

262 

1 

TEST 

256 

295 

1 

TEST 

297 

316 

’1 

TEST 

312 

325 

1 

TEST 

317 

361 ■ 

1 

TEST 

339 

352 

1 

TEST 

342 

366 

1 

TEST 

371 

390 

1 

TEST 

386 

403 

1 

TEST 

386 

404 

1 

TEST 

391 

402 

1 

TEST 

408 

429 

1 

TEST 

431 

455 

1 

TEST 

437 

446 

1 

TEST 

464 

491 

1 

TEST 

466 

475 

1 

TEST 

508 

526 

1 

TEST 

519 

540 

1 

TEST 

575 

578 

1 

TES'T 

575 

602 

1 

TEST 

631 

647 

1 

TEST 

635 

649 

1 

TEST 

657 

677 

1 

TEST 

785 

826 

1 

TEST 

814 

854 

1 

TEST 

858 

870 

1 

TEST 

935 

974 

1 

TEST 

1007 

1066 

1 

TEST 

1079 

1098 

1 

TEST 

1197 

1218 

1 

• TEST 

1354 

1369 

1 

TEST 

3.44 

348 

1 

TEST 

409 

411 

1 

TEST 

413 

420 

1 

TEST 

512 

514 

1 

• TEST 

294 

2 95 

1 

. TEST 

. 834 

837 

1 

TEST 

1191 

1192 

1 

' TEST 

1195 

1200 

1 

TEST 

1206 

1215 

1 

TEST 

1213‘ 

1214 

1 

TEST 

1214 

1227 

1 

TEST 

1216 

1223 

1 

TEST 

1219 

1230 

1 

TEST 

1224 

1238 

1 


124 

128 

1 

SOY 

043 

046 

1 

SOY 

129 

135 

1 

SOY ■ 

078 

'0'90 

1 

SOY 

061 

157 

1 

DEC IDUOUS 

122 

156 

1 

DECIDUOUS 

161 

179 

1 

DECIDUOUS 

113 

137 

1 

DECIDUOUS 

001 

022 

1 

DECIDUOUS 

107 

138 

1 

DECIDUOUS 

160 

178 

1 

DECIDUOUS 

073 

093 

1 

DECIDUOUS 

006 

019 

1 

DECIDUOUS 

115 

128 

1 

DECIDUOUS 

059 

103 

1 

DECIDUOUS 

212 

222 

1 

DECIDUOUS 

123 

160 

1 

DECIDUOUS 

200 

222 

1 

DECIDUOUS 

140 

155 

1 

DECIDUOUS 

177 

199 

1 

DECIDUOUS 

043 

054 

1 

DECIDUOUS 

025 

064 

1 

DECIDUOUS • 

120 

184 

1 

DECIDUOUS 

019 

040 

1 

DECIDUOUS 

166 

192 

1 

DECIDUOUS 

093 

133 

1 

DECIDUOUS 

115 

1.88 

1 

DECIDUOUS 

42 

58 

1 

DECIDUOUS 

060 

068 

1 

DECIDUOUS 

145 

172 

1 

DECIDUOUS 

160 

183 

1 

DECIDUOUS 

141 

148 

1 

DECIDUOUS 

149 

153 

1 

DECIDUOUS 

199 

210 

1 

DECIDUOUS 

001 

027 

1 

DECIDUOUS 

099 

126 

1 

DECIDUOUS 

033 

040 

1 

DECIDUOUS 

177 

181 

1 

DECIDUOUS 

161 

176 

1 

DECIDUOUS 

190 

222 

1 

DECIDUOUS 

058 

073 

1 

DECIDUOUS 

107 

12 6 

1 

DECIDUOUS 

136 

203 

1 

DECIDUOUS 

031 

051 

1 

DECIDUOUS 

180 

222 

1 

DECIDUOUS 

046 

111 

1 

DECIDUOUS 

183 

211 

1 

DECIDUOUS 

099 

130 

1 

DECIDUOUS 

094 

115 

1 

■ DECIDUOUS 

189 

222 

1 

DECIDUOUS 

07 A 

077 

1 

•CONIEEROUS 

145 

150 

1 

CON I EE ROUS 

141 

144 

] 

CONIEEROUS 

.144 

149 

1 

CONIFEROUS 

158 

159 

1 

WATER, POND 

004 

010 

1 

WATER, POND 

164 

167 

1 

WATER, RIVER 

161 

162 

• 1 

WATER, RIVER 

163 

166 

1 

WATER, RIVER 

058 

063 

1 

WATER, RIVER 

007 

007 

1 

'WATER, RIVER 

162 

165 

1 

WATER, RIVER 

083 

086 

1 

WATER, RIVER 

161 

164 

1 

WATER , RIVER 
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TABLF 03, CnMTTMUFn 


71052501 

7105250] 

7105250] 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71 Q 52 501 

71052501 

71Q52501 

71052501 

71052501 

71052501 

7] Q52501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

710 52 501 

71052501 

71052501 

7105250] 

71052501 

71052501 

71052501 

710'52501 

71052501 

7] 052501 

71Q52301 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

710 52.501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

71052501 

7105250] 

7U0 52 5pi 

71052501 

71.0 52 5 01 

71052501 

71052501 

71052501 

71052501 


TFST 

TFST 

TFST 

TFST 

TFST 

TFST 
TFST 
TFST 
TFST 
TFST 
TFST 
TFS'T 
TFST 
TFST 
TFST 
TFST 
TFST 
TFST 
TFST 
TEST 
TEST 
TFST 
TFST 
TFST 
TFST 
TFST 
TFST 
TEST 
TFST 
TEST 
TFST 
TFST 
TFST 
TEST 
TEST 
TFST 
— T-F-ST 
TFST 
TFST 
TFST 
TFST 
TFST 
TFST 
TFST 
TFST 
TEST 
TEST 
TEST 
TFST 
TFST 
TFST 
TFST 
TFST 
TFST 
TEST 
TFST 
TEST 
TFST 
TFST 
TEST 
TFST 
TFST 
TFST 
TEST 
TFST 
TFST 


1252 

1259 

1 

080 

082 

1 261 ■ 

1267 

1 

157 

159 

1'2 79 

1282 

1 

153 

15 5 

1285 

1287 

1 

149 

1 50 

1297 

1299 

1 

121 

12 6 

10 

15 

1 

023 

030 

94 

•99 

1 

185 

192 

121 

142 

1 

215 

22? 

] 35 

140 

1 

163 

174 

1-41 

144 

1 

153 

164 

164 

169 

1 

179 

184 

1 64 

173 

1 

170 

176 

182 

187 

1 

115 

135 

205 

211 

1 

00? 

007 

716 

219 

1 

128 

134 

230 

237 

1 

122 

133 

235 

241 

1 

157 

1 64 

272 

288 

280 

297 

1 

1 

21] 

124 

222 

132 

291 

296 

1‘ 

140 

149 

293 

302 

1 

213 

215 

300 

308' 

1 

153 

15 7 

330 

337 

1 

047 

051 

350 

3 57 

1 

1^5 

1 60 

382 

399 

1 

1^17 

208 

408 

417 

1 

05’1 

058 

420 

431 

1 

208 

215 

426 

429 

1 

02 7 

035 

440 

451 

1 

12? 

131 

447 

462 

1 

167 

172 

458 

467 

1 

113 

119 

463 

469 

1 

175 

182 

471 

48? 

1 

173 

184 

475 

495 

1 

056 

067 

486 

496 

1 

151 

159 

491 

501 

1 

7 6 

94 

499 

508 

1 

023 

028 

499 

509 

1 

0^2 

052 

52 9 

532 

1 

108 

120 

532 

550 

1 

088 

102 

543 

551 

1 

127 

140 

559 

565 

1 

24 

35 

•\ f C 


564 
572 
582 
619 
634 
649 
649 
689 
693 
696 
700 
725 
881 
908 
954 
969 
985 
1016 
1021 
1062 
1075 
1089 
11 23 
1131 


570 
580 
591 
628 
647 
663 
655 
7 04 
704 
703 
703 
733 
B88 
916 
968 
976 
990 
1021 
1039 
1068 
1084 
1098 
112 9 
11-^0 


irJATERjRIVFR 
WATER, RIVER 
WATER, RIVER 
WATER, RIVER 
WATER , RIVER 

FORAGE, PAST 
forage, STUB 
FORAGE , PAST 
FORAGE , PAST 
FORAGE jPAST 
FORAGE, PAST 
FORAGE ,STl)B 
FORAGE, PAST 
FORAGE, PAST 
FORAGE , PAST 
FORAGE jPAST 

FORAGE 5 PAST 
FORAGE , PAST 
FORAGE, PAST 
FORAGE , PAST 
FORAGE, STUB 
FORAGE , PAST 
FORAGE, PAST 
FORAGE , PAST 
FORAGE , PAST 
FORAGE , PAST 
FORAGE, PAST 
FORAGE ,PAST 
FORAGE, STUB 
forage ,STUP 
forage , PAST 
FORAGE, PAST 
FORAGE , PAST 
FORAGE , PAST 
. FORAGE, PAST 
FORAGE jPAST 
FORAGE, PAST 
. FORAGE, PAST 
FORAGE ,STUe 
. FORAGE, PAST 
: FORAGE, PASl 
. FORAGE, PASl 
L FORAGE, PASl 
FORAGE ,PAS; 
L FORAGE, PAS' 
L F0RAG-E,PAS' 
L .FORAGE, PAS’ 
L FORAGE, PAS’ 
1 FORAGE, PAS’ 
L FORAGE, STUf 
1 FORAGE, PAS' 
L forage, STU i 
1 FORAGE, PAS' 
1 FORAGE ,STUI 
1 FORAGE, PAS 
1 forage, PAS' 
1 FORAGE, PAS 
1 forage, PAS 
1 forage, PAS 
1 forage, STU 
3 FORAGEtPAS 
1 FORAGE, STU 
1 forage, PAS 
1 forage, PAS 
1 FORAGE, PAS 
1 FORAGE, PAS 
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TABLF 03, 
7L052501 

COMTINUED 

TEST 

1131 

1142 

1 

047 

058 

1 

FORAGE ,HAY 

71052501 

TEST 

1162 

117.5 

1 

166 

181 

1 

FORAGE , PAST 

71052501 . 

TEST 

1179 

1186 

1 

085 

102 

1 

FORAGE , PAST 

71052501 

TEST 

1196 

1206 

1 

134 

146 

1 

FORAGE, PAST 

71052501 

TEST 

1199 

1209 

1 

173 

183 

1 

FORAGE ,STUR 

71052501 

TEST 

1221 

1252 

1 

061 

069 

1 

FORAGE, STUB 

710 52 5 01 

TEST 

1252 

1274 

1 

109 

125 

1 

FORAGE ,STUB 

71-052501 ' 

TEST 

1313- 

1317 

1 

164 

180 

1 

FORAGE, PAST 

710 52 501 

TEST 

1343 

1356 

1 

050 

059 

1 

FORAGE, PAST 

71052501 

TEST 

1360 

1384 

1 

166 

180 

1 

FORAGE, PAST 

71052501 - 

TEST 

1390 

1400 • 

1 

184 

222 

1 

FORAGE , PAST 

71052501 

TEST 

1402 

1415 

1 

176 

189 

1 

FORAGE , STUB 

71052501 

TEST 

1403 

1416 - 

1 

130 

150 

1 

FORAGE ,HAY 

71052501 

TEST 

1455 

1464 

1 

118 

131 

1 

FORAGE, STUB 

71052501 

TEST 

1474 

1483 

1 

001 

010 

1 

FORAGE ,HAY 

71052501 

TEST 

1482 

148 9 

1 

149 

210 

1 

FORAGE , PAST 

71052501 

TEST 

1496 

1501 

1 

107 

126 

1 

FORAGE ,HA Y 

71052501 

TEST 

1499 

1517 

1 

161 

174 

1 

FORAGE, HAY 

71052501 

TEST 

1510 

1516 

1 

133 

148 

1 

FORAGE, HAY 

71052501 

TEST 

650 

660 

1 

174 

181 

1 

CORN 

71052501 

TEST 

698 

700 

1 

154 

169 

1 

CORN 

71052501 

TEST 

707 

714 

1 

134 

146 

1 

CORN 

71052501 

TEST 

924 

934 

1 

159 

169 

1 

CORN 

71052501 

TEST 

1212 

1248 

1 

172 

193 

1 

CORN 

71052501 

TEST 

1417 

1438 

1 

157 

167 

1 

CORN 

71052501 

TEST 

1482 

1512 

1 

052 

082 

1 

CORN 

71052501 

TEST 

1482 

1492 

1 

085 

103 

1 

CORN 

71052501 

TEST 

1503 

1512 

1 

152 

157. 

1 

CORN 

71052501 

TEST 

1510 

1513 

1 

196 

203 

1 

CORN 

71052501 

TEST 

294 

318 

1 

218 

222 

1 

SOY 

71052501 

TEST 

706 

711 

1 

190 

194 

1 

SOY 

71052501 

TEST • 

921 

926 

1 

109 

124 

1 

SOY 

71052501 

TEST 

1057 

1059 

1 

139 

157 

1 

SOY 

71052501 

TEST 

1101 

1112 

1 

158 

171 

1 

SOY 

71052501 

TEST 

1131 

1136 

1 

032 

043 

1 

SOY 

71052501 

TEST 

1137 

1155 

1 

017 

042 

1 

SOY 

71052501 

TEST 

1140 

1165 

1 

010 

012 

1 

SOY 

71052501 

TEST" 

11 44 

1161 

1 

001 

003 

1 

SOY 

710 525 01 

TEST 

1156 

1165 

1 

017 

023 

1 

SOY 

71052501 

TEST 

1175 

1179 

1 

030 

033 

1 

SOY 

71052501 

TEST 

1186 

1189 

1 

039 

042 

1 

SOY 

71052501 

TEST 

1192 

1196 

1 

013 

026 

1 

SOY 

71052501 

TEST 

1198 

1208 

1 

15 

28 

1 

SOY 

71052501 

TEST 

1215 

1218 

1 . 

046 

051 

1 

SOY 

71052501 

TEST 

1220 

1222 

1 

Oil 

030 

1 

SOY 

71052501 

TEST 

1258 

1264 

1 

017 

029 

1 

SOY 

71052501 

TEST 

1420 

1442 

1 

129 

140 

1 

SOY 

71052501 

TEST 

1420 

1429 

1 

143 

152 

1 

SOY 

71052501 

TEST 

1431 

1443 

1 

143 

150 

1 

SOY 

71052501 

TEST 

1466 

1478 

1 

033 

045 

1 

SOY 

71052501 

JEST 

1477 

1479 

1 

049 

076 

1 

SOY 

71052501 

TEST 

1508 

1514 

1 

107 

122 

1 

SOY 

71052501 

TEST 

1537 

1544 

1 

162 

180 

1 

SOY 
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TABLF D4 TEST AREAS FOR RUN 72di3?R03 


Riki 

ARFA 

FIRST 

LAST 

LIME 

FIRST 

LAST 

COL 

COVER 

MUMRRR DESIGNATION 

LINF 

LINE 

INT 

COL 

COL 

IMT 

TYPE 

72f)3?R03 

135 

268 

271 

1 

1434 

1439 

1 

CORN 

72032B03 

138 

268 

273 

1 

1470 

1486 

1 

Ct )kN 

72032803 

1 32 

272 

279 

1 

1381' 

1389 

1 

CORN 

72032803 

134- 

275 

279 

1 

1427 

1433 

1 

CORN 

72032803 

133 

282 

284 

1 

1383 

1391 

1 

CORN 

720328O3 

139 

287 

292 

1 

1460 

1465 

1 

CORN 

72032803 

141 

297 

302 

1 

1480 

148 f 

1 

CORN 

72032803 

143 

316 

319 

1 

1452 

1456 

1 

CORN 

72032803 

144 

320 

323 

1 

1453 

1457 

1 

CORN 

72032803 

148 

322 

325 

1 

1483 

1489 

1 

COf^N 

72032803 

' 1 5 

329 

336 

1 

1381 

1387 

1 

CORN 

72(!)3?863 

146 

333 

33.8 
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5 96 

598 

1 

1446 

1450 

1 

CORN 

72032803 

125 

599 

602 

1 

1758 

1762 

1 

CORN 

72032803 

24A 

601 

605 

1 

145 6 

1458 

1 

CORN 

72032803 

24B 

602 

606 

1 

1452 

1455 

1 

CORN 

72032803 

24C 

603 

607 

1 

1449 

1451 

1 

CORN 

72032803 

26 

607 

614 

1 

1453 

1461 

1 

CORN 

72032803 

127 

608 

611 

1 

1769 

1772 

1 

CORN 
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TARLR D4, 
72032B03 

CONTINUED 

42 

616 

622 

1 

147 8 

L481 

1 

CORN 

72032803 

lOlA 

616 

618 

1 

1717 

172 0 

1 

COliiN 

72032803 

100 

618 

625 

1 

1705 

1712 

1 

CORN 

72032803 

128 

618 

619 

■ 1 

1773 

1777 

1 

CORN 

72032803 

1016 

620 

622 

1 

1718 

1721 

1 

CORN 

72032803 

28 

62 7 

62 8 

1 

1471 

1476 

1 

CORN 

720 32803- 

27 

628 

630 

1 

1458 

1468 

1 

CORN 

72032803 

88 

631 

634 

1 

1655 

1658 

1 

CORN 

720328O3 

99A 

631 

632 

1 

1694 

1696 

1 

CORN 

72032803 

99B 

633 

634 

1 

1695 

1696 

1 

CORN 

720 32803 

10 3A 

633 

634 

1 

1722 

1723 

1 

CO.RN 

72032803 

93A 

• 635 

63 6 

1 

1670 

1675 

1 

CORN 

72032803 

1038 

6 35 

636 

1 

1723 

1724 

1 

CORN 

72032803 

92 

636 

640 

1 

1665 

1666 

1 

CORN 

72032803 

93B 

637 

640 

1 

1 671 

1673 

1 

CO.RN 

72032803 

103C , 

, 63t7 

63 8 

1 

1724 

1725 

1 

CORN 

7-20 32 803 

89 

’ 6 39 

643 

1 

1645 

1648 

1 

CO,RN 

72032803 

60 • 

641 

642 

. 1 

1631 

1636 

1 

CORN 

72032803 

32 

645 

650 

1 

1479 

1-482 

1 

CORN 

72032803 

58 

645 

648 

1 

1610 

1612 

1 

CORN 

72032803 

31- 

1 648 

661 

1 

1470 

1475 

1 

XORN 

72032803 

97 

649 

651 

1 

1676 

1680 

1 

CORN 

72032803 

53 

650 

656 

1 

1581 

1586 

1 

CORN 

72032803 

52A 

652 

654 

1 

1570 

1574 

1 

CORN 

72032803 

33 

652 

661 

1 

148 3 

1492 

1 

CO'RN 

72032803 ■ 

107 

653 

65 5 

1 

1714 

1719 

1 

CORN 

72032803 

52B 

655 

657 

1 

1571 

1575 

1 

CORN 

72032803 

112A 

656 

660 

1 

1732 

1738 

1 

CORN 

72032803 

51 

. 657 

662 

1 

1544 

1548 

1 

CORN 

72032803 

HOB 

659 

660 

1 

1715 

1717 

1 

CORN 

72032803 

1 lOA 

660 

661 

1 

1711 

1714 

1 

CORN 

72032803 

109A 

661 

662 

1 

1703 

1708 

1 

CORN 

72032803 

112B 

661 

664 

1 

1734 

1740 

1 

CORN 

72032803 

109B . 

662 

663 

1- 

1698 

1702 

1 

CORN 

72032803 

34 

.665 

670 

1 

1487 

1498 

1- 

CORN 

72032803 

L-86 

684 

686 

1 

1340 

1342 

1 

CORN 

72032803 

L“82 

686 

689 

1 

1320 

1323 

1 

CORN 

72032803 

L-81 

687 

690 

1 

1313 

1317 

1 

CORN 

72032803 

L-119 

688 

691 

1 

1342 

1347 

1 . 

CORN 

72032803 

L-80 

688 

690 

1 

1308 

1310 

1 

CORN 

72032803 ' 

L-77 

688 

689 

1 

1285 

1293 

i 

CO’RN 

72032803 

L~79 

690 

692 

1 

1300 

1305 

1 

CO'RN 

72032803 

L-118 

691 

693 

1 

1350 

1353 

1 

CORN 

72032803 

L-122 

693 

698 

1 

1321 

1326 

1 

CORN 

72Q32803 

L-12 3 

695 

6 99 

1 

1317 

1319 

1 

CO’RN 

72032803 

L-91 

696 

700 

1 

1301 

1306 

. 1 

CORN 

72032803 

L-74 

696 

698 , 

1 

1272 

1273 

1 

CO.'RN 

72032803 

L-90 

697 

698 

1 

1293 

1298 

1 

CORN 

72032803 

L-112 

698 

700 

1 

1349 

1356 

1 

COiRM 

72032803 

L-87 

700 

701 

1 

1273 

12 76 

1 

CORN 

72032803 

L-92 

70.1 

704 

1 

1296 

12 99 

1 

CORN 

72032803 

L-lllC 

702 

705 

1 

1352 

1355 

1 

CORN 

7-20 32 803 

L-LllB 

703 

•705 

1 

1347 

1350 

1 

CORN 

72032803 

L-68 

703 

704 

1 

1257 

1262 

1 

CO'RN 

72032803 

L-129 

705 

707 

1 

'1329 

1330 

1 

CO;RN 

72032803 

L-lllA 

705 

707 

1 

1342 

1344 

1 

CORN 

72032803 

L-67 B 

. 705 

706 

1- 

1252 

1254 

1 

CORN 

72032803 

L-51 

70 5 

707 

1 

1201 

1215 

1 

CORN 

72032803 

L-128 

706 

709 

1 

1319 

132 3 

1 

CORN 

72032803 

L-67 A 

706 

707 

1 

1247 

1250 

1 

CORN 

72032803 

L-127 

. 707 

711 

1 

131? 

1315 

1 

CORN 

72032803 

L-114 

707 

708 

1 

1365 

1369 

1 

CORN 

72032803 

L-60 

707 

709 

1 

1'231 

1234 

1 

CORN 

72032803 

L-32 

• 707 

712 

1 

1184 

1193 

1 

CORN 

72032803 

L-93 

707 

709 

1 

1296 

1301 

1 

CORN 

72032803 

L-108C 

708 

711 

1 

1352 

1355 

1 

CORN 

72032803 

L-94 

708 

710 

1 

1291 

1-294 

1 

CORN 
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72032803 

L-108B 

709 

712 

1 

1345 

1349 

1 

CORN 

72032803 

L-54 

709 

710 

1 

1200 

1220 

1 

CORN 

72032803 

L-108A 

711 

713 

1 

1335 

1347 

1 

-CORN 

72032803 

L-52 

711 

713 

1 

1207 

1217 

1 

CORN 

72032803 

L-109C 

712 

715 

1 

1352 

1355 

1 

CORN 

72032803 

L-49 

71? 

714 

1 

1200 

1204 

1 

CORN 

72032803 

L-109B 

713 

717 

1 

1347 

1349 

1 

CORN 

72032803 

L-5 

713 

718 

1 

1102 

1110 

1 

CORN 

72032803 

L-96 A 

714 

716 

1 

1285 

1289 

1 

CORN 

72032803 

L-22 

714 

716 

1 

1158 

1165 

1 

CORN 

72032803 

L-9/S B . 

715 

717 

' 1 

1281 

1283 

1 

CORN 

72032803 

L-109A 

715 

717 

1 

1342 

1346 

1 

CORN 

72032803 

L-73 

716 

717 

1 

1259 

1273 

1 

CORN 

72032803 

L-33 

717 

718 

1 

1174 

1182 

1 

CORN 

72032803 

L-107 

718 

721 

1 

1330 

1334 

1 

CORN 

72032803 

L-103A 

■718 

720 

1 

1292 

1293 

1 

CORN 

72032803 

L-37 

719 

720' 

1 

1187 

1199 

1 

CORN 

72032803 

L-2 

719 

721 

1 

1091 

1097 

1 

CORN 

72032803 

L-ino 

720 

723 

1 

1283 

1285 

1 

CORN 

72032803 

L-103B 

.72 1 

722 

1 

1294 

1295 

1 

CORN 

72032803 

L-102 

722 

725 

1 

1288 

1293 

1 

CORN 

72032803 

L-1 

72 2 

72 8 

1 

1078 

1086 

1 

CORN 

72032803 

L— 103C 

723 

724 

1 

129,6 

1297 

1 

CORN 

72032803 

L-19 

■ 723 

726 

1 

1134 

1143 

1 

CORN 

72032803 

L-38 

724 

728 

1 

1183 

1186 

L 

CORN 

72032803 

L-26 

72 4 

72 8 

1 

1162 

1172 

1. 

CORN 

72032803 

L-40 

7-25 

727 

1 

1190 

1192 

1 

CORN 

72032803 

L-17 

72 5 

728 

1 

1120 

1129 

1 

CORN 

72032803 

L-41 B 

' 726 

728 

1 

1199 

1208 

1 

CORN 

72032803 

L-41 A 

72 7 

729 

1 

1195 

1198 

1 

CORN 

72032803 

L-7 

. 730 

743 

1 

1103 

1106 

1 

CORN 

72032803 

L-44 

732 

741 

1 

1195 

1198 

1 

CORN 

72032803 

L-4 

732 

734 

1 

1078 

1083 

1 

CORN 

72032803 

L-29 

736 

738 

1 

1166 

1173 

1 

CORN 

72032803 

L-25 

740 

751 

1 

1144 

1162 

1 

CORN 

72032803 

L-30 

741 

744 

1 

1169 

1178 

1 

CORN 

72032803 

L-46 

743 

.746 

1 

1187 

1193 

1 

CORN 

72032803 

L-13 C 

744 

744 

1 

1133 

1135 ■ 

1 

CORN 

72032803 

L-13 A 

745 

745 

1 

1126 

1128 

1 

CORN 

72032803 

L-13^B" 

74 5 

745 

1 

1130 

1131 

1 

CORN 

72032803 

L-48 

748 

750 

1 

1199 

1211 

1 

CORN 

72032803 

L-8 

751 

757 

1 

1094 

1104 

1 

CORN 

72032803 

L-11 

751 

761 

1 

1131 

1138 

1 

CORN 

72032803 

L-31 

752 

757 

1 

1153 

1166 

1 

CORN 

72032803 

L-10 

755 

764 

1 

1117 

1124 

1 

CORN 

72032803 

136 

262 

266 

1 

1453 

1463 

1 

SOY 

72032803 

137A 

269 

272 

1 

1462 

1466 

1 

SOY 

72032803 

137B 

273 

276 

1 

1463 

1468 

1 

SOY 

72032803 

142 

307 

313 

1 

1450 

1457 

1 

SOY 

72032803 

150 

332 

335 

1 

1487 

1490 

1 

SOY 

72032803 

147 

338 

341 

1 

1489 

1491 

1 

. SOY 

72032803 

157 

346. 

351 

1 

1383 

1392 

1 

SOY 

72032803 

161 

347 

352 

1 

1437 

1446 

1 

SOY 

72032803 

0-156 

361 

366 

1 

1243 

1246 

1 

SOY 

72032803 

0-152 

362 

368 

1 

1221 

1225 

1 

SOY 

72032803 

163 

363 

368 

1 

1462 

1467 

1 

SOY 

72032803 

0-159 

366 

368 

1 

1254 

1256 

1 

SOY 

72032803 

173 

372 

376 

1 

1407 

1411 

1 

SOY 

72032803 

0-165 

373 

375 

1 

1270 

1272 

1 

SOY 

72032803 

0-139 

374 

376 

1 

1184 

1188 

1 

SOY 

72032803 

0-141 

377 

381 

1 

1194 

1196 

1 

SOY 

72032803 

Q-142. 

381 

383 

1 

1188 

1193 

1 

SOY 

72032803 . 

0-75 

381 

385 

1 

1122 

1124 

1 

SOY 

72032803 

0-143 

382 

383 

1 

1185 

1186 

1 

SOY 

72032803 

0-27 

386 

388 

1 

1108 

1111 

1 

SOY 

72032803 

0-81 A 

387 

390 

1 

1144 

1145 

1 

SOY 
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72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032 803 

72032803 

72032803 

72032803 

72032803 

720-32803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72Q32803 

72032803 

72Q32803 

72032803 

72032803 

720,32803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72d32803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72Q32803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72P32803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 


0-24 
0-172B 
0-172A 
0-1720 
0-81 B 
n-146 
0-17*5 
0-12 A 
0^8 ^ 

□ -12' B 
0-31 
0-182 
0-19 
0-183 
0-187 
0-2 03A 
0-190 
0-9 0 
0-203B 
0-88 
0-95 
0-35 
0-197 
0-102A 
0-62 
0-102B 
0-60 B 
0-40 
0-60 A 
0-115 
0-43 B 
0-43 A 
0-136 
0-50 
0-53 
0-46 
O-llOB 
0-llOA 
0-56 
0-57 
76A 
67B 
116A 
67A 
76B 
73A 
116B 
73B 
73C ■ 
12A 
73D - 
12B 
68 ■ 
118 
14A 
14B 
122A 
80 

122B- 

20A 

131 

20B 

19R 

123A 

19A 

123R 

25A 


388 

389 

389 

390 

391 
39} 

392 
392 
394 
394 
407 
413 
415 
417 
428 
43 5 

436 

437 
439 
443 
448 
450 
4 50 
453 
458 
458 
461 

461 

462 

463 

467 

468 
470 
476 

480 

481 

488 

489 
491 
496 
540 
542 

542 

543 
543 
54 6 
546 
549 
552 

554 

555 
555 
555 
558 
565 
568 

578 

579 
582 
585 
585 

588 

589 

589 

590 
590 
592 


389 

391 

392 

393 
396 

393 

394 
393 
401 

395 
409 
414 
417 
419 
432 
43 9 
436 
440 
442 
444 
450 
452 
452 
457 
464 

460 

461 
467 

462 
470 

46 8 
470 
472 

47 8 
483 
482 

488 

489 
492 
498 
542 

547 
545 

548 
548 

548 

549 
551 
554 
554 
557 
555_ 
561 
561 
567 
569 
581 
581 

585 
590 

586 
590 

593 

590 

594 

591 
594 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 . 

1 

1 

1 

1 

r 

1 

1 

r 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

■ 1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 


1098 
1266 
1269 
12 6? 
1146 
1201 
1254 

1054 
1030- 

1055 
1110 
1222 
10,5 8 
1222 
1204 
1284 
1212 
1146 
1286 
1115- 
]139 
10 34 
1261 
1152 
1122 
1150 
1110 
1042 
110 6 
1181 
1060 
1055 
1190 
1078 
1084 
1057 
11'37 
1132 
1110 
1113 
1729 
1665 

1746 
1660 
1731 
1697 

1747 
1699 
1701 

1450 
1703 
1446 
1645 

1751 

1448 

1449 
1758 
1715 
1760 
1465 
1753 
1456 

1451 
1755 
1444 

1752 
1458 


1098 

1268 

1272 

1265 

1147 
1202 
1256 

1055 
1034 

1056 
1114 
1223 
1062 
,1224 
1206 
1285 
121,6 

1148 
1287 
1-116 
1142 
1040 
1265 

1154 
1125 

1155 
1112 
1045 
1108 
1184 
1065 
1059 
1-195 
1082 
1089 
1057 
1141 
1136 
1114 
1116 
1733 
1671 

1748 
1664 
1733 
1699 

1749 
1701 

1703 
1454 

1704 
1448 
1652 
1754 
1453 

1453 

1761 
1720 

1762 

1459 
1754 

1460 

1454 
1758 
1450 
1754 
1465 


1 

1 

1 

1 - 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


sov 

SOY . 

SOY 

SO.Y 

SOY 

SOY 

SOY 

SOY 

SOY 

SOY 

SO.Y 

SOY 

SOY 

SOiY 

SOY 

SOY 

SOY 

SOY 

SOY 

SOY 

SOY 

SOY 

SOY 

SOY 

SOY 

SOY - 

SOY 

soy 

SOY 
SOY 
SOY 
SOY 
SOY 
SOY 
SOY 
SOY 
SOY 
SOY 
SOV 
soy 
• sov 

SOY 

SOY 

SOY 

SOY 

SOY 

SOY 

so>v 

SOY 

SOY 

SOY 

soV 

SOY 

SOY 

soV 

snV 

SOY 

SOY 

S0Y 

SQY 

SOY 

SOY 

sn’Y 

SOY- 

SOY 

SQ'Y 

SOY 
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TABLE 04, 
72032803 

CONTI MU ED 
25B 

595 

598 

1 

1459 

1466 

1 

SOY 

72032803 

126A 

598 

601 

1 

1766 

1769 

1 

SOY 

72032803 

25C 

599 

603 

1 

1461 

1467 

1 

SOY 

72032803 

84 

602 

609 

1 

1669 

1671 

1 

SOY 

72032803 

126R 

602 

605 

1 

1767 

1769 

1 

SOY 

72032803 

85 

611 

612 

1 

1661 

1667 

1 

SOY 

72032803 

41 ■ 

619 

621 

1 

148 6 

1495 

1 

SOY 

72032803 

40 

62 2 

629 

1 

1493 

1497 

1 

SOY 

72032803 

29 

625 

629 

1 

1479 

1483 

1 

SOY 

72032803 

38 

62 5 

633 

1 

1503 

1507 

1 

SOY 

72032803 

102 

628 

629 

1 

1720 

1723 

1 

SOY 

72032'803 

30 

631 

63 7 

1 

1474 

1485 

1 

SOY 

72032803 

39 

633 

637 

1 

1513 

1516 

1 

SOY 

72032803 

91 

637 

639 

1 

1657 

1661 

'1 

SOY 

72032803 

104 

638 

641 

1 

1719 

1720 

1 

SOY 

72032803 

54 

640 

643 

1 

1585 

1587 

1 

SOY 

72032803 

96 

645 

647 

1 

1674 

1678 

1 

SOY 

72032803 

95C 

64 8 

650 

1 

1668 

1670 

1 

SOY 

72032803 

106 

648 

650 

1 

1720 

1723 

1 

SOY 

72032803 

95D 

651 

652 

1 

1669 

1671 

1 

SOY 

72032803 

108 

653 

6 54 

1 

1723 

1725 

1 

SOY 

.72032803 

35 

655 

657 

r 

1496 

1497 

1 

SOY 

72032803 

50 

658 

663 

1- 

1530 

1533 

1 

SOY 

72032803 

lllA 

658 

659 

1 

1724 

1726 

1 

SOY 

72032803 

36 

659 

665 

1 

1520 

1524 

1 

SOY 

72032803 

1116 

660 

661 

1 

1725 

1727 

1 

SOY 

72032803 

lllC 

662 

663 

1 

1726 

1727 

1 

SOY 

72032803 

113 

670 

673 

1 

1729 

1733 

1 

SOY 

72032803 

L-85 

684 

687 

1 

1333 

1337 

1 

SOY 

72032803 

L-84 

685 

687 

1 

1329 

1331 

1 

SOY 

72032803 

L-120 

690 

692 

1 

1335 

1339 

1 

SOY 

72032803 

L-115 

692 

695 

1 

1358 

1361 

1 

SOY 

72032803 

L-75 

695 

696 

1 

1276 

1277 

1 

SOY 

72032803 

L-124 

697 

699 

1 

1313 

1314 

1 

SOY 

72032803 

L-B9 B 

697 

700 

1 

1286 

1290 

1 

SOY 

72032803 

L-89 A 

698 

700 

1 

1280 

1285 

1 

SOY 

72032803 

L-65 

698 

700 

1 

1250 

1251 

1 

SOY 

72032803 

L-113 

699 

702 

1 

1361 

1364 

1 

SOY 

72032803 ■ 

L-69 

699 

6.99 

1 

1264 

1268 

1 

SOY 

72032803 

L-63 

699 

702 

1 

1242 

1244 

1 

SOY 

72032803 

L-59 

701 

703 

1 

1231 

1240 

1 

SOY 

72032803 

L-53 

703 

705 

1 

1218 

1221 

1 

SOY 

72032803 

L-57 

703 

706 

1 

122 5 

1225 

1 

SOY 

72032803 

L-50 

704 

708 

1 

1197 

1197 

1 

SOY 

72032803 

L-71 

706 

707 

1 

1259 

1269 

1 

SOY 

72032803 

L“61 

706 

708 

1 

1237 

1238 

1 

SOY 

72032803 

L-126 

7 08 

712 

1 

1305 

1308 

1 

SOY 

72032803 

L-106B 

712 

717 

1 

1325 

1328 

1 

SOY 

72032803 

L-106C 

712 

716 

1 

1330 

1333 

1 

SOY 

72032803 

L-95 B 

712 

713 

1 

1299 

1302 

1 

SOY 

72032803 

L-97 

713 

714 

1 

1277 

1279 

1 

SOY 

72032803 

L-95 A 

713 

715 

1 

1292 

1298 

1 

SOY 

72032803 

L-.36 

714 

716 

1 

1191 

1197 

1 

SOY 

72032803 

L-106A 

714 

719 

1 

1316 

1322 

1 

SOY 

72032803 

L-105 

715 

721 

1 

1308 

1312 

1 

SOY 

72032803 

L-104A 

717 

719 

1 

1303 

1304 

1 

SOY 

72032803 

L-101 

719 

720 

1 

1287 

1291 

1 

SOY 

72032803 

■ L-104B 

720 

722 

1 

1304 

1306 

1 

SOY 

72032803 

L-99 A 

721 

723 

1 

1280 

1281 

1 

SOY 

72032803 

L-99 B 

■ 72 4 

725 

1 

1281 

1282 

1 

SOY 

72032803 

L-3 

724 

726 

1 

1.092 

1102 

1 

SOY 

72032803 

L-42 

730 

730 

1 

1184 

1188 

• 1 

SOY 

72032803 

L-27 B 

730 

730 

1 

1171 

1175 

1 

SOY 

72032803 

• L-27 A 

731 

731 ■ 

1 

1165 

1170 

1 

SOY 

72032803 

L-16 

-731 - 

732 

1 

1121 

1130 

1 

SOY 

72032803 

L“6 

731 

740 

1 

1092 

1097 

1 

SOY 

72032803 

L-43 

733 

739 

1 

1186 

1189 

1 

SOY 



183 


TARLF D4, 
72032803 

COMTIMUED 

L-28 

733 

733 • 

1 

1165 

1172 

1 

SOY 

72032803 

L-14 

738 

742 

1 

1126 

1133 

1 

SOY 

72032803' 

L-47 

743 

745 

1 

1197 

1206 

1 

SOY. 

72032 803 

L-12 

747 

748 

1 

1127 

1135 

1 

SOY 

72032803 

L“9 

748 

752 

1 

1109 

1121 

1 

SOY. 

72032803 
72 032 803 

0-164 

374 

376 

1 

1258 

1259 

1 

OTHER, ALFAL 

0-120 

376 

376 

1 

1178 

1129 / 

1 

OT HE R t A L FA L 

72032803 

0-78 

388 

390 

1 

1126 

1130 • 

1 

OTHER ,,ALFAL 

72032803 

0-30 

390 . 

390 

1 

1116 

1117 

1 

OT HE R , A L FA L 

72032803 

0-10 A 

391 

395 

1 

1037 

1046 

1 

OTHER » AL F AL 

72032803 
7 20 32 803 

0-14 

396 

396 

1 

1062 

1062 

1 

OT HER , A L FA L 

0-6 A 

397 

398- 

1 

1022 

1022 

1 

OTHER , ALFAL 

t L— y* y-' 

72032803 

72032803 

6-10 R 

39 8 

400 

1 

10 39 

1048 

1 

OTHER, ALFAL 

0-6 R 

399 

400 

1 

1023 

1023 

1 

OTHER ALFAL 

72032803 

0-214 

404 

405 

1. 

1279 

1284 

1 

OTHER, ALFAL 

72032803 

72032803 

6- 1 

406 

406 . 

1 

1001 

1002 

1 

OTHER ,AL F'A L 

0-213 

407 

408 

1 

1275 

1285 

1 

OTHER, ALFAL 

720 32803 

0-212 

410 

412 

1 

1278 

1286 

1 

OTHER , ALFA L 

72032803 

72032803 

0-1 83 

42 3 

42 6 

1 

1224 

1226 

1 

OT HE R ,.A L FA L 

' 0-201B 

445 

448 

1 

1282 

1286 

1 

OTHER , ALFAL 

72032803 

0-201A 

44 6 

449 

1 

1278 

1281 

1 

OTHER , A L FAL 

72032 8 03 

0-35 

467 

468 

1 

1096 

1102 

1 

OTHER, ALFAL 

72032803 

71 A 

541 

543 

1 

1676 

1680 

1 

OTHER,, ALFAL 

1 VC mC V ' vy M^y 

72032803 

71B 

544 

546 

1 

1677 

1681 

1 

OTHER , AL FAL 

72032803 

81A 

584 

585 

1 

1716 

1720 

1 

OTHER, ALFAL 

( c»v'Myt»WV2wC 

72032803 

81B 

5 86 

587 

1 

1717 

1721 

1 

OTHER , ALFAL 

72032803 

90A 

63 8 

639 

1 

1651 

1654 

1 

OTHER, ALFAL 

72032803 

90R 

640 

641 

1 

1652 

1655 

1 

OTHER ,AL FAL 

72032803 

59R 

642 

643' 

1 

162 3 

1627 

1 

OTHER t.ALFAL 

72032803 

59A 

643 

644 

1 

1619 

1622 

1 

OTHER , ALFAL 

72032803 

105 

644 

644 

1 

1719 

1722 

1 

OTHER , ALFAL 

72032803 

98 

660 

662 

1 

1674 

1674 

1 

OTHER , ALFAL 

72032803 

64A 

54 6 

'549 . 

1 

1628 

1630 

1 

OT HE R , G RASS 

1 ^vCmCc— ‘'vy -y 

72032803 

64R 

550 

554 

1 

1630 

1631 

■1 

OTHER , GRASS 

72032803 

119 

563 

566 

1 

1753 

1753 

1 

OTHER , GRASS 

72032803 

15A 

571 

572 

1 

145 0 

1451 

1 

OTHER, GRASS 

72032803 

1 5B 

57 3 

574 

1 

1451 

1452 

1 

OTHER, GRASS 

72032803 

16A 

576 

578 

1 

1452 

1453 

1 

OTHER, GRASS 

72032 80 3 

1 6R 

57 9 

580 

1 

1453 

1454 

1 

OTHER, GRASS 

72032803 

55B 

642 

644 

1 

1594 

1597 

1 

OTHER, GRASS 

72032803 

55A • 

643 

645 

1 

1590 

1593 

1 

OTHER , GRASS 

72032803 

56B 

'647 

649 

1 

1594 

1599 

1 

OTHER , GRASS 

72032803 

56A 

64 8 

651 

1 

1590 

1593 

1 

OTHER, GRASS 

72032803 

57 

650 

654 

1 

1604 

1608 

1 

OTHER, GRASS 

72032803 

0-161 

364 

366 

1 

1267 

1269 

1 

OTHER, AGRIC 

72032803 

0-162 

369 

372 

1 

] 261 

1270 

1 

OTHER , AGRIC 

72032 803 

0-77 

381 

381 

1 

1140 

1144 

1 

OTHER, AGRIC 

72032803 ’ 

0-121 

382 

385 

1 

1158 

1162 

1 

OTHER, AGRIC 

72032803 

0-144 

388 

393 

1 

1193 

1197 

1 

OTHER, AGRIC 

72032803 

0-177 

395 

397 

1 

1246 

1250 

1 

OTHERtAGRIC 

72032803 

0-2 

40 1 

405 

1 

1004 

1007 

1 

OTHER, AGRIC 

72032803 

0-84 B 

402 

■406 

1 

1127 

1131 

i 

OTHERvAGRIC 

72032R03 

0-84 A 

40 3 

407 

1 

1123 

1126 

1 

OTHERtAGRIC 

72032803 

0-179A 

408 

410 

1 

1254 

1256 

1 

OTHERtAGRIC 

72032803 

0-179B 

411 

412 

1 

] 255 

1257 

1 

OTHER t AGRIC 

72032803 

0-128 

427 

433 

1 

1191 

1199 

1 

OTHER , A6R IC 

72032803 

0-88 B 

42 8 

430 

1 

1169 

1171- 

1 

OT HER, AGRIC 

72032803 

0-88 A 

429 

431 

1 

1166 

1168 

1 

OTHERtAGRIC 

72032803 

0-73 

433 

43 5 

1 

1122 

1128 

1 

OTHERtAGRIC 

72032803 

0-129 

435 

437 

1 

1191 

1193 

1 

OTHER', AGR IC 

1 Cm vy i^y Cm V 1 VI M^ 

72032803 

0-63 

444 

444 

1 

1087 

1089 

1 

OTHER , AGRIC 

72032803 

0-133B 

453 

456 

1 

1183 

1184 ■ 

1 

OTHER ,AGRIC 

7203? 803 

0-61 

459 

460 

1 

1114 

1118 

1 

OTHER ,AGRIC 

1 Cm vy m' Cm vyvy .y 

72032803 

0-116 

473 

478 

1 

1185 

1187 

1 

OTHER ,AGR IC 

f Cm vy ay Cm 1 y vy -y 

72032803 

0-52 

482 

484 

!■ 

1072 

1074 

1 

OTHER ,AGRIC 

1 vy M*y Cm V y vy 

72032803 

0-47 

485 

487 

1 

1053 

1055 

1 

OTHER , AGRIC 



184 


TARLF DA? 

72032R03 

72032R03 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032 803' 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 ' 

72032803 

72032803 

72032803 

72032803 

72032803 

72032803 


CONTINUED 


46 

534 

536 

1 

1517 

1522 

1 

OTHER?AGRIC 

78P 

552 

554 

1 

1744 

1745 

1 

OTHER?AGRIC 

78A 

553 

555 

1 

1740 

1743 

1 

OTHERjAGRIC 

2 3B 

599 

600 

1 

1453 

1457 

1 

OTHERtAGRIC 

23A 

600 

601 

1 

1446 

1452 

1 

OTHERtAGRIC 

129A 

62 0 

621 

1 

1768 

1769 

1 

OTHER, AGRIC 

129B 

621 

621 

1 

1765 

1767 

1 

OTHERtAGRIC 

1 30 

632 

634 

1 

1778 

1781 

1 

OTHER, AGRIC 

L- 11 6 

686 

689 

1 

1356 

1359 

1 

OTHER, AGRIC 

L-83 

686 

688 

1 

1326 

1326 

1 

OTHERjAGRIC 

L~117 

689 

690 

1 

1350 

1352 

1 

OTHER, AGRIC 

L-78 

692 

694 

1 

1289 

1294 

1 

OTHER, AGRIC 

L-76 

693 

695 

1 

1280 

1284 

1 

OTHER, AGRIC 

L-125 

697 

699 

1 

1310 

1310 

1 

OTHERjAGRIC 

L-88 

699 

701 

1 

1278 

1278 

1 

OTHERjAGRIC 

L-64 

699 

701 

1 

1246 

1248 

1 

OTHERjAGRIC 

L-70 

701 

703 

1 

1265 

1269 

1 

OTHERjAGRIC 

L-58 

702 

704 

1 

1228 

1228 

1 

OTHERjAGRIC 

L-66 B 

703 

703 

1 

1253 

1255 

1 

OTHERjAGRIC 

L-66 A 

704 

704 

1 

1247 

1251 

1 

OTHERjAGRIC 

L-62 

705 

707 

1 

1242 

1242 

1 

OTHERjAGRIC 

L-2 0 

70.6 

717 

1 

1143 

1150 

1 

OTHERjAGRIC 

L-56 

707 

709 

1 

122 6 

1226 

1 

OTHERjAGRIC 

L-110 

707 

708 

1 

1334 

1339 

1 

OTHERjAGRIC 

L-55 

708 

710 

1 

122 3 

1223 

1 

OTHERjAGRIC 

L-98 

70 8 

711 

1 

1276 

1278 

1 

OTHERjAGRIC 

L-7? 

709 

711 

1 

1271 

1273 

1 

OTHERjAGRIC 

L-35 

715 

717 

1 

1185 

1188 

1 

OTHERjAGRIC 

L-23 

719 

721 

1 

1160 

1169 . 

1 

OTHERjAGRIC 

L-i 8 

719 

719 

1 

1138 

1141 

1 

OTHERjAGRIC 

L-21 

■ 721 

724 

1 

1148 

1155 

1 

OTHERjAGRIC 

L-34 

72 1 

721 

1 

1176 

1180 

1 

OTHERjAGRIC 

L-39 

722 

723 

1 

1189 

1200 

1 

OTHERjAGRIC 

L-24 

726 

729 

1 

1154 

1158 

1 

OTHERjAGRIC 

L-45 

741 

741 

1 

1188 

1191 

1 

OTHERjAGRIC 

0-82 A 

383 

385 

1 

1150 

1154 

1 

OTHER, PAST 

0-82 B 

386 

388 

1 

1151 

1155 

1 

OTHER, PAST 

0-82 C 

389 

389 

1 

1152 

1156 

1 

OTHERjPAST 

0-196 

452 

455 

1 

1245 

1251 

1 

OTHER, PAST 

69A 

551 

554 

1 

1663 

1665 

1 

OTHERjPAST 

69B 

555 

557 

1 

1664 

1666 

1 

OTHERjPAST 

69C 

558 

560 

1 

1666 

1668 

1 

OTHERjPAST 

86A 

611 

613 

1 

1677 

1681 

1 

OTHERjPAST 

1 40 

295 

297 

1 

1462 

1468 

1 

OTHER jWHEAT 

155 

319 

32 5 

1 

138'7 

1393 

1 

OTHERjWHEAT 

149 

32 7 

330 

1 

1485 

1489 

1 

OTHERjWHEAT 

171 

386 

389 

1 

1398 

1406 

1 

OTHERjWHEAT 

0-2 00 

44 8 

451 

1 

1267 

1273 

1 

OTHERjWHEAT 

87A 

629 

630 

1 

1662 

1665 

1 

OTHERjWHEAT 

87B 

631 

63 3 

1 

1663 

1665 

1 

OTHERjWHEAT 

95A 

642 

644 

1 

1 666 

1668 

1 

OTHERjWHEAT 

958 

64 5 

647 

1 

1667 

1669 

1 

OTHERjWHEAT 

OGLE-WTR 

338 

344 

1 

1105 

1109 

1 

WOODS 

OGLE-WTR 

343 

348 

1 

1121 

1130 

1 

WOODS 

OGLE-WTR 

352 

358 

' 1 

1131 

1135 

i 

WOODS 

SYCAMORE 

351 

361 

1 

1595 

1605 

1 

TOWN 

DEKALB 

445 

455 

1 

1545 

1555 

1 

TOWN 

DEKALB-M 

455 

465. 

1 

1525 

1535 

1 

TOWN 

HINCKLEY 

.643 

650 

1 

1786 

1795 

1 

TOWN 




6 



